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FOREWORD 


Competition has pushed aircraft gas turbine engine 
designs to ever higher levels of performance. The resulting 
higher pressures and temperatures have led to a decrease in 
engine durability. In hot section components - the combus- 
tor and turbine - the more hostile environments have acceler- 
ated the damage and wear out of parts, thus, dramatically 
increasing maintenance costs. While recent technology 
advancs have been aimed primarily at improving performance, 
minimal efforts have been made to assure required durability. 

The activities of the NASA Turbine Engine Hot Section 
Technology (HOST) Project are directed toward durability 
needs, as defined by industry, and toward a more balanced 
approach to engine design. The HOST efforts will improve the 
understanding of aerothermal environments and material 
structural responses by focused analytical and experimental 
research activities. The Project supports both contract/ 
grant and Lewis in-house research activities. The overall 
approach is to assess existing analysis methods for strengths 
and deficiencies, conduct supporting analytical and experi- 
mental research to rectify the deficiencies, incorporate 
state-of-the-art improvements into the analysis methods, and 
finally verify the improvements by systematic test programs. 

To provide representatives from government, industry, 
and universities with the latest findings and progress toward 
improved aircraft turbine engine durability, a two-day 
workshop was held in October 1983. This publication contains 
the papers presented at the workshop. 


i i i 


Daniel E. Sokol owski 
Manager, HOST Project 
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TURBINE ENGINE HOT SECTION TECHNOLOGY (HOST) PROJECT 


Daniel E. Sokolowski 
and 

C. Robert Ensign 


National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, OH 44135 


The Hot Section Technology (HOST) Project is a NASA-sponsored endeavor to 
improve the durability of advanced gas turbine engines for commercial and 
military aircraft. Through improvements in the analytical models and life 
prediction systems, designs for future hot section components — the combustor 
and turbine -- will be more accurately analyzed and will incorporate features 
required for longer life in the more hostile operating environment of high 
performance engines. 

Started in 1981, the project has activities presently planned through 1989 
with an estimated total cost of $50-million. While the focused research 
activities are necessarily analytical in nature, significant experimental 
testing is required for benchmark quality assessments as -well as model 
verifications. The efforts are being conducted in-house at the NASA Lewis 
Research Center, under contracts at major domestic turbine engine 
manufacturers, and under grants to qualified universities. Currently, the 
contract and grant funding equals about 60-percent of the total. 

At NASA Lewis Research Center, the HOST project is the major element in a 
project management and coordination office and serves as the focal point for 
advocacy, funding, technical coordination, and information exchange. This 
workshop serves as the primary vehicle for this last function; that is, to 
disseminate information and elicit the exchange of ideas among participants. 

The activities of the HOST project are accomplished within six disciplines as 
shown in Figure 1. Management of the project's efforts is delegated to six 
Subproject Managers (Figure 2) in the following areas: Structural Analysis, 
Fatigue & Fracture, Surface Protection, Combustion, Turbine Heat Transfer, and 
High Temperature Instrumentation. Structural Analysis includes research into 
thermal mechanical load models, component geometry specific models, and 3-D 
inelastic analysis methods development. Fatigue and Fracture includes 
constitutive model development for both isotropic and anisotropic materials, 
including single crystal and directionally solidified forms. It also includes 
research in life prediction methods for creep-fatigue interactions, and 
elastoplastic crack propagation. The Surface Protection research includes 
studies of corrosion and oxidation phenomena, environmental mechanics models, 
and metallic and thermal barrier coating analysis method developments. The 



Combustion work includes aerothermal model assessment and development, 
dilution jet modeling, high pressure flame radiation/heat flux testing, and 
development of a thermal structural cyclic test facility. The Turbine Heat 
Transfer area is studying 2-D and 3-0 flow and heat transfer on airfoil 
external boundaries, emphasizing boundary layer transition and viscous 
modeling. It also investigates coolant passage heat transfer, including 
midchord jet impingement cooling and rotational passage effects. 
Instrumentation is being developed to obtain high temperature, benchmark 
quality data to develop and verify the analysis methods. These include flow 
sensors (LDV), heat flux sensors (thin film), strain sensors (1800°F static 
thin film), gas temperature sensors (frequency compensated), and hot section 
optical viewing systems. 

Schedules for activities under each discipline are shown in Figures 3-8. 
Because of the nature of the problem of durability, interdisciplinary 
cooperation is often required and the synergistic effects have proven to be 
beneficial. For example, researchers in combustors, high temperature 
instrumentation, and structural analysis are jointly investigating combustor 
liner behavior, modeling, and life prediction. 

Workshop publications and most contractor final reports carry the label "For 
Early Domestic Dissemination" (FEDD) to protect national interests and, thus, 
are available only to qualified U.S. citizens. While several contractor final 
reports have been published recently, they represent initial phases of 
multi-phased work. Thus, aside from this annual workshop report, only a few 
reports are available at this time. 
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Figure 3. - Structural analysis. 
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Figure 5. - Surface protection. 




Figure 6. - Combustion. 
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Figure 7. - Turbine heat transfer. 




Figure 8. - Instrumentation. 
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HOST INSTRUMENTATION R&D PROGRAM 


OVERVIEW 


D. R. Englund 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


The HOST Instrumentation R&D program is focused on two categories of 
instrumentation. One category is that required to characterize the 
environment imposed on the hot section components of turbine engines. 

This category includes instruments for measuring gas flow, gas tem- 
perature, and heat flux. The second category is that for measuring 
the effect of the environment on the hot section components. This 
category includes strain measuring instruments and an optical system 
for viewing the interior of an operating combustor to detect cracks, 
buckling, carbon buildup, etc. 

The HOST Instrumentation R&D program was formulated to concentrate on 
the critical measurements that could not be made with commercially 
available instruments or with instruments already under development via 
NASA- or DOD-funded efforts, or in IR&D programs. Over the past year 
we have not added any new measurements to the program. 

The HOST Instrumentation R&D program schedule is included in the accom- 
panying figures. The program schedule shows all HOST-funded efforts 
plus non-HOST-funded efforts that were initiated prior to HOST and which 
have HOST-related goals. Each line represents a separate effort, either 
contract, grant, or in-house. In two of these efforts, the turbine 
blade/vane static strain gage and the dynamic gas temperature measurement 
system, follow-on work is shown beyond the original contract. 

The heaviest resources are concentrated on the measurements of strain and 
gas flow because these measurements are critical to the success of the 
HOST program and because the HOST requirements differ from the current 
state of the art by a considerable margin. Follow-up and complementary 
efforts not shown in the schedule are being planned for the strain mea- 
surement area. 
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HOST INSTRUMENTATION R&D PROGRAM 


GENERAL GOALS: 

• DEVELOP INSTRUMENTATION FOR CHARACTERIZING THE ENVIRONMENT 

AROUND TURBINE ENGINE COMPONENTS 

• DEVELOP INSTRUMENTATION FOR CHARACTERIZING THE EFFECT OF THE 

ENVIRONMENT ON THE TURBINE ENGINE COMPONENTS 

CS-83-2793 
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I 1 HOST FUNDED V////A NON-HOST FUNDED 


GOAL 


HIGH-TEMPERATURE SYSTEM FOR VIEWING 
COMBUSTOR INTERIORS DURING OPERATION 


DYNAMIC GAS TEMPERATURE MEASUREMENT 
SYSTEM WITH 1 kHz RESPONSE 


TURBINE BLADE/VANE STATIC STRAIN GAGE 

BURNER LINER STATIC STRAIN MEASUREMENT 
SYSTEM 

BURNER LINER TOTAL HEAT FLUX SENSOR 

TURBINE BLADE/VANE TOTAL HEAT FLUX SENSOR 

LASER ANEMOMETER SYSTEM FOR HIGH- 
PRESSURE HIGH-TEMPERATURE FLOWS 


OPTIMIZATION OF LASER ANEMOMETER 
SYSTEMS 


CS-83-2795 
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HOT SECTION VIEWING SYSTEM 


William W. Morey 

United Technologies Research Center 
East Hartford, Connecticut 06108 

The objective of the hot section viewing program is to develop a prototype 
optical system for viewing the interior of a gas turbine combustor during high 
temperature, high pressure operation in order to produce a visual record of some 
causes of premature hot section failures. The program began by identifying and 
analyzing system designs that would provide clearest images while being able to 
survive the hostile environment inside the combustion chamber. Different illumi- 
nation methods and computer techniques for image enhancement and analysis were 
examined during a preliminary test phase. In the final phase of the program, which 
will be discussed in this presentation, the prototype system was designed and 
fabricated and is currently being tested on a high pressure combustor rig. 

The combustor viewing system consists of a 12.7 mm diameter viewing probe 
with an actuator unit that mounts onto a rig or engine, an optical interface unit 
to transfer images to video and film cameras and provide illumination light, and 
two control chassis for remote operation of the probe and the optics interface 
unit. Probe actuation consists of a ±180 rotation from a downstream view and 
insertion radially along the combustor axis with a depth range of 7.6 cm. In order 
to cover a wide range of views, two viewing probes were constructed. One probe has 
a wide view, the other a narrow view. Both probes have their direction of view 
angled off the probe axis so that the view can be scanned on rotation. In addi- 
tion, each probe has a choice of two lenses. The wide field probe has 90 and 60° 
f ield-of-view lenses, and the narrow field probe has 35 and 13 f ield-of-view 
lenses. For the wide field probe the view is oriented at 45 to the probe axis by 
bending the fiber optic image bundle. A cross section view of this probe is shown 
in figure 1. Since the lens focal lengths are much longer for a narrow f ield-of- 
view, the narrow field probe had to use a mirror to turn the view from the probe 
axis. Figure 2 shows a cross section of the narrow field probe. 

Both probes are made of copper and have water-cooled walls. Four baffles in 
the probe walls guide cooling water to the distal end of the probe. Hot, soot- 
ladened combustor gas is kept from the probe interior and from the fiber and lens 
exposed surfaces with a flow of purge gas. The gas flows down the center of the 
probes and exits around the illuminating fibers and out of the lens aperture or 
across the mirror in the case of the narrow field probe and into the combustor. 

The probes were designed to handle the high heat fluxes expected in a primary 
combustion zone. A finite element analysis indicated that the probe tip tempera- 
tures were well within material limits and that cooling water at the probe tip 
should be below the boiling point. Measured temperatures were about 50 C higher 
than calculated. 

The image carrying fiber element is 3 mm in diameter and consists of a fused 
bundle of about 75,000, 10 p.m diameter fibers. Each fiber corresponds to a picture 
element. The fused fiber bundle is 33 cm long and transfers the image from the 
viewing lens at the distal end of the probe to the opposite end of the probe where 
it is butt-coupled to a 3 m long flexible fiber bundle. The flexible bundle 
transfers the image to the optics interface unit for recording by video or film 
cameras. One millimeter diameter plastic clad fused quartz fibers, shown in the 



figures, were used for transferring the illumination light source. Two fibers were 
used in each probe. There are also two thermocouples mounted in each probe tip. 

The viewing probe mounted in its actuator along with the remote control 
chassis is shown in figure 3. The purge gas line connects to the upright fitting 
on the left-hand side of the probe and the four bent tubes with connectors are for 
the probe water cooling. A separate line is used for cooling the actuator case and 
mounting plate. The remote control unit has a probe tip temperature indicator 
along with water, gas, and temperature interlocks. Probe depth is indicated with 
an LED bargraph meter. Rotation is indicated in degrees with a digital panel 
meter. The rotational speed can be varied. Another view of the probe and actuator 
with covers removed is shown in figure 4. The flexible fiber bundle can be seen 
leaving in the upper left-hand corner of the figure, and the actuator motors for 
depth and rotation can be easily identified along with their position indicating 
potentiometers . 

The optics interface unit is shown in figure 5. The 3 m long fiber optic 
image bundle from the probe comes into a filter wheel in the upper left-hand corner 
of the figure. A series of eight different optical filters can be used such as 
neutral density, color, and interference filters to vary exposure and spectral 
•range. A pneumatically-controlled slide positions a mirror that directs the image 
from the filter wheel to a video camera. When the mirror is positioned away from 
the line-of-view, the image can be seen with a film camera. A mercury arc lamp, 
its power supply, and the two illuminating fibers are on the right-hand side of the 
board. A second pneumatically-operated mirror directs the collimated and focused 
arc lamp light to the fiber ends. With the mirror out of position, a laser light 
beam can be focused into the fibers. There is also a control chassis that permits 
remote operation of the filter wheel, mirrors, and cameras with position indi- 
cators . 

Tests are currently underway with the viewing system in a high pressure 
combustor rig. A modified JT-12 combustor liner is used in the rig as illustrated 
in figure 6. Three viewing ports on the side of the rig are available for mounting 
the probe. Two of the ports are shown in the figure with the probe attached to 
port 1. The depth range of the probe (7.6 cm) in relation to the combustor is also 
shown in the figure. Both film and video data of the combustor liner surfaces will 
be recorded along with probe temperature and operational conditions. 


12 



WIDE FIELD PROBE 


VIEWING LENS 



Figure 1 

NARROW FIELD PROBE 



Figure 2 


13 



wwM 




VIEWING PROBE AND ACTUATOR 





Page intentionally left blank 


Dynamic Gas Temperature Measurement System* 

Denny L. Elmore, Woodrow W. Robinson and William B. Watkins 
Pratt and Whitney Aircraft Engineering Division 


SUMMARY 


The objective of this effort was to develop a gas temperature measurement 
system with compensated frequency response of 1 KHz and capability to operate in the 
exhaust of a gas turbine combustor (Figure 1). Results of the initial portions of 
this effort were reported in the first Hot Section Technology Conference (Reference 
lj. Further progress in this development effort is summarized in this presen- 
tation. Environmental guidelines for this measurement are presented, followed by 
a preliminary design of the selected measurement method. Transient thermal 
conduction effects were identified as important; a preliminary finite-element 
conduction model quantified the errors expected by neglecting conduction. A 
compensation method was developed to account for effects of conduction and convec- 
tion. This method was verified in analog electrical simulations, and used to 
compensate dynamic temperature data from a laboratory combustor and a gas turbine 
engine. Detailed data compensations are presented. Analysis of error sources in 
the method were done to derive confidence levels for the compensated data. 

Environmental Guidelines and Sensor Selection 


The sensor design and environmental guidelines for this effort are listed in 
Figure 2. Environmental parameters are representive of a modern gas turbine engine 
combustor exhaust, and sensor life, accuracy, spatial resolution and vibration 
capability are nominal goals for experimentation. 

The sensor approach chosen for development is shown in Figure 3. The thermo- 
couple employs two wires of different sizes to obtain data necessary for evaluation 
of the time constant for each wire. If two thermocouples are positioned in close 
proximity such that both are exposed to the same instantaneous temperature and 
velocity, the difference in their thermal responses will be governed by their 
relative diameters. These responses can then be used to obtain time constants for 
compensating the thermocouples. The unique feature of these wire thermocouples is 
a beadless laser butt-welded element, which allows the geometry to be modelled as 
a cylinder in crossflow. 

A detailed review of the platinum/ rhodium alloys tensile strengths, melting 
temperatures, emf outputs, thermal conductivities, stresses to rupture, and specific 
heats was made. Tensile and stress-rupture values are higher for increasing rhodium 
content, indicating the best thermocouple should have a high rhodium content. The 
Type B (platinum - 6% rhodium/platinum - 30% rhodium) was selected over the other 
commercially available thermocouple materials based on its higher emf output, 
availability, and known fabrication characteristics. Detailed structural analysis 
revealed that allowable yield stresses contrained design length-to-diameter ratios 
to less than 6.5 for the support wires and less than 15.5 for the thermocouple 
elements. 


* Work performed under contract NAS3-23154 



Probe Thermal Analysis and Compensation Method 


A preliminary thermal model of the thermocouple probe wires was generated to 
evaluate the effects of radiation losses and end conduction losses. Governing 
equations and the model nodal breakup are shown in Figures 4 and 5 respectively. 

The physical thermocouple junction, the small wire, and the larger support wires 
were simulated by a finite difference model. The gas temperatures were assumed to 
be periodic fluctuations. The worst case gas temperature was assumed to be a mean 
gas temperature of 1400K (2520°R) with fluctuations of + 500K (+ 900°R). The 
more realistic case defined in the data acquisition analysis, +_ 65K (+ 117°R) was 
also evaluated. These cases were evaluated at frequencies of 20, 100, and 1000 Hz. 

Results of the transient simulations (shown in Figure 6) indicate that radia- 
tion heat losses are not significant (less than 10K) , but the conduction losses are 
too large (maximum of 67K) to be ignored and because of structural requirements 
previously discussed, cannot be reduced by making the wires longer. This discovery 
required that analysis of the measured temperature data be performed with a second 
order equation containing both conduction and convection terms, rather than the 
simplified first order equation containing only convection terms. The second order 
energy equation is a linear equation in time and space, whereas the first order 
equation is only time dependent. 

A compensation method which accounts for both conduction and convection 
effects was developed. A second finite difference model using nine nodes (Figure 
10) was used for the remainder of the calculations in this effort. The method is 
outlined in Figure 7 and explained below. 

1. The theoretical transfer functions between the 76 }Jim (3 mil) t/c signal 6 ^ 
and the gas stream signal a and the 254 /im (10 mil) t/c signal 0 ? and 
the gas stream are computed from the finite difference solution of the 
differential equations for various values of an aerodynamic parameter F, 
at a number of discrete frequencies falling between the corner frequencies 
of the two t/c's. These data are then used to compute the theoretical 
transfer function H t (f) between the 254jUm (10 mil) t/c and the 76pm 
(3 mil) t/c ( Bd 6 w for the corresponding values of f and frequency. 

These curves (^/pj) will be used to determine the insitu value of F 
from the measured transfer function of B^f 0-^. The process is as follows: 

a. The following parameters are input or already stored in the computer. 
For the thermocouple wire - L, 1, D, d,Pw, Kw, Cp w , and cx w . For 
the gas stream - p g , K g , c P g >7 g >p g and Pr g . 

b. The average or mean conditions for the test data for the following 
variables are entered into the computer: 

T = mean gas temperature 

P = mean gas pressure 

F/A = fuel air ratio 

f -i — f = frequencies of f at which transfer functions will be 

1 x v j. j n 

evaluated 

Mn = Mach No. 

c. The program computes an estimated value of F based on the test 
conditions. 
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The program then computes t q , the transfer function between the wire 
thermocouples and the gas stream, for the 76/im (0.003 in.) and the 
254 jjLm (0.010 in.) thermocouple from 0.5 F to 1.5Tin steps of 0.1 F at 

frequencies f, , f which are user selected to fall in between 

estimated valaes of the corner frequencies of the two t/c's. The 
equations are evaluated until steady state conditions are reached. The 
computer code determines the sampling interval for each frequency 
evaluated to ensure mathematical stability of the finite element model 
and minimize computation time. The normalized ratio of the magnitude 
of the temperature fluctuation in the wire to the temperature fluctu- 
ation of the gas stream ( = 0 /a ) at frequency f is determined 

by locating the maximum peak amplitude after the model has iterated to 
steady-state conditions. The phase shift ( 7] ) of the temperature 
fluctuation in the wire is determined by locating the time at which 
the crossed zero going positive at the beginning of the period in 
which the model reached steady-state conditions. 


e. The data from (d) are then used to compute the theoretical transfer 
function 0 2r / 6 i n f rom 0,5 ^ t0 1,5 ^ at frequencies of fj through f x . 

2, Thermocouple test data are digitized into the Fourier system computer, 
typically 32 to 120 records each of the 76jUm t/c dynamic signal and the 
254jum t/c dynamic and dc signals. Each record contains 2048 samples of 
data. These data are then converted from millivolts to temperature 
utilizing NBS calibration curve coefficients for type B thermocouples. 

The 254/zm dc channel is utilized as the mean for both dynamic data 
channels in converting the nonlinear t/c mv singals to linearized temp- 
erature. These data records are then saved for recall for additional 
processing or plotting. 


3. An ensemble averaged FFT (Fast Fourier Transform) transfer function 

analysis is then performed on x number of records of the dynamic records 
to yield the measured value of 02 n / $in as a f unc ti° n °f frequency. The 
transfer function is computed as tne Fft crosspower divided by the FFT 
autopower of the 76 jim signal: 

H (f) = 01n 02n ^ 02n 

d in Q In 6 In 


where * = complex conjugate multiplication. 

In conjunction with the computation of the measured transfer function, the 
coherence function is computed and utilized to assess the quality of the 
measurement. For the 2048 time sample data record lengths utilized 1024 
line FFT's are produced. For the typical sampling rate of 4096 Hz, the 
FFT analyses yield spectral information dc to 2048 Hz in 2 Hz intervals. 

A standard Hewlett Packard windowing function (P301) is utilized prior 
to computation of the FFT's. This window is characterized by excellent 
spectral amplitude accuracy ( +0.1%). Side lobe suppression is >- 70 db 
at +_ spectral lines and the effective noise bandwidth is 3.4 spectral 
1 i nes . 

4. Each measured value of 0 2n / $ in at frequencies f n = fj - f are used in 
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conjunction with the theoretical curves of 0 2 n /0 ln vs F to determine a 
measured value of T (the program interpolates^Between the 0.1 r increments 
computed in (1) above). The arithmetic average obtained for each 
frequency is taken as the insitu measured value. 


5. Using the measured insitu average value of r obtained in (4), £ g , the 
normalized transfer function (gain 0, /a and phase 77 , ) of the 76 //m 
thermocouple with respect to the gas 'stream temperature! s then computed 
at all frequencies from the 1st spectral line of the FFT spectrum to the 
Nyquist frequency of the FFT for each discrete frequency contained in the 
FFT. This is typically from 2 Hz to 2048 Hz in 2 Hz increments. This 
spectrum is then used to compensate the 76fj.m t/c data as follows. 


6. To compute the compensated ensemble averaged power spectral density 

function, the ensemble averaged auto power spectrum of the 76jum (3 mil) 
t/c obtained in (3) above is divided by the auto power of its compensation 
spectrum: 
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where 


* = 


complex conjugate multiplication. 


Scaling factors for effective noise bandwidth and FFT Fourier symmetry are 
applied. 


7. To compute the compensated instantaneous time waveform, an FFT spectrum is 
made on a specific data record. This spectrum is then divided by the 
compensation spectrum. The compensated spectrum is then inverse Fourier 
transformed to yield the compensated instantaneous time waveform. The 
software contains information on specific techniques employed to prevent 
time waveform distortions associated with the inverse Fourier transform. 

A threshold, in relative db, is applied to the frequency spectrum of the 
data signal prior to division by the compensation spectrum to prevent 
errors where the signal to noise ratio is too low. 
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Sensor Test Program 

The test program summarized in the following section was based on three test 
series consisting of (1) System Shakedown and Compensation Verification Lab Tests, 
(2) Laboratory Burner Tests, and (3) Full Scale F100 Engine Tests (Figure 8). 

These tests allowed for a step by step checkout and optimization of the various 
components of the system while allowing for the experimental evaluation and sub- 
stantiation of sensor guidelines and development effort goals. 
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Equations were derived for the electrical analog equivalent, passive RC net- 
work of the nine node finite difference thermocouple model for use in lab eval- 
uations of the dual t/c approach. Two breadboard nine node RC element network 
models, one for the 76jJLm (3 mil) and one for the 254 pim (10 mil) thermocouple were 
fabricated. The circuits were modeled for the F100 probe geometry operating at the 
following conditions: 

T t = 1200°K (1700°F) 

Pj = 19.7 atm (290 psia) 

1^ = 0.233 

F/A = 0.02. 

Design values of resistance and capacitance for these simulators are shown in 
Figure 11 . Potentiometers were used for the resistive components and values 
were set by measurement with a digital ohmmeter (+ 1%). Values of capacitive 
components were combinations of standard value capacitors (+ 10%) and were not 
measured. Schedule constraints precluded obtaining more precise capacitance 
values. 

The purpose of this test was to determine the compensation spectrum (gain and 
phase as a function of frequency) of the 76 jum (3 mil) analog circuit (and the 
254 jum analog circuit) from the measured transfer function taken between the out- 
puts of the analog circuits of the 76 jlm (3 mil) t/c and the 254 jUm (10 mil) 
circuit (Figure 9 ) utilizing the HOST software and compare it with the known 
frequency response spectrum (measured single input - single output using conven- 
tional FFT techniques). Wideband (1250 Hz BW) random noise was input in these 
measurements. Because of the lack of precision in the analog circuit components, 
thermocouple length and diameter values input to the HOST program compensation 
software were varied until best agreement between the HOST program results and the 
known compensation spectra were obtained. The conventionally measured spectra 
were generated by measuring the transfer function between the input and output of 
the analog circuit using the 256 ensemble averaged FFT transfer function: 

H(f ) = Gainz. phase (f)= G in^ f ^ G out ^ 

where * = complex conjugate multiplication 

The best match for the 76 fJLm (3 mil) t/c was obtained when it was modeled with 
an element diameter of 74.9jUm (2.95 mils) and an element length of 711.2/im (28.00 
mils). The analog circuits were designed to model the t/c as 76.2/im (3.00 mils) 
diameter with a length of 609.6 fjim (24.00 mils). Gain match was within + 15% and 
phase match within + 4.2 to - 0.8 deg. over the frequency range of 4Hz to 1000Hz. 

Similarly, for the 254 /im (0.010 in.) thermocouple case, the following para- 
meters produced the best match of compensation spectra: the diameter was 260.4/im 

(10.25 mils) and the length was 1587. 5 jLlm (62.5 mils) as opposed to design values 
of diameter equal to 254jUm (10 mils) and length equal to 1587 . 5 jjtm (62,5 mils). 
Gain disagreement for this case was less than + 5% out to 100Hz. The compensation 
program developed under this contract uses the 254/im (10 mil) data only out to 
about 40Hz. The agreement between the two methods is satisfactory despite lack of 
precision in the capacitance values in the analog model. Experiments performed 
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with the analog circuits verified that the finite difference models of the thermo- 
couples work as required. Minor discrepancies, due to lack of precision in elec- 
trical component values, could easily be corrected, given adequate time and 
resources . 

Subscale combustor probes were fabricated based on the detailed designs shown 
in Figure 14, which also shows the completed probe. The subscale combustor probe 
element junctions were inspected to determine the junction's actual dimensions 
and weld uniformity. Each thermocouple junction was inspected under a Nikon Meta- 
phot microscope at 100X or 200X magnification* Figure 15 . Photomicrographs of 
each junction and a graduated scale were made to measure the junction diameter and 
uniformity of the wire cross section. Two angular orientations of each wire were 
photographed 90 degrees apart, A number of different readings of each photograph 
were obtained to define the mean value and precision of each sample. Measurement 
uncertainty of the technique is + 2 am (9xl0 -5 in.) for the 76|im elements and 9/Xm 
(4x10"^) for the 254 fim elements. Effects of these uncertainties on temperature 
uncertainty were later determined in sensitivity analyses using the completed data 
reduction method. 

The probe was installed and tested in the subscale combustor rig (Figure 16 ), 
The probe and a platinum tube used for measuring gas stream total pressure were 
mounted downstream of the combustor exit plane. The mounting system was retractable 
so the probe was not in the gas stream except when data were recorded. 

Dynamic temperature data and ambient temperature data were recorded on FM tape. 
The output of each Preston amplifier was double recorded on the FM tape recorder 
(0.5V rms and 1.0V rms full-scale record level channels). The signal levels were 
monitored on an oscilloscope and the fixed gain setting (s) of the Preston 
amplifier(s) were set as necessary to maintain optimum signal levels on the 0.5V 
rms record channel. Gain settings were correlated with test points to maintain 
calibration. Data were recorded at ten steady-state combustor test conditions 
ranging from the minimum conditions obtainable with the combustor facility. Average 
gas temperatures varied from 1048K (1887°R) to 1797K (3235 Q R). The gas stream 
Mach number and total pressure was relatively constant at 0.23 and 1.0 atms 9 
respectively. 

Figure 17 are uncompensated and compensated dynamic temperature waveforms 
for the 3 mil thermocouple and test point 10. The uncompensated waveform has max- 
imum fluctuations of approximately 400°F peak to peak with an rms value of §4°F. 

The compensated waveform has peak to peak fluctuations of approximately 120u°F and 
an rms value of 281°F. Figure 18 are uncompensated and compensated power spectral 
density plots for 120 record averages. Signal to noise ratios for compensated 

signal to compensated ambient noise are >18db up to 500Hz* > 13db up to 1000Hz and 

>10db out to 2000Hz. A noticeable break in the compensated PSD can be seen at 
1200Hz*, this results from use of an anti-aliasing low-pass filter, and for this 

reason, signal to noise ratio from 1200 to 2000Hz is actually greater than lOdb. 

Note that a combustor resonance is indicated by the peak at 145Hz. 

The full scale test were conducted as planned in an F100 jet engine which 
was undergoing testing at P&WA’s Government Products Division (GPD) test facility. 
Test conditions for the data which follows were: P = 158 psia, T = 1790°F, 

M = 0.355. Figure 19 shows the uncompensated and compensated dynamic temperature 
waveforms. The uncompensated waveform has peak to peak characteristics of about 
400°F and rms value of 73°F, The compensated waveform has peak to peak charac- 
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teri sties of greater than 2000°F, and rtns value of 393°F, Figure 20 are power 
spectral density plots of compensated and uncompensated 120 record average data. 
Signal to noise ratio is >26db to 500Hz,>20db to 1000Hz and >14db to 2000Hz. The 
signal is basically wide band random , with no evidence of resonance. 

Measurement Uncertainty 

The overall compensation uncertainty was obtained by room-sum-squaring errors 
associated with finite-windowing effects in the Fast Fourier Transform, thermo- 
couple element diameter measurement, temperature waveform ratios, and data acquisi- 
tion/reproduce noise. These errors are shown in Figure 21 for two input waveforms. 
The 27°F peak to peak//Hz wideband random waveform and corresponding errors are 
representative of the signals observed experimentally. 




L 

l 

D 

d 


Pw 

C pw 

Pg 


Mg 

7g 

^ln 

^2n 

a n 

A 


"w 


Pw C pw 


LIST OF SYMBOLS 

thermal diffusivity of the wire 


length of thermocouple support wire 
one half of length of smaller thermocouple wire 
diameter of thermocouple support wire 
diameter of smaller thermocouple wire 
density of thermocouple wire 

thermal conductivity of thermocouple wire 
specific heat of ther thermocouple wire 
density of the gas stream 


specific heat of the gas stream 


Prandtl number of gas stream 


velocity of the gas stream 

viscosity of the gas stream 

ratio of specific heats of gas stream 

amplitude of the 75 jum (0.003 in.) thermocouple at frequency f n 

amplitude of the 254 jjtm (0.010 in.) thermocouple at frequency f n 

amplitude of the gas temperature at frequency f n 

phase shift of the gas temperature with respect to arbitrary time t Q 
at frequency f n 
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^2n 

\n 

^2n 
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phase shift of 76fj.m (0.003 in.) thermocouple with respect to gas 
temperature at frequency f 

phase shift of 254 //m( 0.010 in.) thermocouple with respect to gas 
temperature at frequency f 


phase shift of 76/j.m (0.003 in.) thermocouple at frequency f with 
respect to arbitrary time t Q 

phase shift of 254 jum (0.010 in.) thermocouple at frequency f with 
respect to arbitrary time t Q n 


0.48k p 1 / 3 ^ 2 
irg q 


Pa \ 1/2 


(« 


P c 

'w pw 


Aerodynamic parameter 


J = Theoretical transfer function (gain and phase) of the wire thermo- 
couple x with respect to the gas temperature at frequency f 

H f (f) = Theoretical transfer function between large and small thermocouples 

H e (f ) = Experimental transfer function between large and small thermocouple 
signals 

Kg = thermal conductivity of the gas stream 
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DYNAMIC GAS TEMPERATURE 
MEASUREMENT SYSTEM 


Objective - Develop a temperature measurement system with: 

• Compensated frequency response of 1 KHz 

• Operation at jet engine combustor exit 

Talk Outline - • Environmental guidelines and sensor selection 

• Transient conduction effects and compensation method 

• Test program 

- Compensation checkout 

- Lab burner 

- Engine 

• Error analysis 

• Summary FIGURE 1 AV238751 

SENSOR DESIGN AND 
ENVIRONMENT GUIDELINES 


Geometry: Annular combustor, 2 cm < H < 8 cm 

Temperature: T ~ 1400° K (2060° F); T' ~ 500° K (900° F) - 

Frequency response: 1 kHz 

Pressure: 10 < P < 20 ATM 

Flow: v ~ 150 m/s; v' ~ 50 m/s - 

Gas composition: Fuel (nominal jet A) and air 

Sensor life: 5 hr minimum 

Accuracy: Temperature uncertainty < 5% for f < 200 Hz 

Temperature uncertainty 10% for 200 Hz < f < 1 kHz 
Spatial resolution: D < 0.5 cm 
Vibration: lOg FIGURE 2 
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• Ratio of T/C responses determines hg 


0.003 in. Type B— v 0.080 in. 


■ Laser butt weld 
(no bead) 

0.010 in. Type B 
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Alumina 
probe body 
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0.0875 in. 
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Sec B-B 


FIGURE 3 


FIRST-ORDER AND SECOND-ORDER 
MODELS COMPARED 


First-order model 

Rate of energy change = Q con vection + Qradiation 
— = - (Tg - T) + - 4 — - (l 8 4 - T 4 ) 

dt PyyCpyyD PwQpV/Q 

Second-order model 

Rate Of energy = Qconvection "F Qconduction "F Qradiation 
dT 4h /T „ , 5 2 T , 4 ere /T . 4 

dt Pw^pwD Pw c pwD 


FIGURE 4 
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FIGURE 5 AV238776 


MAXIMUM DEVIATION IN PREDICTED 
JUNCTION TEMPERATURE FROM A 1ST 


ORDER SYSTEM 


60 

50 

40 

30 

20 

Deviation in 10 

temperature 0 

- °K -10 

-20 
-30 
-40 
-50 
-60 



Conduction effects must be included to meet accuracy goals 


FIGURE 


AV238785 


27 



COMPENSATION METHOD 


• Compute theoretical response (76/nm and 254jum T/C’s) vs heat transfer 

coefficient (finite element conduction effects included) over 
frequency range 

• Measure (data) response of 76^m and 254/um T/C’s over frequency range 

(using FFT techniques) 

« Determine actual heat transfer coefficient from computed and 
measured response 

• Generate theoretical response of 76jum T/C for actual heat transfer 

coefficient for frequency range 

• Compensate 76/um T/C data in frequency domain 

• Inverse Fourier transform to time domain 

AV250399 

FIGURE 7 

DYNAMIC GAS TEMPERATURE 
MEASUREMENT SYSTEM 

Test program 


• Compensation method checkout 

• Lab burner 

• Engine 

FIGURE 8 AV260427 
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SIMULATION OF DUAL WIRE T/C’S 



analyzer 


• RC networks modeled as F100 probe 
1200° K (1700°F) 

19.7 ATMS (290 psia) 

Mn = 0.2231 

F/A= 0.02 FIGURE 9 an^t 

FINITE DIFFERENCE THERMAL MODEL 

Used in computer compensation program 
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RC ANALOG NETWORK SIMULATION 

Full scale test conditions 



Thermocouple element diameter 
Where: 



76fim 

254^m 

Ro 

7592 fi 

46.2 Kfi 

Ri 

1739 fi 

11.8 Kfi 

r 2 

217.4 fi 

5.66 Kfi 

r 3 

3796 fi 

23.11 Kfi 

r 4 

1049 fi 

12.447 Kfi 

r 5 

543.5 fi 

8.515 Kfi 

r 6 

1087 fi 

17.031 Kfi 

Co 

0.5 /ufd 

0.5 jifd 

C , 

1.0 M fd 

1.0 /ifd 

c 2 

2.45 M fd 

4.34 fihi 

c 3 

78 /tfd 

7.67 fifd 


Note: 


Values listed 
correspond to FI 00 
probe operating at: 

T t = 1200°K (1700°F) 

P T = 19.7 atm (290 psia) 
M n = 0.2331 
F/A = 0.02 


AV2S7666 630902 2493B 


FIGURE 11 


SIMULATION OF DUAL WIRE TC’S 


76 (3 mil) compensation spectrum - measured single 

input/single output 


Gain (db) 


7 


2 



0 (deg) 


FIGURE 12 
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76 fim (3 mil) compensation spectrum - measured 
dual wire 



4 > (deg) 
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SENSOR GEOMETRY 








Photomicrograph 
of 0.003 in. 
element (200X) 


Photomicrograph 
of 0.010 in. 
element (100 x) 


Measured Measured 

diameter = 0.0029 diameter = 0.0099 


FIGURE 15 
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SUBSCALE COMBUSTOR RIG TEST DATA 



F/A = 0.025 (est’d) 

• Ambient - probe retracted/rig 


running 


FIGURE 16 
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SUBSCALE COMBUSTOR 
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FIGURE 18 












DYNAMIC GAS TEMPERATURE 
MEASUREMENT SYSTEM 


Overall accuracy - instantaneous time waveform 
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Overall accuracy - averaged frequency spectrum 
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FIGURE 21 
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HIGH TEMPERATURE STATIC STRAIN 
SENSOR DEVELOPMENT PROGRAM 


Charles Hulse and Richard Bailey 
United Technologies Research Center 
and 

Howard Grant 

Pratt & Whitney Aircraft Group 

The purpose of this program is to develop electrical resistance strain gages 
useful for static strain measurements on nickel or cobalt superalloy parts inside a 
gas turbine engine on a test stand. Measurements of this type are of great impor- 
tance in meeting the goals of the Host Program because, without reliable knowledge 
of the stresses and strains which exist in specific components, it will be difficult 
to fully appreciate where improvements in design and materials can be implemented. 

The first part of the effort has consisted of a strain gage alloy development program 
which is to be followed by an investigation of complete strain gage systems which 
will use the best of the alloys developed together with other system improvements. 

The specific goal for the complete system is to make measurements to 2,000ye 
with error of only 1 10% over a 50 hour period. In addition to simple survival 
and stability, attaining a low thermal coefficient to resistivity, of order 100 
ppm/K or less, is also a major goal. This need results from the presently unavoid- 
able uncertainty in measurements of the exact temperatures in the turbine. The 
size and thickness requirements to avoid aerodynamic effects suggests the use of 
the sputtering technique as the best system fabrication approach. The results from 
the first year of this effort resulted in the identification of an FeCrAl alloy 
and the Pd-Cr alloy systems as the basis for further alloy modifications and 
development. Alloy candidates are evaluated and compared using a grading system 
consisting of the product of the following factors with their total weight potential 
given in parenthesis: Repeatability (20), Oxidation (18), Resistivity (16), Thermal 

Coefficient of Resistivity (14), Elastic Range (12), Differential Thermal Expansion 
(10) and Miscellaneous Judgements (10) . 

Although the final strain gage system fabrication procedure will probably use 
the sputtering technique, this is too expensive and too slow a process for alloy 
development work. A drop-casting technique using repeatedly arc-melted buttons 
was therefore developed. The cast rods produced are ground and polished using 
slotted plates to finally produce long, thin strips of material about 17cm long 
suitable for subsequent testing. Other techniques for producing the strain gage 
elements, such as wire drawing and melt spinning, are also to be examined in this 
program. 


Work done under NASA contract NAS3-23722. 



The specially constructed thermal cycling apparatus developed earlier in this 
work to make resistivity measurements has been further developed to make more 
accurate measurements. This system employs a split tube metal heater which can 
be cycled or held at constant temperature under program control. The test sample 
is positioned axially in the center of this tube with Platinum leads for voltage 
measurements and thermocouple wires attached by spot welding. Special provisions 
have been added to better accomplish the faster (250K/min) heating and cooling 
rate testing. Additional computer program development has also been accomplished 
to improve both the data measurement and the data reduction process. 

The FeCrAl alloy system studies have defined a ternary area whose edges define 
compositions whose average thermal coefficients of resistance are zero between room 
temperature and 1250K. Compositions along these edges provide significant improve- 
ments over the previously known best alloy (Kanthal A-l). Measurements of the 
resistances to oxidation of these alloys with different amounts of Hf, Y, Sc, Co 
and additional amounts of A1 and Cr have shown that, although all show protective 
behavior, none are significantly better than the base FeCrAl alloy. 

The Pd-Cr alloys develop a thin coherent, continuous coating of Cr 203 . A 
variety of other alloys containing oxide formers, especially the rare earth metals, 
are being examined to find compositions with even better resistance to oxidation 
and lower thermal coefficients of resistance. 
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GENERAL RULES CONCERNING p AND a 


Matthieson: p = p 0 + p y 

dp 

Dellinger: p a = Constant a = ~ 

Linde: A p oc V X 2 

Vs valence 
X= alloying element 

Norbury: pa(N - Nx) 2 

N = atomic number 
X= alloying element 


Figure 1 


HIGH-SPEED THERMAL CYCLE/ 
RESISTIVITY MEASUREMENT APPARATUS 


Springs 


Sample 
Thermocouple 


Clamps 



Figure 2 


Adjustable 
furnace mount 
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KANTHAL A-1 AND FeCrAI MOD #3 
AFTER 2 HRS. AT 1153°K 
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Figure 3 
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DEMONSTRATION TEST OF BURNER LINER STRAIN 
MEASUREMENT SYSTEMS: INTERIM RESULTS* 


Karl A. Stetson 

United Technologies Research Center 
East Hartford, Connecticut 


Howard P . Grant 
Pratt and Whitney Aircraft 
East Hartford, Connecticut 


Abstract 

Work is in progress to demonstrate two techniques for static strain measure- 
ments on a jet engine burner liner. Measurements are being made with a set of 
resistance strain gages made from Kanthal A-l wire and via heterodyne speckle 
photogrammetry . The background of the program is presented along with current 
results . 

Introduction 

A previous program (NAS3-22126) has resulted in identifying two potentially 
applicable methods for static strain measurement of the thermal deformations of jet 
engine burner liners. The first employs resistance strain gages made from Kanthal 
A-l wire and the second is a form of photogrammetry based upon laser speckle photo- 
graphs ( speck legrams) of the surface under study. The photographs are made by 
illuminating the object with laser light and making photographic recordings on high 
resolution plates through a telecentric lens system. The coherence of laser light 
creates high contrast speckles in the image that move as if attached to the object 
surface. Strain is revealed as image magnification by photogrammetric comparison 
of pairs of specklegrams . The accuracy of this comparison is greatly enhanced by 
use of a heterodyne interferometer. 

Both methods have features that are attractive as well as problematic. The 
wire gages provide continuous output in a conventional format and, because they are 
bonded to the actual object surface, are not influenced by rigid-body tilts and 
translations. On the other hand they are subject to failure, they can exhibit large 
and variable values of apparent strain, they have temperature limitations, they 
only measure strain at specific locations, and they can influence the mechanical 
behavior of the object. By comparison, speckle photogrammetry offers virtually no 
temperature restriction, allows continuous geometrical mapping of strain distri- 
butions, has an apparent strain equal to thermal expansion, and is noncontacting. 
Nonetheless, it only samples the object at discrete times, it can be disturbed by 
rigid-body motions, it requires optical access to the object surface, and it can be 
influenced by thermal inhomogeneities in the gas through which the object is 
viewed . 

*This work is being performed under contract NAS3-23690. 
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Program Description 


The current program combines both methods in a prototype study of combustor 
behavior. A JT12D combustor (see Fig. 1) has been modified for operation in the 
temperature range desired (Up to 950°K for the Kanthal gages and up to 1140°K for 
the photogrammetry) . A section has been cut from the burner and instrumented with 
ten gages and seven thermocouples as shown in Fig. 2. A hole has been provided so 
that the louver lip may be photographed and included in the photogrammetric strain 
measurement. The gaged section has been put through temperature calibration tests 
to predetermine the apparent strains of the gages. Reference gages have also been 
tested on a constant strain bar to determine the gage factor as a function of 
temperature. The section of the burner was then welded back into place and the 
burner installed in the test stand. 

A window was provided in the test cell as shown in Fig. 3 through which the 
burner is illuminated and photographed. Figure 4 shows the layout of the optical 
components. A pulsed laser provides illumination of about 1 msec duration with 
adequate energy to expose the plates. The system is operated remotely from the 
control room. 

During the actual test runs , specklegram data is recorded in the form of 
photographic exposures at designated times, while strain gage data is recorded 
continuously as electrical signals. To maintain correlation between the speckle 
patterns of successive photographs, the burner is heated isothermally by increasing 
the inlet air temperature until spontaneous ignition is possible. After ignition, 
various cycles can be simulated by variation of the fuel flow. After testing is 
complete, both strain gage and specklegram data are processed to obtain strain 
information. 

Processing for the strain gage data consists primarily of numerical correc- 
tion for apparent strain and gage factors. Processing for the specklegram data 
consists of developing the photographic plates and evaluating them on an inter- 
ferometric comparator. The number of locations on the photographs at which two- 
dimensional strain data may be obtained is limited by the time available for data 
processing. Strains are to be measured, therefore, only at the ten gage locations 
and at ten additional locations including the louver lip. Thermocouple readings 
will be used to correct for thermal expansion and provide thermomechanical strain. 
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FIG. 3 
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DEVELOPMENT OF HEAT FLUX SENSORS FOR TURBINE AIRFOILS AND COMBUSTOR LINERS 

William H. Atkinson 
Pratt & Whitney Aircraft Group 
Engineering Division 

This paper briefly describes the work performed under two contracts: NAS3-23529, 
"Turbine Blade and Vane Heat Flux Sensor Development"; and a pre-HOST effort, 
NAS3-22133, "Advanced High Temperature Heat Flux Sensor Development". The objective 
of each of these contracts is to develop heat flux sensors for use on gas turbine 
engine hot section components. 

The design of durable turbine airfoils that use a minimum amount of cooling air 
requires knowledge of the heat loads on the airfoils during engine operation. 
Measurement of these heat loads will permit the verification or modification of the 
analytical models used in the design process and will improve the ability to predict 
and confirm the thermal performance of turbine airfoil designs. Heat flux sensors 
for turbine blades and vanes must be compatible with the cast nickel -base and 
cobalt-base materials used in their fabrication and will need to operate in a 
hostile environment with regard to temperature, pressure and thermal cycling. There 
is also a need to miniaturize the sensors to obtain measurements without perturbing 
the heat flows that are to be measured. 

At the start of the turbine blade and vane program, a literature search and a 
survey of commercially available sensors were conducted to determine the current 
state of the art for heat flux sensors. These investigations revealed that there 
were no existing sensors which met the geometrical and environmental requirements of 
the blade and vane heat flux measurements. An effort was then initiated to develop 
feasible sensor designs for turbine blade and vane applications. Several concepts 
were chosen for further consideration: 

o embedded thermocouple sensors using the airfoil wall as a thermal barrier, 

o Gardon gage sensors fabricated into the airfoil wall, 

o several types of transient sensors. 

Sensor concepts that were only suitable for transient testing were subsequently 
dropped from consideration because of the complexity of introducing sufficiently 
rapid transient conditions into the hot side heat transfer coefficient or into the 
airfoil coolant during actual engine testing. 

Figure 1 is a schematic of the embedded thermocouple sensor. This type of one 
dimensional steady state sensor determines the heat flux by measuring the 
temperature drop across a thermal barrier. In this case, the airfoil wall itself 
acts as the thermal barrier. The Gardon gage sensor is shown schematically in figure 
2. This type of steady state sensor determines the heat flux from the temperature 
rise due to radial conduction from the insulated section of the hot side of the 
airfoil wall. Both these sensors have been optimized for steady state conditions but 
may also be used as transient heat flux sensors. 

A series of thermoelectric tests were conducted on blade and vane alloys to 
determine the suitability of these materials as thermoelectric elements. The use of 
the airfoil material in the thermoelectric circuit gives a direct differential 
output and reduces the number of wires required, minimizing the potential for 
temperature perturbations caused by the presence of the wires. All of the 
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superalloys tested were found to produce positive thermoelectric outputs with 
temperature when paired with platinum. These thermoelectric outputs were stable 
throughout four thermal cycles and several hours of thermal soaking. Alumel was 

selected as the preferred thermoelectric material because of the relatively high 
voltage output, availability and low cost. Millivolt outputs versus temperature of 
the various superalloys when paired with Alumel are presented in figure 3. 

The sensor designs were analyzed for thermal performance using a three 
dimensional finite difference heat transfer program (TCAL) developed at Pratt & 
Whitney Aircraft. This program determined the thermal perturbations caused by the 
presence of the sensor and the theoretical outputs to be expected from the sensor. 
This information was used to calculate the anticipated error in the heat flux 
indicated by the sensor compared with the heat flux that existed in the unperturbed 
airfoil. 

Both the embedded thermocouple and Gardon gage sensor designs were fabricated 
into first stage turbine blade halves. For the embedded thermocouple sensors, 
grooves were eloxed into the inner and outer airfoil surfaces. The thermocouples 
were installed in the grooves using a wedge wire technique, and the surface of the 
installation was smoothed off. Figures 4 and 5 show the interior and exterior walls 
of a typical installation. The blade halves were then welded together and the 
assembled blade is shown in figure 6. 

For the Gardon gage sensors, a cavity was eloxed into the interior wall of the 
blade half and a groove was eloxed to route the leadwire away from the sensor. To 
minimize the number of grooves required, the Gardon gage sensor used a three wire 
swaged lead. The thermocouple junctions were made in the cavity by resistance 
welding, and the leadwire was installed in the groove using a wedge wire technique. 
The cavity was then filled with a ceramic cement to provide oxidation and mechanical 
protection for the thermocouple wires and to restore the aerodynamic integrity of 
the internal blade wall. Figure 7 shows the completed Gardon gage sensor on the 
interior blade surface. 

A calibration fixture, shown in figure 8, was designed and fabricated for the 
calibration of a irfo i 1 -mounted heat flux sensors. This fixture is mounted below a 
quartz lamp bank heat source, shown in figure 9. The fixture allows positioning of 
the airfoil so that the surface of the heat flux sensor is normal to the incident 
radiation. The surface of the airfoil is coated with material having a known and 
stable absorptance and emittance. Hie incident radiation is measured with a 
reference Hy-Cal asymptotic calorimeter. This same test apparatus is used for 
thermal cycle tests ana thermal soak tests as well as the calibration tests. 

Data from the calibration of a typical embedded thermocouple sensor are shown in 
figures 10 and 11. Figure 11 shows that sensor sensitivity (sensor output/heat flux 
transmitted) varies with temperature. This would be expected because of variation in 
the thermal conductivity of the blade material and the thermoelectric output with 
temperature. If these factors are analytically accounted for, the sensor output can 
be normalized to an 1150 K sensor reference temperature. The normalized data for 
this sensor are presented in figures 12 and 13. Similarly, data for a typical Gardon 
gage sensor are presented in figures 14 and 15. These data normalized to an 1150 K 
sensor reference temperature are shown in figures 16 and 17. 

Several sensors were subjected to thermal cycling tests, in which the sensors 
were rapidly heated to near maximum temperature and heat flux, held at this 
condition for one minute, cooled rapidly to near ambient temperature and held there 


46 



for three minutes. The normal test program was 50 thermal cycles. Calibration data 
for a typical sensor before and after the thermal cycle test are shown in figure 18. 
Thermal soak tests were also conducted on several sensors in which the sensor was 
held at near maximum temperature and heat flux for a period of 10 hours. Calibration 
data for a typical sensor before and after the thermal soak test are shown in figure 
19. 


The sensors developed under the first phase of the turbine blade and vane heat 
flux sensor program have been shown to be capable of measuring the heat loads on 
turbine airfoils with an accuracy of +10 percent, and have withstood thermal cycling 
and thermal soak conditions that are to be expected to be representative of a gas 
turbine engine environment. A second phase of this program will demonstrate a wide 
variety of heat transfer measurement methods on a simple test piece in an 
atmospheric combustor rig. This work, just recently initiated, will consist of 
running a cylinder in cross flow behind a laboratory combustor rig. The first 
portion of this work will be to characterize the combustor exit as fully as possible 
with respect to temperature, pressure and velocity, including the spatial and 
temporal variations of these quantities. Various types of heat flux sensors and 
other heat flux measurement instrumentation will be mounted in cylindrical test 
sections and tested behind the combustor. This instrumentation will include both 
transient and steady state sensors as well as instrumentation for measuring dynamic 
gas and metal temperatures. The results of the various measurement techniques will 
be correlated and compared with the analytically predicted heat flux levels. 

Under a pre-HOST NASA contract (Contract NAS3-221 33, "Advanced High Temperature 
Heat Flux Sensor Development", from the NASA LeRC Instrumentation R&D Branch) 
described in Reference 1, high temperature steady state heat flux sensors were 
designed and fabricated for use in combustor liners. There were three sensor designs 
generated: 

o the laminated sensor, 

o the Gardon gage sensor, 

o the embedded thermocouple sensor. 

Figure 20 shows a schematic representation of each of these sensors. The sensors 
were subjected to laboratory testing consisting of calibrations, thermal cycle tests 
and thermal soak tests. The test results on these sensors were presented at the 1982 
HOST Contractor's Workshop. 

In a continuation of that contractual effort, a combustor liner has been 
instrumented with eight heat flux sensors and is ready to be tested in a high 
pressure combustor facility. There are three Gardon gage sensors, three embedded 
thermocouple sensors and two laminated sensors mounted in the combustor liner. An 
overall view of the cold side of the liner is shown in figure 21, and a close up 
view of the cold side of two of the sensors is shown in figure 22. The hot side of 
one of the sensors, presented in figure 23, shows the resistance weld around the 
periphery of the sensor that was used to fasten the sensor into the liner wall. 
After the sensors were installed into the liner, a functional check was performed to 
verify that they were still operational. The liner will be run as a part of the test 
sequence on NASA Contract NAS3-22392, "Broad Specification Fuels Program". 

These high pressure combustor tests will also include two radiometer probes for 
measuring the radiative portion of the total heat load. These two radiometers, a 
porous plug radiometer probe developed at Pratt & Whitney Aircraft and a 
commercially available Medtherm radiometer probe, are shown in figure 24. The porous 
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plug radiometer probe has a field of view nearly 180 degrees while the Medtherm 
radiometer has a 50 degree field of view. During the test, the two radiometers will 

be mounted on opposite walls facing each other as shown in figure 25. In this 
configuration, both probes will view the same volume and, hence, the signals should 

correlate. It is anticipated that the two probe locations will be interchanged 
during the test program to eliminate any rig bias. It is hoped that the data from 
these tests will help clarify some of the wide angle and narrow angle radiometer 
data taken in other test programs. 

As a result of the two contracts for heat flux sensor development, high 
temperature steady state heat flux sensors have been fabricated for use on turbine 
blades and vanes and on combustor liners. The use of these sensors in upcoming test 
programs will provide detailed information on the heat loads on hot section 
components, and will represent a significant advance in the modeling and prediction 
of the heat transfer characteristics of gas turbine engines. 
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Figure 1 Schematic of Embedded 
Thermocouple Sensor 


Figure 2 Schematic of Gardon Gage 
Sensor 



Figure 3 Plot of Thermoelectric Output of Superalloys versus Alumel 
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Figure 4 


View of Interior Wall of 
Embedded Thermocouple Sensor 


Figure 5 View of Exterior Wall of 
Embedded Thermocouple Sensor 




Figure 7 View of Interior Wall of 
Gardon Gage Sensor 


Figure 6 View of Assembled Embedded 

Thermocouple Instrumented Blade 





Figure 8 


View of Calibration Fixture 


Figure 9 View of Calibration Setup 
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Figure 10 Plot of Output versus Heat 
Flux for Embedded 
Thermocouple Sensor 


Figure 11 Plot of Sensitivity versus 

Temperature for Embedded 
Thermocouple Sensor 



Figure 12 Plot of Normalized Output Figure 13 Plot of Normal ized Sensitivity 

versus Heat Flux for versus Temperature for Embedded 

Embedded Thermocouple Thermocouple Sensor 

Sensor 



Figure 14 Plot of Output versus Heat Figure 15 Plot of Sensitivity versus 
Flux for Gardon Gage Sensor Temperature for Gardon 

Gage Sensor 
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Figure 16 Plot of Normalized Output 
versus Heat Flux for 
Gardon Gage Sensor 



Figure 18 Plot of Data before and 
after Thermal Cycle Test 
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Figure 17 Plot of Normalized 
Sensitivity versus 
Temperature for Gardon 
Gage Sensor 
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Figure 19 Plot of Data before and 
after Thermal Soak Test 
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Figure 21 View of Cold Side of 
Instrumented Combustor 
Liner 


Close up View of Cold Side 
of Two Combustor Liner 
Heat Flux Sensors 








Figure 23 Close up View of Hot Side Figure 24 View of Porous Plug and 

of Combustor Liner Heat Medtherm Radiometer Probes 

Flux Sensor 
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Figure 25 Schematic Representation of Radiometer Installation 
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LASER ANEMOMETRY FOR HOT SECTION APPLICATIONS 


Richard G. Seasholtz, Lawrence G. Oberle, and Donald H. Weikle 
National Aeronautics and Space Administration 
Lewis Research Center 


INTRODUCTION 

The objective of this in-house program is to develop laser anemometers 
(LA's) for use in the study of the hot section components of turbomachinery. 
Specifically, laser anemometers are being developed for use in the 50.8-cm 
(20-in.) diameter warm turbine and hi gh-pressure turbine (HPT) facilities at 
Lewis. This paper presents a brief review of the status of the program along 
with some preliminary data taken in an open-jet burner. 

A stepwise approach is being followed which builds on the experience 
gained in the development of laser anemometer systems for cold-flow turbo- 
machinery facilities (refs. 1 and 2). The procedure being followed for the 
development of the hot section laser anemometer systems is 

(1) establish measurement objectives for each facility — parameters to be 
measured and accuracies 

(2) develop an analytical model for the complete LA systems 

(3) perform necessary experiments to provide data not otherwise available 
e.g., seed properties, surface relectivity, and window contamination 

(4) experimentally verify the model using the open-jet burner facility 

(5) optimize LA system design to give desired accuracies in minimum run 
time. 


MEASUREMENT OBJECTIVES 

The measurement objectives for both the warm turbine and the high pres- 
sure turbine facilities are to measure the following parameters with the given 
accuracies: mean velocity, 1 percent; flow angle, 1 ; turbulence intensity, 

10 percent; and turbulence scale, 20 percent. Also, it is desirable that 
these measurements be made to within 2 mm of the end walls. Finally, because 
of the high cost of operating these facilities, the LA systems must be de- 
signed to take these data in a minimum of experimental run time. 


OPEN JET BURNER 

An open-jet burner facility is being used as a test bed for the develop- 
ment and testing of all components of the laser anemometer system — optics, 
traversing hardware, seed generators, and windows. Furthermore, the extensive 
software needed for control and data acquisition will be written and tested in 
this facility before the system is installed in the turbine facilities, thus 
minimizing the amount of expensive debugging in the turbine facilities. 

The open-jet burner facility consists of a general purpose laboratory 
combustor (fig. 1) and associated controls. The mixing and flow in this com- 
bustor closely approximate those of many conventional gas turbine combustors 
and the critical features of the combustion chemistry are reproduced. The 


57 



fuel (Jet A) is introduced at the front of the combustor by a pressure atom- 
izing fuel nozzle. The combustion is stabilized in the forward section of the 
combustor by developing a strong, swirl stabilized recirculation zone. The 
hot combustion gases are cooled by dilution jets which are located at the back 
of the combustor. The primary and secondary jets have independent air sup- 
plies to permit adjustment of the velocity and temperature profiles at the 
exit. The ignitor, a flame monitor, and instrumentation ports are mounted in 
an instrumentation ring located between the primary and secondary sections. 
Seed particles are introduced through a tube mounted in one of the instru- 
mentation ports. The 50.8-mm (2-in.) diameter exhaust can be operated over 
a temperature range of 1000 to 1600 K (1300° to 2500° F) at Mach numbers of 
0.2 to 0.8. 


OPTICAL CONFIGURATION 

Several types of optical designs are being considered. The primary 
design is a conventional single component dual-beam fringe-type configuration 
as illustrated in figure 2. The other types of laser anemometers under con- 
sideration are 

(1) a two-spot time-of-fl ight LA, which offers the ability to measure 
smaller particles than a fringe LA but suffers from low data rates 
in highly turbulent flow. 

(2) a new design time-of-fl ight system developed under Grant NAS 3-2 at 
Case Western Reserve University (ref. 3). This configuration using 
four elliptical spots should give data rates similiar to a fringe 
system, but should be significantly more sensitive. 

To illustrate the effect that various optical system parameters have on 
the performance of a fringe-type LA we can write the expression for the lower 
bound of the variance of the velocity measurement for a single particle 
(ref. 4) as 


2 8 / he \ r o v w o m 

mv ff l/2 \nXoflj p 2 + ig/^ 2 

where r 0 represents the background light flux; v, the velocity; 2w 0 , the 
diameter of the probe volume; P, the laser power; N, the number of fringes; 
h, Planck's constant; c, the velocity of light; x, the wavelength of the laser 
light; n, the quantum efficiency of the photodetector; n, the solid angle of 
the light collection optics; and cr, the light-scattering cross section of the 
particle. Equation (1) shows that to reduce the error for a given measurement 
it is desirable to 

(1) have particles with a large light-scattering cross section (yet small 
and light enough to follow the flow, (e.g. about 1 micrometer diameter 
for most turbomachinery flows) 

(2) have small f number optics (i.e. larger n)that captures as much 
of the light scattered from the particle as possible 

(3) have a small probe volume with many fringes 

(4) minimize the amount of stray background light (flare) 

Note that the fringe spacing for a given velocity is constrained by the 
frequency responses of the photodetector and the signal processor (typically 
about 200 MHz). For example, with a 500-m/sec flow velocity resulting in a 
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150-MHz Doppler frequency, the fringe spacing is v/f = 3.3 micrometers. The 
minimum probe volume diameter 2w 0 is then set by the number of cycles re- 
quired by the signal processor. If 10 cycles are required, the minimum probe 
volume diameter would be 33 micrometers for this example. 

The final focusing lens is a critical element of the optical system 
because it not only determines the minimum size of the probe volume, but it 
also determines how well the scattered light can be focused at the field 
stop. The field stop (a pinhole aperture usually located in front of the 
photodetector) acts as a spatial filter that prevents light scattered from 
surfaces near the probe volume (flare) from reaching the photodetector, thus 
reducing the r 0 term in equation (1). This ultimately determines how close 
to the end walls and blades that measurements can be made. The lens selected 
for the warm turbine system is a three-element, air-spaced, f/2.5 design, with 
a 100-mm clear aperture and a 235-mm back focal length. 

Another consideration in determining the time required to measure tur- 
bulent flow pararmeters to the required accuracy is that velocity measurements 
taken within the correlation time of the flow are not independent. As shown 
in reference 5 the variance in the measurement of the mean velocity is 



where a^ v is the variance of the individual particle velocity measurements 
(given by eq. (1)), T is the time over which the data is acquired, R is the 
mean data rate, is the variance of the flow fluctuations (the square of 
turbulence intensity), and x c is the flow correlation time. Thus, in addi- 
tion to the data rate, the time scale and turbulence intensity of the flow must 
be known in order to estimate the data acquisition time. Note that increasing 
the data rate will not necessarily reduce the data acquisition time. 

The selection of the type of optical configuration will be based on 
experimental evaluations to be made in the open-jet burner facility. The time 
required to measure the flow parameters given in the Measurement Objectives 
section along with the minimum distance that these measurements can be made 
from solid surfaces are the criteria that will be used to select the optical 
configuration. 


TRAVERSING SYSTEM 

The traversing system for the open jet burner facility is constructed 
using standard, commercially available linear actuators for XYZ motion of the 
optics package. A similar traversing system has been designed for the warm 
turbine system. A two-axis goniometer stage will be used to position a mirror 
located after the final focusing lens. This will permit the optical axis to 
be pointed along a radial line as well as slightly off the radial line to 
minimize shadowing caused by the shape of the turbine blades. A commercial 
controller is to be used for the actuators of the traversing systems; both 
local control and remote computer control via an RS-232 serial interface are 
available. 
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SEED PARTICLES 


Particle characteristics necessary for hot section laser anemometry are 
primarily the same as ambient temperature laser anemometry with the exception 
that hot section particles must retain those characteristics at the high tem- 
peratures. A selected grade of alumina sized at 1.0 micrometer is being used 
initially. Other particles on hand for later trials are titanium dioxide, 
silicon carbide, hollow phenolic spheres, and Georgia clay. 

A commercial aerodynamic particle sizing system is being used to pre- 
sample the particles in use before a test and will be used to sample the 
exhaust products. This system consists of the particle sizer (ref. 6) along 
with a microcomputer, disk drive, and printer. The included software for this 
system is sufficient to analyze and graphically portray all the needed infor- 
mation on particles in the size range that we are using. The system can easily 
be calibrated with a selection of known particles. Samples of the alumina 
have been taken with the particle sizing system and the size distribution 
looks good for our use (fig. 3). No samples of this material have yet been 
acquired after going through the burner. 

The particle generation system consists of a commercial, high volume, 
fluidized bed generator (ref. 7) as shown in figure 4. This system was 
selected so that we could be assured of high particle concentrations for pre- 
sent and future known applications. This system operates by injecting the 
alumina or other particles into a bed of brass beads (approx, diam 100 micro- 
meters). As the bed of beads is fluidized and mixed the particles are de- 
agglomerated and ejected from a 20-mm diameter nozzle. 

Since we are restricted in the diameter of a penetration into the burner, 
we cannot achieve a sufficient flow rate to fluidize the brass beads without 
bypassing some of the particles. The concentration is sufficient using a 3-mm 
inlet into the burner when the particle generator is operating properly. We 
have had clogging problems and two failures in the mixer in approximately 
40 hr of usage. The clogging problem may be caused by moisture in the alumina. 

A seed injection probe has been designed for the warm turbine that can be 
remotely moved in the radial and circumferential directions. The range of 
movement of the seed injection probe takes into account the large turning 
angle (75°) of the stator. 


WINDOWS 

Gaining optical access to the stator and rotor passages of a turbine 
without modifying the flow is a challenging design problem even for cold flow 
at low pressures. The warm turbine facility will use two windows. One will 
cover the stator passage and the other the rotor passage. Only one window 
will be installed at a time. Each window will be curved to match the tip 
radius and will be mounted flush with the inside surface so the flow will not 
be disturbed. The window material is synthetic sapphire, selected for its 
strength at high temperature and good optical properties. Because of their 
curvature the windows will act as cylindrical lenses; the aberration intro- 
duced is minimized by making the windows relatively thin (about 3 mm). 
Compensating optics will be used to correct the remaining aberration. 
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One of the major concerns in the design of the LA for the turbine facili- 
ties is the unknown rate of contamination of the windows. It is expected that 
combustion products, lubricating oil, and seed material will be deposited on 
the window. Tests are being conducted in the open-jet burner facility to 
establish the contamination rate. Various means of keeping the window clean, 
such as using an air purge, will be examined. Window cleaning procedures will 
be studied to determine a practical means of cleaning the window without re- 
moving it from the turbine. A sapphire test window has been mounted on the 
traversing table and we will be watching the change in data rate during burner 
startup and as a function of time at various distances from the hot gases. 

The window will be examined after exposure to determine the type of contami- 
nants. The type of contaminants and rate of buildup data can then be used to 
determine if a high-temperature, on-line window cleaning method is necessary 
and possible. 


DATA PROCESSING AND EXPERIMENT CONTROL 

In the laser anemometry program the computer has become a necessary tool, 
chiefly due to the requirement that data be obtained at rates exceeding 5 kHz 
and that the data be reduced, at least partially, on line. In addition to the 
data acquisition requirements, the computer will become the experiment con- 
troller, performing such housekeeping procedures as deciding when to move the 
optics into position, moving the optics, how long to take data, and which data 
to keep. 

In order to meet these requirements, a computer with a large memory, a 
fast mass storage device, and control capabilities was purchased. The com- 
puter, a PDP 11/44, is now installed in the open-jet burner facility. Using 
the RSX-11M multitasking operating system, the PDP 11/44 will be able to run 
several on-line programs 'simultaneously': acquiring data, reducing and 
storing the acquired data, and presenting a preliminary data summary. In 
addition to the data acquisition and experiment control functions, the PDP 
11/44 will be able to produce high quality on-line graphics, both on the 
terminal and on a graphics printer. A block diagram of this computer system 
is shown in figure 5. 

Until the PDP 11/44 is operational, a small laboratory minicomputer is 
being used to take data with a program which controls the traversing system, 
takes up to 4000 data points, and draws a frequency (or velocity) histogram of 
the acquired data (figs. 6(a) and 6(d)). The data comes in from a counter- 
type signal processor, on a direct memory access (DMA) line, and is stored in 
memory as two 16-bit words per data point. The time-of-flight through the 
required number of fringes (time data) and the time between the present and 
most recent data point (TBD data) are both converted to a FORTRAN readable 
form. If desired, this data can then be stored on disk. Finally the program 
performs the necessary calculations to plot a frequency histogram. 

Even though the on-line data reduction provides a measure of the average 
data rate, in order to properly evaluate the autocorrelation outputs (figs. 
6(c) and 6(f)) a probability histogram of the time between data points is 
necessary (figs. 6(b) and 6(e)). Because the flow is randomly time-sampled, 
data points must occur relatively close together temporally (and therefore 
spatially) for the autocorrelation to be meaningful. If this occurs, as shown 
in both examples, data can then be autocorrelated to obtain a measure of the 
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turbulence scale. (The autocorrelation is calculated using the procedure 
given in ref. 8.) The data shown in figures 6(a) to (c) were taken in the 
potential core of a small air jet (turbulence intensity about 1.6 percent) and 
the data in figures 6(d) to (f) were taken in the mixing region of the open 
jet burner exhaust (cold flow), which had a much higher turbulence intensity 
(about 30 percent). 

The on-line program has recently been modified to take velocity surveys; 
that is, to move the optics under program control, take data, and plot a pro- 
file of the velocity component chosen by the experimenter. Examples of radial 
surveys of mean axial velocity and turbulence intensity taken in the exhaust 
of the open-jet burner are shown in figure 7(a) for cold flow and figure 7(b) 
for hot flow. 

Future modifications will allow the program to take data between the 
blades of a rotating turbine, and keep track of not only the position of the 
optics, but also the position of the turbine rotor with reference to a fixed 
point. 


PLANS 

We plan to install the laser anemometer system in the warm turbine 
facility during 1984 and make it operational in early 1985. Preliminary de- 
sign work for the high-pressure turbine facility (HPT) window is now underway. 


REFERENCES 

1. Goldman, L. J. and Seasholtz, R. G., "Laser Anemometer Measurements in an 
Annular Cascade of Core Turbine Vanes and Comparison with Theory," 

NASA TP-2018, 1982. 

2. Powell, J. A., Strazisar, A. J. and Seasholtz, R. G., "High-Speed Laser 
Anemometer System for Intrarotor Flow Mapping in Turbomachinery," 

NASA TP-2018, 1982. 

3. Edwards, R. V., "Laser Anemometer Optimization," Presented at Turbine 
Engine Hot Section Technology (HOST) Workshop, NASA TM-83022, Cleveland, 
Ohio, October 19-20, 1982. 

4. Lading, L., "The Time-of-Fl ight verses the LDA," Proceedings of the Third 
International Workshop on Laser Velocimetry, Purdue University, 1978. 

5. Edwards, R. V. and Jensen, A. S., "Particle Sampling Statistics in Laser 
Anemometers: Sample and Hold Systems and Saturable Systems," J. Fluid 
Mech., vol . 133, 1983, pp. 397-411. 

6. TSI, Inc. Model APS 33 Aerodynamic Particle Sizer Instruction Manual . 

7. TSI, Inc. Model 9310 Fluidized Bed Aerosol Generator Instruction Manual, 

p. 2. 

8. Smith, D. M. and Meadows, D. M., "Power Spectra from Random-Time Samples 
for Turbulence Measurements with a Laser Velocimeter, " Proceedings of the 
Second International Workshop on Laser Velocimetry, Vol. 1, H. D. Thompson 
and W. H. Stevenson, eds., Bull. No. 144, Purdue Univ., 1974, pp. 27-44. 


62 




DUAL-BEAM FRINGE-TYPE USER ANEMOMETER 






0 



AERODYNAMIC DIAMETER, |im 


Figure 3 

FLUIDIZED BED PARTICLE GENERATOR 


AEROSOL OUT 



7 7 

L MOTOR ^OTOR 


Figure 4 


64 





Figure 5 


65 















CRT DISPLAYS 




CT> 

Ch 


CA-C) LOW TURBULENCE CASE 



(D-F) HIGH TURBULENCE CASE 


CS- -83-2951 


Figure 6 



MEAN AXIAL VELOCITY AND TURBULENCE INTENSITY 


COLD FLOW 



HOT FLOW 


TURBULENCE 

INTENSITY, 

% 



Figure 7 


Page intentionally left blank 



TIME-OF-FLIGHT ANEMOMETER FOR 
HOT SECTION APPLICATIONS 

Robert V. Edwards 
Chemical Engineering Department 
Case Western Reserve University 

The purpose of this section of the project is to design, construct, and test 
laser anemometer configurations for Hot Section velocity measurements. Optimiz- 
ing the laser anemometer system necessarily included the data processing algorithms 
used. It is felt that the requirements here are too demanding for standard laser 
anemometer systems. 


Relevant Hot Section Properties 

1) High temperature with possibility of large background radiation 

2) Difficult optical access 

3) Large flow velocity variation - especially in the rotating sections 

4) Presence of solid surfaces that generate spurious reflections 

5) Low seed particle density 

In the past few years, the laser scattering group at Ris0, Denmark, under the 
direction of Lars Lading, and the laser scattering group at Case Western Reserve 
University under Robert V. Edwards, have worked together to develop procedures for 
the optimal design of laser anemometry systems. The principles derived are being 
used to design the system for Hot Section measurements. 

The system decided on is a so-called ti me- of-f light anemometer with elliptical 
spots. The version of the time-of-fl i ght designed for this project contains two 
new features: 1) Elliptical spots - this gives the wide flow angle acceptance 

characteristics of a "fringe" anemometer combined with the superior spatial reso- 
lution of a time-of-flight anemometer. 2) The prototype for the Hot Section 
measurements uses a unique optical coding to transform the pulse into the optimal 
form for pulse position sensing. Heretofore, this required rather complex and 
inflexible electronic circuitry. This optical processor is intrinsically free from 
some of the errors to which the electronic circuits were prone. 

Figure 1 is a photograph of the prototype optical system. The system is 
presently undergoing tests for accuracy and spatial resolution. The preliminary 
indications are that the system is capable of measurement within 300 ym of a 
surface. A better estimate of the spatial resolution will be possible only after 
we mechanically stiffen the system. It requires very tight tolerances on the optics. 

A prototype has been built of an electronic signal processor for the anemom- 
eter. Figure 2 shows the signals generated to detect the pulse position. The zero 
crossing of the s-shaped pulses are used. A velocity histogram obtained using this 
system on a seeded flow in a small wind tunnel is shown in Figure 3. 

The prototype system is being evaluated in terms of scale-up to a system 
capable of the desired hot stage measurements. In particular, the mechanical and 
optical requirements of the system are being evaluated. 
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Figure 1: Prototype of Time-of-Fli ght Anemometer. It is constructed of 

commercial optical "breadboard" components. 
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Figure 2; Oscilloscope trace of difference signal from the two photo- 
detectors . 



Figure 3: Velocity histogram obtained using new anemometer on air flow in 

a small wind tunnel. The pulse separation is converted to a pulse 
height and then stored on a pulse height analyzer. 
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TURBINE HEAT TRANSFER 


John E. Rohde 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


Improved turbine durability and performance and reduced development cost 
will all result from improved methods of predicting turbine metal temperatures. 
As you know, better metal temperature prediction methods require improvements 
in the methods of predicting the hot gas flow over the turbine airfoils and 
the cooling air flow inside the airfoil and in the methods of predicting the 
heat transfer rates on both the hot gas-side and coolant-side of the airfoil. 
The overall HOST Turbine Heat Transfer effort is directed at improving all 
four of these areas of concern. 

Achievement of these improvements requires a rigorous and systematic 
research effort from both the experimental and analytical sides. The experi- 
mental approach being pursued starts with fundamental experiments using simple 
shapes and flat plates; progresses on to more realistic cold, warm, or hot 
cascades using airfoils; continues to progress on to more realistic warm tur- 
bine, large low-speed turbine, or transient turbine tests; and finally com- 
bines all the interactive effects in real-engine or real-engine type turbine 
tests. Analytical approaches being pursued also start with relatively simple 
mathematical models and progress to more realistic cases that include more 
interactive effects, and finally combines all the interactive effects of the 
turbine operating in the real engine environment. 

The HOST Turbine Heat Transfer Subproject schedule showing the current 
activities is included in the attached figures. Each line on tne figure re- 
presents a separate contract, grant, or NASA Lewis in-house effort. The dotted 
blocks indicate future potential contract or grant activities. Currently, 
contract and grant efforts are being conducted to obtain fundamental experi- 
mental data and to develop and/or assess analytical tools in all four areas of 
concern. These contract and grant activities will be discussed in detail 
later in this meeting by the respective principal investigators. 

NASA Lewis in-house turbine research efforts are being pursued to obtain 
more realistic warm turbine and real-engine type turbine experimental data. 

The major turbine research parameters of interest that will be measured or 
determined to provide this better understanding of the aerothermal loads on 
air-cooled airfoils are the following; 

1) local hot gas flow velocities and secondary flows, 

2) local hot gas recovery temperature along the airfoil surfaces and other 
hot gas path surfaces, 

3) local total and radiation heat fluxes to the airfoils and other hot gas 
path surfaces, and 

4) local airfoil wall temperatures. 

Very little or none of this type experimental research information exists with 
controlled warm or real-engine type conditions and known boundary conditions. 
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These in-house turbine research efforts will be conducted using the best 
available analyses to help define the test configurations, the types of 
research measurements, and for comparisons with the measured research results. 
Initially, analytical efforts will use the best available flow and heat trans- 
fer codes such as a twoor three-dimensional inviscid flow code and a two- or 
three-dimensional boundary layer heat transfer code. These analyses will be 
applied at the mid-span section and possibly at the hub and tip sections or 
other local zones of the hot gas flow path. More sophisticated three- 
dimensional viscous codes and three-dimensional viscous codes with boundary 
layer resolution will be used as they become available. These analytical 
efforts will be conducted using the best available source or sources in-house 
and on contract with industry and universities. 

An Airfoil Metal Temperature Runoff effort is planned to start in 1985. 
This Airfoil Metal Temperature Runoff effort will be initiated using an all- 
impingement cooled static airfoil configuration and a multipass cooled rotating 
airfoil. Both of these air-cooled airfoils will use NASA conceived profiles 
and no film cooling initially. This program will provide engine manufacturers 
with an opportunity, through contracts and/or cooperative agreements, to assess 
and improve their analytical tools used to predict airfoil metal temperatures. 
The results of these assessments will be accumulated as an industry uncertainty 
band from the experimental data with the individual contributors predictions 
remaining known only to the individual contributor. Contracted contributors 
would be required to identify the analytical tools used and the uncertainty 
levels obtained. 
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AIRFOIL METAL TEMPERATURE RUNOFF 

GOAL: ASSESS AND improve analytical tools used to predict airfoil 
METAL TEMPERATURES 


• ASSESS METAL TEMPERATURE PREDICTION 
USING AN ALL-IMPINGEMENT COOLED VANE 
WITH NO FILM COOLING INITIALLY 

• REVISE ANALYTICAL TOOLS BASED ON 
REAL-ENGINE TYPE MEASUREMENTS OF 
THE LOCAL HOT GAS TEMPERATURES AND 
HEAT FLUXES AND REASSESS VANE METAL 
TEMPERATURE 

• ASSESS METAL TEMPERATURE PREDICTION 
USING A MULTI-PASS COOLED BLADE WITH 
NO FILM COOLING INITIALLY 


ANALYTICAL TOOLS USED TO PREDICT AIRFOIL 
METAL TEMPERATURES 



mU. AWO*. METAl TEMPERATURES 


• REVISE ANALYTICAL TOOLS BASED ON REAL-ENGINE TYPE MEASUREMENTS OF THE 
LOCAL HOT GAS TEMPERATURES AND HEAT FLUXES AND REASSESS BLADE METAL 
TEMPERATURES 

STATUS : mulitple contracts and/ or cooperative agreements to be 

AWARDED IN 1985 
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GAS FLOW ENVIRONMENTAL AND HEAT TRANSFER 
NONROTATING 3D PROGRAM 


T. Gei 1 and J. Steinhoff 
University of Tennessee Space Institute 

OBJECTIVES 

The experimental contract objective is to provide a complete set of "benchmark" 
quality data for the flow and heat transfer within a "large" rectangular turning 
duct. These data are to be used to evaluate and verify three-dimensional internal 
viscous flow models and computational codes. The analytical contract objective is to 
select such a computational code and define the capabilities of this code to predict 
the experimental results. Details of the proper code operation will be defined and 
improvements to the code modeling capabilities will be formulated. 

The experimental and analytical efforts are being conducted under a coordinated 
multiphase contract. Phase one, the current work, is the study of internal flow in 
a rectangular, square cross-sectioned, 90° bend turning duct, and is planned as a 28 
month investigation which started in April, 1982. Phase one is divided into the 
following five tasks: I. Design and Fabrication, II. Experimental Velocity Mea- 

surements, III. Experimental Heat Transfer, IV. Theoretical Analysis and Data 
Comparison, V. Reporting and Technical Data. Future work to be performed at NASA's 
option includes the investigation of flow over an airfoil cascade, with and without 
film cooling, inside the turning radius of the duct. 

Separate but coordinated experimental and analytic approaches are in progress 
to attain the contract objectives. 

APPROACH, EXPERIMENTAL 

The experimental facility design features modular tunnel components which allow 
flow measurements every 15° in the 90° bend. The 25.4 cm (10 in) square cross sec- 
tion tunnel is designed with a 13 to 1 area ratio bell mouth contoured to provide 
uniform flow velocity and is powered by a variable speed, six-bladed fan. (See 
Figures 1 and 2) . The tunnel is designed for incompressible flow and will produce 
Reynolds numbers of 0.2 to 2.0 X 10^ at the entrance of the 90° bend for tunnel velo- 
cities of 6 to 30 m/sec (20 to 100 ft/sec) . These two flow conditions provide 
laminar and fully turbulent boundary layer profiles at the entrance to the 90° bend. 
The facility is also designed for adiabatic wall testing with large thermal gradients 
in the air stream. Heated air will be provided by electric resistance heaters sized 
so that 100 KW will produce a minimum of 110°C (200°F) temperature increase in the 
air stream. 

The primary instrumentation is designed for non-intrusive flow measurements 
utilizing a three-dimensional, laser velocimeter (LV) and wall static pressure and 
heat flux gages. The LV utilizes two color beams and Bragg diffraction beam split- 
ting/frequency shifting to separate the three simultaneous, orthogonal, vector 
velocity components. The LV signal processors determine the digital values of velo- 
city from the seed particles crossing the laser beam probe volume. To improve and 
speed up digital data acquisition, the LV processors are designed around an S-100 


Work done under NASA contract NSA3-23278. 



bus Z-80 microprocessor which provides on-line, near-real time data reduction. This 
on-line data reduction capability will be used to assess the adequacy and precision 
of the data as it is acquired and recorded for off-line detailed analysis. Typical 
data output is shown in Figure 3. To qualify the measurements as "benchmark" data, 
the LV data will be compared with both pitot probe and hot wire anemometer measure- 
ments for flow conditions which permit comparisons. 

CURRENT RESULTS AND PLANS, EXPERIMENTAL 

All experimental hardware and instrumentation systems have been designed, fabri- 
cated, assembled, installed and checked out. Only the air heater remains to be 
fabricated, and only the heat flux measurement system remains to be checked . Exten- 
sive tunnel flow quality surveys have been completed in the inlet section following 
"benchmark" calibration of the LV system, pitot-static, and hot wire anemometer. All 
three measurement systems were calibrated against reference standards and were found 
to agree within i 1% on the entrance section velocity. The LV system was calibrated 
with a spinning disk reference velocity at 20 m/sec. The pitot-static pressures were 
read on a precision slant manometer with 13mm of H 2 O full scale. The hot-wire system 
was calibrated in a reference nozzle flow at 30.5 m/sec. Velocity surveys taken 50.8 
cm (20 in) behind the bell mouth exit showed flat velocity profiles i 0.5% mean velo- 
city outside of the boundary layer. Laminar and turbulent boundary layers were 
measured at velocities of 4.5 m/ sec and 19.4 m/sec, respectively as shown in Figure 
4. Most significant finding was that the LV and hot-wire turbulence intensity mea- 
surements agreed within 10% or nominally 1% turbulence intensity. The mill bed 
traverse system has demonstrated repeated accuracy of i 0.1mm on all three axis of 
movement. Both LV and probe positions are controlled by the computer driven mill 
bed. 


Development of the flow seeding system was successfully completed during the LV 
system checkout. Phenolic micro-ballons of 2-5 micron size are sprayed in a slurry 
of alcohol and water into the air stream in front of the bell mouth. Uniform seed 
distribution was obtained in the test section with this system. 

The experimental effort is progressing toward completion of flow measurements 
at the seven stations in the 90° duct for the unheated flow. A minimum of 300 
spatial points in the duct half-plane are being surveyed for mean velocity, unsteady 
velocity and total pressure at each station. After completion of the unheated flow 
surveys at two Reynolds numbers, scheduled for 1 December 1983, the air heater and 
tunnel insulation will be installed for the adiabatic wall testing during the winter 
months when lower ambient temperatures will reduce tunnel wall temperatures . The 
experimental phase of this effort is scheduled for completion in May, 1983. 

APPROACH, ANALYTICAL 

The analytical approach involves, first, selecting a computer code capable of 
solving the Navier Stokes equations with turbulence models for three dimensional 
internal flow, and adapting it to the experimental geometry and flow conditions. 

After this, calculations are to be made for laminar flow conditions for unheated 
flow. Analysis of these calculations will define the grid size and stretching fac- 
tors required for adequate resolution as well as the values of time steps and 
smoothing factors required for convergence. Also, any output and graphics capability 
required for comparison with data is to be developed during this phase. The 
adequacy of the code with respect to the differencing scheme, adaptability and con- 
vergence will be decided in this phase, by comparisons with published experimental 
and computational data and by grid sensitivity studies. 
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The next phase of the effort involves computation of heated laminar flows with 
adiabatic wall conditions. This calculation is important for determining maximum 
wall temperatures to be expected in the experiment. It will also be important in 
determining whether for the temperatures planned, the velocities will deviate suffi- 
ciently from the unheated case to provide an interesting flow. 

The third phase involves detailed comparison with experiment for laminar flows 
with and without heat added. The resolution of the computer code will be the main 
aspect tested here. The heated case will involve adiabatic wall conditions. 

The fourth phase involves computing turbulent flows corresponding to actual 
experimental conditions, with and without heat added. Particular attention will be 
paid to the treatment of heat transport in the turbulence model. Detailed compari- 
sons with experiment will be made. 

A final phase involves selective grid refinement in regions where high resolu- 
tion is required so that a measure of the numerical truncation error can be obtained. 

CURRENT RESULTS AND PLANS, ANALYTICAL 

The code selected for the analytical study is a version of the Beam-Warming 
algorithm adapted to generalized coordinate systems by P. D. Thomas. It is very well 
documented and relatively easy to use. The original plan to use the "MINT" code 
developed by Briley and McDonald has not been followed due to our inability to obtain 
the code. Detailed comparisons of the results of our code, for a laminar flow case, 
with published results of the MINT code have been made. Eor comparable grids the 
codes gave similar results. These results have also been compared to published 
experimental results as shown in Figure 5. This comparison indicated a discrepency 
which appeared to be due to lack of adequate resolution in the computed results for 
both codes. This discrepency has been investigated by grid refinement studies. 
Further, the effect of the grid stretching factors on the solution has been deter- 
mined. The time step size and smoothing factors have been determined which 
approximately optimize convergence for this case. A graphical display program has 
been written and used to simulate streamlines in the flow as shown in Figure 6. 

These results have been qualitatively compared with preliminary experimental results 
in a 1/3 scale duct. 

A heated laminar flow case with adiabatic wall conditions was then computed. It 
was determined that the wall temperature at a point near the end of the bend would be 
close to the maximum inlet temperature, which occurs at the center of the duct. 
Computed total temperature profiles are shown in Figure 7. This result is important 
in determining the maximum allowable inlet temperature. For the temperature profile 
considered, which was close to that planned in. the experiment, it was also determined 
that the velocity profiles would differ significantly from those of the unheated 
flow. Thus, the heated flow experiment would result in a qualitatively different 
flow, for the planned inlet temperatures, and the experiment would give useful 
results. 

Currently, we are running unheated turbulent flow cases. Preliminary results 
are being studied to determine the sensitivity of the results to the turbulence 
model parameters. Detailed comparisons with experiment will be made for the laminar 
cases, with and without heating, as soon as they are available. We will then 
proceed with the comparisons for turbulent flow and additional grid refinement 
studies. 
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GAS SIDE HEAT TRANSFER 


Larry D. Hylton 
Allison Gas Turbine Operations 
General Motors Corporation 

Improvements in methods for predicting heat transfer rates on the hot gas side 
of turbine airfoils are necessary for improved turbine durability and performance. 

The development and verification of improved analytical models requires a systematic, 
closely coupled experimental and analytical program. 

Work in this area has been performed under two contracts. The first, NAS 3-22761, 
addressed the problem of the prediction of hot gas side heat transfer rates to 
internally cooled (non-film cooled) airfoils. This effort was completed this year. 

The second contract, NAS 3-23695, is currently underway and will investigate the 
effect of leading edge "showerhead" film cooling on downstream heat transfer rates. 

The objectives of the first contract were to assess the capability of currently 
available modeling techniques to predict airfoil surface heat transfer distributions 
in a 2-D flow field, acquire experimental data as required for model verification, 
and to make and verify improvements in the analytical models. 

Two data sets, Turner (ref. 1) and Lander (ref. 2), were selected from the 
literature for use in evaluating existing models. Two additional airfoils were 
chosen for cascade testing under this contract. These airfoils, designated the 
Mark II and C3X, are representative of highly loaded, low solidity airfoils currently 
being designed. The aerodynamic configurations of the two vanes were carefully 
selected to emphasize fundamental differences in the character of the suction 
surface pressure distributions and the consequent effect on surface heat transfer 
distributions. Cross sections of the four airfoils and the grid used to make 
inviscid flow predictions for each airfoil are shown in figure 1. Note the 
significant variation in airfoil geometry. This variation was intended to provide 
a significant test of the analytical models. Predicted surface pressure distributions 
for the four airfoils are shown in figure 2. 

The two heat transfer cascades tested were run in the Allison Aerothermodynamic 
Cascade Facility (ACF) . The facility, described in figure 3, provides the capability 
of obtaining both heat transfer and aerodynamic measurements at simulated engine 
conditions. The experimental measurements were made in moderate- temperature , three- 
vane cascades under steady-state conditions. The principal independent parameters 
(Mach number, Reynolds number, turbulence intensity, and wall-to-gas temperature ratio) 
were varied over ranges consistent with actual engine operation, and the test matrix 
was structured to provide an assessment of the independent influence of each parameter. 
The test matrix over which both cascades were operated is shown in figure 4. Data 
was obtained at two exit Mach numbers, 0.9 and 1.05, and over a range of exit Reynolds 
numbers from 1.5X10^ to 2.5X10^. The inlet turbulence intensity, Tu, and wall-to-gas 
temperature ratio, T w /Tg, were also varied. 

The method employed in the facility to obtain airfoil surface heat transfer 
measurements is shown schematically in figure 5. Basically, the exterior of the 
airfoil is instrumented with grooved surface thermocouples with this data serving 



as the exterior boundary condition input to a finite element analysis. The internal 
boundary conditions are calculated heat transfer coefficients in the internal cooling 
holes. These values are calculated from measurements of the coolant temperature, 
pressure and flowrate. The instrumented airfoil also contains static pressure taps, 
thus permitting simultaneous measurement of the surface pressure and heat transfer 
distributions. A typical set of data for the C3X cascade illustrating the effects 
of exit Reynolds number on the heat transfer distribution is shown in figure 6. 

Note that the heat transfer measurement technique provides sufficient detail to 
clearly see the effect of Reynolds number on transition on the suction surface. 

Data from the two cascades, coupled with that from the literature cases, provide a 
data base covering a wide range of operating conditions and geometries and thus 
present a significant test for the predictive capabilities of the analytical methods. 

The analytical methods development program consisted of two separate phases. 

In the first phase, the literature was reviewed to establish general candidate methods 
that were characteristic of current methodology incorporated within actual gas turbine 
preliminary design systems. As a result of this survey, three 2-D boundary layer 
methods were chosen: an integral method, a finite difference (differential) method 

with a zero-equation mixing length hypothesis turbulence model, and the same differ- 
ential method with a two-equation turbulence model. The literature was thoroughly 
reviewed to obtain relevant airfoil heat transfer experimental data to use in a 
general evaluation of the three selected boundary layer methods. The data sets were 
selected based on relevance to realistic gas turbine environments (i.e., Reynolds 
number effects, free-stream turbulence effects, strong pressure gradient effects, 
etc.). Analytical/numerical solutions were compared with experimental results. 

Based on the findings of this first phase general methods evaluation process, the 
differential method with zeroth order turbulence modeling was selected for the 
second phase of the analytical program. The literature was further reviewed for 
models that had the potential of treating the airfoil heat transfer problem more 
realistically. A number of transition process models, free-stream turbulence 
augmentation models, and a single explicit longitudinal surface curvature correction 
model were selected for evaluation. Various single and/or combined model solutions 
were evaluated using data from four different airfoil experiments. This evaluation 
process eventually led to a final "gas turbine airfoil specific" modeling effort 
which resulted in an effective viscosity formulation that, when implemented, gave 
better overall solutions than any literature modeling approach tested previously. 

Figures 7a and 7b respectively show the unmodified and modified suction surface 
heat transfer predictions for three different operating conditions using the STAN5 
boundary layer code for the Lander airfoil. Increasing run numbers correspond to 
increased inlet or exit Reynolds number and free-stream turbulence intensity. The 
experimental data are represented as symbols. Lander's data are important in that 
they illustrated nominally laminar heat transfer augmentation attributed to free- 
stream turbulence effects, as well as Reynolds number effects related to transition 
origin. As shown in figure 7, the augmentation phenomenon is predicted significantly 
better by the final model. 

Figures 8a and 8b show the unmodified and modified solutions compared with the 
data of Turner. The significance given to Turner's data was that they isolated the 
effects of free-stream turbulence. Figure 8a shows only one solution because the 
original unmodified method did not account for the effects of free-stream turbulence. 
As can be seen in figure 8b, the modified solutions give a very good representation 
of the pressure surface experimental data. The modified suction surface solutions 
give reasonable trends up to the point where a transition process is indicated by the 
experimental data. Overall, the modified solutions are a significant improvement 
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over the unmodified solution, represent the pressure surface data very well, and provide 
qualitatively good trends for the suction surface. 

Unmodified and modified predictions of the characteristic Reynolds number effects 
in the Mark II cascade are compared with the data in figures 9a and 9b, respectively. 

It should be pointed out that the analytically predicted stagnation point was dis- 
placed approximately 5% (0.05) of pressure surface distance toward the pressure 
surface away from the extreme forward point on the airfoil, which was used as the 
datum (0) in these figures. The stagnation point corresponds to the predicted 
inviscid flow solution zero velocity location on the pressure surface. Note that 
this does not correspond to the highest local value of measured heat transfer in the 
leading edge region. Both the modified and unmodified solutions reflect the proper 
trends moving away from the stagnation point. The absence of solutions beyond 0.2 
normalized surface distance on the suction surface indicates that all solutions 
encountered separation due to the presence of a suction surface shock at the 
location. No attempt was made to restart the solutions downstream of the shock. 

Overall, the modified solutions are able to qualitatively and quantitatively predict 
the pressure surface data reasonably well and yield much better predictions than 
the unmodified solutions, which predicted pressure surface transition. 

In a manner similar to the Mark II comparative studies, the experimental results 
for the C3X cascade were also simulated numerically and the predictions are shown 
in figure 10. Figures 10a and 10b show both unmodified and modified solutions at 
three different Reynolds number conditions. Qualitatively, the modified pressure 
surface solutions represent a substantial improvement over the original (unmodified) 
approach. However the quantitative predictions (using the modified procedure) begin 
to deviate significantly from the data along the aft portions of the surface. The 
suction surface predictions of both the unmodified, figure 10a, and the modified, 
figure 10b approaches yield quantitatively acceptable results for some of the cases, 
but the indicated suction surface transition process (i.e., gradual transition) is 
better represented by the modified solutions. 

Finally, in response to the objectives of this program, a recommended procedure 
was developed for constructing a viable, 2-D airfoil external convective heat 
transfer method for gas turbine design systems, including the specification of 
boundary conditions, initial conditions, and preferred definitions of effective 
viscosity determined here to be most suitable for gas turbine preliminary design 
applications . 

The analytical and experimental work performed under contract NAS 3-22761, 
including the recommended design procedure, are reported in detail in NASA CR 168015 
(ref. 3) which was published in May, 1983. 

Work on the second contract, NAS 3-23695, began earlier this year. This effort 
is intended to extend the work performed under the first contract in two respects. 

First, the analytical boundary layer analysis and experimental cascade studies of the 
first contract will be extended to include a leading edge showerhead film cooling 
array. Secondly, recognizing the long term limitations of the boundary layer approach, 
an analytical effort was added to investigate the application of a Navier-Stokes 
solver to the turbine cascade problem. This effort was subcontracted to Scientific 
Research Associates, Inc. (SRA) and will utilize their MINT code. 

The boundary layer efforts are structured similar to the first contract. Namely 
the STAN5/STANC00L type approach to making heat transfer predictions for airfoils 
with leading edge film cooling will be evaluated and modified as required to improve 
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heat transfer predictions. These modifications will then be validated by comparisons 
with experimental data acquired under the contract. 

The experimental program utilizes the C3X cascade from the original contract. 

The center, instrumented airfoil was replaced with one containing five rows of film 
cooling holes in the leading edge region. Figure 11 illustrates the original 
airfoil with the radial cooling holes and finite element grid shown. The new 
film cooled airfoil was modified in the leading edge region as shown in figure 12. 
Downstream of the region shown in figure 12 the cooling hole geometry and instru- 
mentation is identical to the non-film cooled airfoil tested under the first 
contract. Heat transfer measurements will be made downstream of the film cooling 
array. No heat transfer measurements will be made within the array. By using a 
modified airfoil of the same profile as one of the two tested under the original 
contract, a good non-film cooled baseline will have been established, and fully 
qualified cascade hardware and experimental techniques will be employed to maximum 
advantage . 

The analytical efforts on the Navier-Stokes MINT code are intended to provide 
a major step toward developing an analytical tool capable of predicting the flow 
and heat transfer in a full 3-D turbine cascade. Initial efforts on the MINT code 
will be performed on a 2-D version. The "C" grid generator previously used in the 
program will be replaced with an "0" grid generator. Due to its construction, the 
"C" grid requires a cusped trailing edge. In solving the transonic turbine problem, 
the trailing edge geometry appears to be very important. In preliminary calculations 
with the "C" grid, discrepancies between predictions and experimental data in the 
trailing edge region are thought to be a result of the cusped trailing edge approxi- 
mation. The addition of the "0" grid generator to the program should enhance its 
capabilities. 

The 2-D version of the code with the "0" grid generator will be used to make 
predictions for the cascades tested under the first contract. The comparison of 
the predictions with the data will serve as verification of the code modifications. 
Following these comparisons the code will be extended to handle the full 3-D case 
and a sample calculation for a 3-D cascade will be made. The final modification 
scheduled to be made to the code is the incorporation of a film cooling-trailing 
edge blowing capability into the code. Following completion of all modifications 
to the code, SRA will assist NASA personnel in running demonstration cases with the 
code on the NASA Lewis computer. 

Progress on the second contract has been made on both analytical efforts and on 
the experimental effort. The boundary layer analytical efforts have concentrated 
on evaluating the STAN5/STANC00L programs abilities to predict heat transfer to 
circular cylinders. This is the initial step in predicting the airfoil leading edge 
film cooling problem. Comparisons of the predictions for both solid cylinders and 
film cooled cylinders were made with the data in reference 4. Based on these results 
efforts are currently underway to make modifications to the models to improve their 
capabilities . 

Progress on the MINT code includes installation and checkout of the program on 
the NASA Lewis computer. In addition, the "0" grid generator has been developed and 
sample "0" grids have been constructed. Efforts on verifying the capabilities of 
the program with the "0" grid generator by comparing predictions with cascade data 
sets are currently underway. 

Under the experimental phase of the contract the new C3X airfoil containing 
the leading edge "showerhead" film cooling array, has been designed, fabricated, 
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instrumented and installed in the cascade. Preliminary testing of the cascade 
has started. Technical efforts on this program are scheduled to be completed 
by May, 1984. 
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Figure 2. Airfoil predicted pressure distribution. 
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OPERATIONAL SPECIFICATIONS 
AIR SUPPLY 

INLET PRESSURE 
STAGNATION TEMPERATURE 
PRIMARY FUEL 


I — CONTROL CONSOLE 

9.5 LB/SEC AT 105 PSIA OR 
5.0 LBM/SEC AT 245 PSIA 

20 PSIA TO 245 PSIA 
400° F TO 3200°F 
NATURAL GAS 


OPERATIONAL PHILOSOPHY 
o TWO-DIMENSIONAL LINEAR CASCADE 

o PROVIDES HEAT TRANSFER AND AERO- 
DYNAMIC DATA SIMULTANEOUSLY 

o OPERATED AT SIMULATED ENGINE 
CONDITIONS 

o REDUCED TEMPERATURE 
o REDUCED PRESSURE 
o SCALED-UP AIRFOIL GEOMETRY 
o HIGH FREESTREAM TURBULENCE 
o WIDE OPERATING RANGE 

o REYNOLDS NUMBER CONTROL 
o EXIT MACH NUMBER CONTROL 

o WALL-TO-GAS TEMPERATURE RATIO 
CONTROL 

o INLET TURBULENCE INTENSITY 
CONTROL 

o HIGH DENSITY INSTRUMENTATION 
o UP TO 300 TEMPERATURES 
o UP TO 288 PRESSURES 
o DEDICATED FACILITY COMPUTER 

o COMPUTER CONTROLLED DATA 
ACQUISITION 

o ONLINE DATA ANALYSIS 


Figure 3. Allison Aerothermodynamic Cascade Facility 




Figure 4. Cascade test matrix. 
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Figure 5. Heat transfer measurement technique. 
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Figure 6. Reynolds number effect on surface heat transfer 
distribution for C3X airfoil. 
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(b) Modified STAN5 results 


Figure 7. STAN5 solutions compared with Lander airfoil suction surface 
experimental heat transfer coefficient data illustrating the combined 
effects of varying Reynolds number and free-stream turbulence intensity. 
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(b) Modified STAN 5 results 


Figure 9 . STAN5 solutions compared with Turner airfoil experimental 
heat transfer coefficient data illustrating the effects 
of varying free-stream turbulence intensity. 


96 



?ua<ii 4snc*0* 
HA^kTI 441»'43» 
MARk.II «3fl(46) 


RUN 40 4$» I 

RJ<4 43 441 t 

RUN 48 43« f 


HO* 1135 watts/M /K 
(200 Btu/hrffl 2 /*F) 


> , s 

•*f\ .• .;••• 

A'W"-.;. 

\ ♦V+X.* +4 ! st I 


' .-obb* . | + 

rfA*? % ; 

a> ~Xkk / 


?°? + ? JE 


S'; /\. ' 

V v < 


B ** • 3* ♦ 

O 04 _ D$ 

00 B, + • s 

+* B 

t * B + 0 »« 

1 4 B * 4 _ 


SURFACE DISTANCE S/ARC 


TC82-6331A 


(a) Unmodified STAN5 results 


JIARKIJ 4SMC40) 
HAR<TI 44 M f 43' 
NARtJI 43XC46) 


RUN 40 4Slf 

RUN 43 44 11 

RUN 46 431 1 


HO -1135 watts/M /K 

(200Btu/hr/ft 2 rFI 


s. 

* 3 fl A V 


s 

♦ a * 

4 BB 

4 * + 4 


Li Vi 

k 


♦ ♦ *4 Q 


*♦ * 

9 0 4 

<W 0 B \ B 

V t B * B b«4 *• 

+ B* + 


.4 

SUCTION 


SURFACE DISTANCE S/ARC 


(b) Modified STAN5 results 


TE82-6331B 


Figure 9. STAN5 solutions compared with Mark II airfoil experimental 
heat transfer coefficient data illustrating the effects 
of varying exit Reynolds number. 
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Assessment of a 3-D Boundary Layer Code to Predict 
Heat Transfer and Flow Losses in a Turbine* 

Veer N. Vatsa 

United Technologies Research Center 

The prediction of the complete flow field in a turbine passage is an extremely 
difficult task due to the complex three-dimensional pattern which contains separa- 
tion and attachment lines, a saddle point and horseshoe vortex (Fig. 1). Whereas, 
In principle such a problem can be solved using full Navier-Stokes equations. In 
reality methods based on a Navier-Stokes solution procedure encounter difficulty In 
accurately predicting surface quantities (e.g. heat transfer) due to grid limita- 
tions imposed by the speed and size of the existing computers. On the other hand 
the overall problem Is strongly three-dimensional and too complex to be analyzed by 
the current design methods based on inviscid and/or viscous strip theories. Thus 
there is a strong need for enhancing the current prediction techniques through 
inclusion of 3-D viscous effects. A potentially simple and cost effective way to 
achieve this is to use a prediction method based on three-dimensional boundary layer 
(3-DBL) theory. The major objective of this program is to assess the applicability 
of such a 3-DBL approach for the prediction of heat loads, boundary layer growth, 
pressure losses and streamline skewing in critical areas of a turbine passage. A 
brief discussion of the physical problem addressed here along with the overall 
approach is presented in the following paragraphs. 

In the present investigation, zonal concepts are utilized to delineate regions 
of application of 3-DBL theory — these being the endwall, suction and pressure sur- 
faces. Each of the regions selected for this investigation has some unique 
features. For example, the experimental data of Ref. 1 for the surface streamline 
pattern (Fig. 2a) and the corresponding Stanton (St) number distribution (Fig. 2b) 
for the endwall region of a planar cascade displays strong three-dimensional effects 
due to sweeping of the boundary layer across the passage from the pressure to the 
suction surface. A modified version of the 3-DBL code of Ref. 2, named "TABLET" 
(Ihree- Dimensional Algorithm for Boundary-Layer .Equations in Iurbulent Flows) will 
be used to analyze the viscous flow downstream of the attachment line between the 
pressure and suction surface junctures in the endwal I region. 

The second region of interest is the turbine blade suction surface. As shown 
in Fig. 3a, the growth of endwall boundary layers produces Inward deflection of the 
streamlines along the suction surface. Thus strong crossflow velocities are induced 
from the endwall region towards midspan and the flow becomes progressively more 
three-dimensional as the trai I i ng-edge is approached. The effect of cross flow 
Induced three-dimensionality is also clear from the measured St number distributions 
(Ref. 1) shown in Fig. 3b. The "TABLET" code will be used to analyze the flow and 
predict the effect of streamline convergence on St number distributions in the 
region which lies downstream of the leading-edge and between the separation line and 
midspan. 


*This work was sponsored by NASA-Lewis Research Center under Contract 
NAS3-23716. 
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The final region of interest is the pressure side of the blade. Here, the 
three-dimensionality In the flow is induced by a different mechanism altogethei — 
namely, that of blade rotation. A typical streamline pattern on the pressure 
surface of a rotating blade obtained from Ref. 3 is presented in Fig. 4, where a 
strong outward radial flow was encountered. Due to the lower flow velocities In the 
viscous- layer, the flow close to the wall is skewed outwards more than the inviscid 
flow. The flow downstream of the leading-edge and away from the hub will be 
analyzed to assess the "TABLET" code for predicting radial outflow of surface 
streamlines due to blade rotation effects. 

The computer code "TABLET" being used In this Investigation solves the finite- 
difference form of the compressible 3-DBL equations (including the energy equation) 
In a nonorthogonal surface coordinate system. An efficient. Implicit, fully coupled 
finite-difference solution procedure is employed. Boundary conditions are obtained 
from experimental data to eliminate errors associated with inviscid approximations. 
Starting solutions along two inflow boundaries (selected from close examination of 
data) are obtained by solving the appropriate limiting form of the 3-DBL equations. 
Sample solutions obtained to date will be described in the next few paragraphs. 

The first test case for this investigation was selected to establish the 
accuracy of the present 3-DBL code for analyzing passage flows. For this reason, 
the flow In a 60 deg curved duct, for which an extensive set of data is available in 
Ref. 4, was computed using the "TABLET" code. The schematic of the flow problem 
along with the nomenclature and locations of the measuring stations is shown In Fig. 
5. The two inflow boundaries for this problem were selected to be along line A and 
along the radial line at Station 1. Streamwise velocity profiles along these Inflow 
boundaries were generated from Whitfield’s wall-wake correlation (Ref. 5) which uses 
experimentally measured values of the skin friction coefficient, momentum-th i ckness- 
Reynolds number and shape factor as input. Cross-flow velocity profiles along the 
inflow boundaries were generated from Mager’s (Ref. 6) cross-flow representation 
using measured skewing of surface streamlines. The "TABLET" code has been used to 
compute the viscous flow on the upper wail of the curved duct from Station 1 to 
Station 15 in the streamwise direction and from the pressure side A to the suction 
side E In the spanwise direction. The predicted values for all the integral 
properties and the skin friction coefficient compare very well with the measured 
data. A typical comparison from this study is presented in Fig, 6 where the surface 
streamline skewing angle, is shown. This figure clearly indicates that 
relatively large skewing of surface streamlines, typically encountered In turbine 
passage flows, can be accurately predicted using the present 3-DBL analysis. 

The next test case chosen for this study Is the flow in the endwall region of a 
turbine passage. The surface streamline pattern and St number distribution for this 
case have already been shown In Figs. 2(a) and (b). Figure 7a shows the schematic 
diagram of the computational mesh which consists of boundary fitted lines at con- 
stant percentage pitch locations and vertical lines at constant percentage axial 
locations. The upstream Inflow boundry is selected to be the line at 10 percent 
chord distance downstream of the leading edge and the pressure side Inflow boundary 
is taken to be along the Intersection of endwall and pressure surfaces (Fig. 7a). 
Locally similar solutions are used along these inflow boundaries to generate the 
starting profiles. Since the location and extent of the transition zone is not 
known for this problem, a preliminary case has been run by assuming that the flow is 
fully laminar in the leading-edge region and it transits to turbulent flow instan- 
taneously at the second streamwise mesh point (x/Bx = 0.105). The computed Stanton 
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(St) number distribution at x/Bx = 0.2, 0.4, 0.6 and 0.8 is shown In Fig. 7b along 
with the measured data. The predicted values of St number display the correct 
qualitative variation when traversing from pressure side to suction side over the 
region considered here. Keeping in mind the fact that no attempt has been made to 
modify the turbulence model, the present results are very encouraging. A sensitl- 
vlty study will be conducted to determine the effect of changes In transition loca- 
tion and turbulence model on the St number distribution in the near future. 

It Is planned to use the "TABLET" code for predicting the viscous flow on the 
suction surface of the blade to study the effect of streamline convergence on heat 
transfer. Finally, this code will be used to predict the skewing of surface 
streamlines on the pressure surface of the rotating turbine blade of Ref. 3, to 
complete the assessment of the applicability of the 3-DBL theory for analyzing 
viscous flow in a turbine passage. 
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Fig. 5. Schematic of Experimental Setup 
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Fig. 6. Wall Cross Flow Angle Distribution 
for Flow in a Curved Duct 



Fig. 7a. Schematic Diagram of Endwall Fig. 7b. Comparison of Endwall Stanton 
Computational Region Number Distributions 
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GAS-SIDE HEAT TRANSFER WITH ROTATION 


Robert P. Dring 

United Technologies Research Center 

The primary basis for heat transfer analysis of turbine blades is experimental 
data obtained in linear cascades. These data have been very valuable in identifying 
the major heat transfer and fluid flow features of a turbine airfoil. The question 
of major interest is how well all of these data translate to the rotating turbine 
blade. It is known from the work of Lokay and Trushin (ref. 1) that average heat 
transfer coefficients on the rotor may be as much as 40 percent above the values 
measured on the same blades non-rotating. Recent work by Dunn and Holt (ref. 2.) sup- 
ports the Russian conclusion. What is lacking is a set of data from a rotating 
system which is of sufficient detail as to make careful local comparisons between 
static cascade and rotor blade heat transfer. In addition, data is needed in a ro- 
tating system in which there is sufficient documentation of the flow field to support 
the computer analyses being developed today. A second major question is the influ- 
ence, if any, of the first stator row on the heat transfer of the second stator row 
after the flow has passed through the rotor. 

The first objective of the present program, is to obtain a detailed set of heat 
transfer coefficients along the midspan of a blade in a rotating turbine. The tur- 
bine stage to be examined is that of Dring et al. (ref. 3), as shown in figure 1. 
Aerodynamic testing of this stage has already been conducted to document, among other 
things, the impact of rotor-stator interaction on the airfoil unsteady pressure dis- 
tributions. The time averaged and unsteady pressure envelopes are shown in figure 2. 
The rotor heat transfer data are to be such that they can be compared directly with 
data taken in the static cascade experiment of Graziani et al. (ref. 4), as illus- 
trated in figure 3. The data are also to be compared to some standard analysis of 
blade boundary layer heat transfer which is in use today. In addition to providing 
this all-important comparison between rotating and stationary data, this experiment 
should provide important insight to the more elaborate full three-dimensional pro- 
grams being proposed for future research. A second program objective is to obtain 
a detailed set of heat transfer coefficients along the midspan of a stator vane loca- 
ted in the wake of a full upstream turbine stage. Particular focus here is on the 
relative circumferential location of the first and second stators. The final program 
objective is to improve the analytical capability to predict the experimental data. 

The program is divided into two phases. The first phase is designed to provide 
a high density of high accuracy heat flux measurements on both the pressure and 
suction surfaces of a rotating turbine blade, figure 1. The program will include two 
rotor-stator spacings to assess the influence of rotor-stator interaction as in the 
experiment of Dring et al, (ref. 3). It will also include both positive and negative 
incidence as well as the design incidence. Finally, the first phase program will 
include comparison of the data with static cascade data, figure 3, and with a common 
boundary layer heat transfer analysis. The second phase will be designed to provide 
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a high density of high accuracy heat flux measurements on both the pressure and suc- 
tion surfaces of a second stator row turbine vane. The program will be run at one of 
the rotor-stator spacings used in the first phase. It will include, at minimum, four 
different circumferential positionings of the first stator row relative to the second 
stator row. In addition, the second phase will include an analysis to describe the 
experimental phenomena of both phases, which will improve the current ability to pre- 
dict turbine midspan airfoil heat transfer. 
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OBJECTIVES 

• Turbine rotor heat transfer data 
e Rotor-stator interaction {AX} 

e Comparison with cascade data 
e Second stator heat transfer data 

• Stator-stator interaction (ACirc) 

• Quasi-steady and/or unsteady 
boundary layer analysis 
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FIGURE 1. TURBINE STAGE AT 15% GAP (Ref. 3) 


15% gap, <C X /U) = 0.78 



FIGURE 2. UNSTEADY PRESSURE ENVELOPES (Ref. 3) 



FIGURE 3. MEASURED AND PREDICTED MIDSPAN HEAT 
TRANSFER DISTRIBUTIONS (Ref. 4) 
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COOLANT PASSAGE HEAT TRANSFER WITH ROTATION 
A Progress Report on the Computational Aspects 


L. D. Aceto, Program Manager 
G. J. Sturgess, Presenter 

Pratt & Whitney Engineering, Pratt & Whitney Aircraft Group 


Turbine airfoils are subjected to increasingly higher heat loads which escalate 
the cooling requirements in order to satisfy life goals for the component materials. 
If turbine efficiency is to be maintained, however, cooling requirements should be as 
low as possible. To keep the quantity of cooling air bounded, a more efficient 

internal cooling scheme must be developed. One approach is to employ airfoils with 
multi-pass cooling passages that contain devices to augment internal heat transfer 
while limiting pressure drop. (fig. 1). 

Design experience with multi-pass cooling passage airfoils has shown that a 
surplus of cooling air must be provided as a margin of safety. This increased cooling 
air leads to a performance penalty. Reliable methods for predicting the internal 
thermal and aerodynamic performance of multi-pass cooling passage airfoils would 
reduce or eliminate the need for the safety margin of surplus cooling air. 

The objective of the program is to develop and verify improved analytical methods 
that will form the basis for design technology which will result in efficient turbine 
components with improved durability without sacrificing performance. The objective 
will be met by: 1) establishing a comprehensive experimental database that can form 
the basis of an empirical design system. 2) developing computational fluid dynamic 
techniques, and 3) analyzing the information in the database with both 
phenomenological modeling and mathematical modeling to derive a suitable design and 
analysis procedure. 

Currently, a design phase has been completed which defined a rotating experiment 
to simulate blade passage cooling. Models of blade internal geometries will be 
studied. The first model with instrumented smooth wall legs has been designed and is 
being fabricated. It is scheduled for testing in December, 1983. The subsequent 
models will have "rough" walls and will address passage aspect ratio changes. These 
tests will generate a database for the development of a design system. 

The prediction of local coolant side heat transfer coefficients and coolant 
temperature rise and pressure drop in the cooling passages is difficult because of 
the three dimensional, elliptic nature of the flow, the effects of rotation, and the 
complex turbulence promoting devices. The analytical approach under consideration is 
based on the 3D-TEACH computer code being developed for gas turbine combustor 
applications by Pratt & Whitney Aircraft, (ref. 1). 

The 3D-TEACH code is a generalized aerothermal fluid dynamic solver for 
three-dimensional, elliptic, turbulent, steady flows. The approach stays within the 
framework of continuum mechanics and uses a statistical description of turbulence, 
coupled with the accepted Eulerian description provided by the Navier-Strokes 
equations of motion. Closure to the resulting time-mean equations is provided by 
turbulence modeling of the eddy viscosity type. The modeled partial differential 




equations are manipulated into a general form that permits a single solution 
algorithm to be used for a numerical procedure. A hybird (upwind /central) finite 
differencing scheme is used to discretize the equations. 

A version of the 3D-TEACH code has been modified to calculate fluid flow and heat 
transfer in rotating passages. In order to be able to apply the code with confidence 
to the experimental data generated in this program, it is essential to assess the 
suitability of the physical modeling used. To do this, the code is being exercised in 
flow situations similar to the turbine airfoil application, but simpler in nature. 
Such a procedure enables the weaknesses, if any, to be identified easily. 

To test the behavior of the modified momentum equations flow in a rotating 
multi-pass passage (figure 1), was calculated. Flow patterns are made visible by the 
use of streak lines. Figure 2 shows flow at 600 RPM (0*174 Rossby Number) in two 
cross-sections of the passage. The development of a vortex pair in the outflow (away 
from axis of rotation) leg due to Coriolis forces can be seen clearly. The action of 
the sharp turns is to coalesce these vortices into a single vortex that dominates 
flow in the inflow leg. The observed behavior appears sensible. 

The quantitative behavior in rotation was examined through reference 2. For 

example, a comparison of velocity profiles for the duct rotating at 165 R.P.M. is 
given in figure 3. Streakline flow visualization shows also the development of 

Coriolis vortices in the duct. 

The ability of the code to calculate heat transfer in a rectangular duct with a 
sharp 180 degree bend was examined through reference 3. A comparison of calculated 
and measured Nusselt numbers for the various zones through the passage is given in 
figure 4 at two Reynolds numbers. 

The behavior of the code thus far is pleasing, and the results of the comparisons 
with the verification experiments is very encouraging. Further investigations of 
these and other suitable experiments will continue to build confidence in the code. 
Use of more sophisticated turbulence modeling will also be explored. The developed 
code will then be applied to the measurements to be made as part of this program. 
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Figure 4 Heat Transfer In 180° Sharp Bend For Metzger's Experiment 
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JET ARRAY IMPINGEMENT HEAT TRANSFER CHARACTERISTICS 


L.W. Florschuetz and D.E. Metzger 
Department of Mechanical and Aerospace Engineering 
Arizona State University 

Two-dimensional arrays of circular air jets impinging on a heat transfer surface 
parallel to the jet orifice plate are considered. The jet flow, after impingement, 
is constrained to exit in a single direction along the channel formed by the jet 
orifice plate and the heat transfer surface. In addition to the crossflow which 
originates from the jets following impingement, an initial crossflow is present 
which approaches the array through an upstream extension of the channel (fig. 1) . 
The configurations considered are intended to model the impingement cooled midchord 
region of gas turbine airfoils in cases where an initial crossflow is also present 
(fig. 2). A major objective is determination of the effect of initial crossflow air 
temperature relative to jet array air temperature on impingement surface heat 
fluxes . 

Earlier work in this NASA sponsored project was directed at modeling uniform 
arrays in cases where an initial crossflow is not present. In those cases, referred 
to here as noninitial crossflow geometries, there was an endwall or upstream edge to 
the channel positioned one-half a streamwise hole spacing upstream of the first 
spanwise row of holes in the array. Note that with the initial crossflow geometry 
(fig. 1), the special case of zero initial crossflow (m c = 0) is identical to the 
noninitial crossflow case except for the presence of the upstream channel endwall in 
the latter case. 

In this extended abstract, nomenclature and definition of parameters is first 
indicated for the more general initial crossflow case. Then the main features of 
the prior noninitial crossflow studies are briefly summarized. The current status 
of the initial crossflow work is then discussed. 


Consider a uniform rectangular array of circular jet orifices with an inline 
hole pattern (figs. 1 and 3). Consider further the heat flux, q, at the impingement 
surface, averaged across the span, but resolved in the streamwise direction to a 
region of width x n , centered opposite an arbitrary spanwise row of holes within the 
array at location x. It is desired to express q as a function of parameters 
associated with the entire array. Assume constant fluid properties, negligible 
viscous dissipation, a uniform impingement surface temperature, T s , an adiabatic jet 
orifice plate surface at the jet exit plane, and a fully developed channel flow at 
the entrance to the array (x = 0) with a mixed-mean temperature T c . Then, at the 
specified x 

q — f (mj, m c , Tj , Tq, Tg, x^, y^» z, d, L, p., k, Cp) 
flow temperature geometry fluid 

or in one possible dimensionless form, at x/L, 

= f[m c /mj, Rej , (T c -Tj) /(T s -Tj) , x n /d,y n /d, z/d, L/x n , Pr] 


k(Ts-Tj) 
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In this form the dependent parameter, which has the form of a Nusselt number, 
incorporates a heat transfer coefficient defined as q/(T s -Tj) with the jet 
temperature Tj as the characteristic fluid temperature. With this definition the 

Nusselt number is not independent of the temperature differences of the problem, 
since it still depends on the initial crossflow temperature ratio (T c -Tj ) / (T s -Tj ) . 
In spite of the restrictions imposed, there are still eight independent 
dimensionless parameters, nine including the specified x/L. Note that one of these 
is L/x n which is just the number of spanwise rows in the array. 

Still considering q in terms of parameters associated with the entire array, an 
alternate formulation is 

q - f (mj , m c , Taw 1 T a , ^n* yn j 1, ep) 

where 

Taw = ^®j* ®c * Tj , Tq , x^, y n , z, d, L, g, k, Cp) 

Then in dimensionless form 

Nu = ^ _ = f(m 0 /mj. Re j , x n /d, y n /d, z/d, L/x n , Pr) (1) 

k(Tg-T aw ) 

and 

q = Taw = f(m c /mj, Rej, x n /d, y n /d, z/d, L/x n , Pr) (2) 

T c-Tj 

This formulation, employing the adiabatic wall temperature, has the advantage that 
the heat transfer coefficient, q/(T s -T aw ), and Nusselt number are rendered 
independent of the temperature differences of the problem. But, unless T c =Tj , the 
additional dependent parameter q is also needed. This Nusselt number and the 
corresponding q are a function of seven independent dimensionless parameters, eight 
including the specified x/L. 

At this point the possibility of representing q as a function of parameters 
associated only with a given spanwise row is considered (fig. 4). Here, it is 
assumed that the dependence of q on the velocity and temperature distribution at the 
control surface (channel cross-section) immediately upstream of the row may be 
represented in terms of a single velocity and a single temperature parameter, 
denoted by G c and T c . Then in dimensionless form one may write 

Nu = q d = f(G c /Gj, Rej, x n /d, y n /d, z/d, Pr) (3) 

k(T s -T aw ) 


and 


Taw = f(G c /Gj, Rej, x n /d, y n /d, z/d, Pr) (4) 

T c -Tj 


In general, of course, q will depend on the velocity and temperature distributions 
over such a control surface, which, in turn, depend on the history of the mixing jet 
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flows and crossflows upstream of the control surface. However, the above 
approximation, if shown to be adequate , admits the possibility of applying 
measurements obtained for uniform array geometries to heat transfer opposite 
individual rows of nonuniform array geometries. 

Results based on an extensive series of tests for a range of uniform array 
geometries in noninitial crossflow configurations using air as the working fluid are 
reported in tabular and graphical forms in references 1, 2 and 3 in terms of the 
parameters of equation (1), except for m c /mj which is not a relevant parameter in 
noninitial crossflow cases. Results were obtained for every geometry tested with Nu 
resolved in the streamwise direction to at least x n , and in most cases to x n /3. 
Correlations, in algebraic form, for Nu (resolved to x n ) expressed in terms of 
parameters associated with a single spanwise row as in equation (3), are reported in 
references 4 and 5. In these references, a theoretically based flow distribution 
model is also presented and validated by comparison with measured flow distributions 
for these noninitial crossflow cases. The validated flow model was used in 
developing the correlations referred to above. Results for nonuniform arrays in 
noninitial crossflow configurations have also been obtained and compared with the 
uniform array data and the correlation based on that data (ref. 6) . Examples are 
shown in figure 5. 

One of the more extensive prior heat transfer studies of two-dimensional arrays 
of circular jets was reported in reference 7. The studies completed under the 
present project go substantially beyond prior work, and provide a much more complete 
understanding of the problem as well as needed design information. The main 
features of the present project which had not been addressed in prior work may be 
summarized as follows: (1) streamwise spatial resolution of at least x n for all 
geometries tested; (2) observations of ''damped'' periodic streamwise variations of 
heat transfer coefficients along entire array with spatial resolution of x n /3j 
(3) detailed verification of validity of geometric scaling; (4) rectangular arrays, 
i.e., not restricted to Xn/yn = 1; (5) staggered as well as inline hole patterns; 

(6) correlation for streamwise local application developed in algebraic form; 

(7) nonuniform as well as uniform array geometries tested; (8) experimentally 
validated theoretical flow distribution model developed for both uniform and 
nonuniform array geometries; (9) applicability of uniform array data and correlation 
to nonuniform arrays examined; (10) effect of confined crossflow on jet orifice 
discharge coefficients; (11) effects of initial crossflow rate and temperature, 
currently in progress, discussed further below. 

An example of heat transfer characteristics with initial crossflow, specified in 
terms of array parameters as in equations (1) and (2), is shown in figure 6 for the 
B(5,4,3)I geometry (i.e., B-size plenum, x n /d=5, y n / d=4 , z/d=3. Inline hole 
pattern). The corresponding flow distribution is shown in figure 7. Additional 
such results are reported in reference 6. 

Nusselt numbers specified in terms of parameters at an individual spanwise row, 
as in equation (3), are illustrated in figure 8. The presence of initial crossflow 
extends the data point range of G c /Gj beyond the maximum values existing for the 
noninitial crossflow (or essentially identical zero initial crossflow) cases. The 
B(5,4,2)I geometry (upper plot in fig. 8) clearly shows that a minimum value of Nu 
may occur at a certain G c /Gj. A point is reached where the jet no longer 
effectively impinges on the surface, but the Nusselt number increases with the 
increasing crossflow velocity. The prior noninitial crossflow correlation (refs. 4 
and 5) though shown extrapolated beyond the range of G c /Gj on which it was based. 
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appears to be reasonably consistent with the data points until the increase in 
Nusselt number begins, except for several data points from the first and second 
spanwise rows of the array. 

Dimensionless adiabatic wall temperatures specified in terms of parameters at an 
individual spanwise row, as in equation (4), are illustrated in figures 9 and 10. 
In figure 9, T m>11 , the mixed-mean total temperature, was used as the characteristic 
crossflow temperature, T c , immediately upstream of a given row. Figure 10 shows 
data for the same cases, but with T aw >n -i, the adiabatic wall temperature opposite 
the row immediately upstream of the row in question, used to represent T c . It is 
not clear at this writing whether either of these approaches will prove to be of 
general utility over a range of geometries. This question is currently being 
studied. In some instances the dimensionless adiabatic wall temperature at the 
first row or two, like the Nusselt number, clearly does not follow the general trend 
of the data points. This is because the relationship between q and T s at the 
impingement surface opposite a given row clearly must depend on the details of the 
velocity and temperature distributions over the channel cross-section immediately 
upstream of the row. These distributions at downstream rows may differ considerably 
from those at the first several rows of an initial crossflow array, and their 
effects cannot necessarily be represented accurately by a single crossflow velocity 
and a single temperature parameter. 
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NOMENCLATURE (includes symbols not explicitly defined in text or figures) 


G c 

k 

m c 

m j 

Pr 

<1 

Re 4 


l aw 



jet plate discharge coefficient 
constant pressure specific heat 
friction coefficient 

crossflow mass velocity based on channel cross-sectional area 
jet mass velocity based on jet hole area 

heat transfer coefficient at impingement surface, q/(T s -T aw ) 

thermal conductivity 

initial crossflow rate 

total jet flow rate 

dynamic viscosity 

Prandtl number 

heat flux at impingement surface 
jet Reynolds number, Gjd/p 

adiabatic wall temperature (with subscript n, denotes location opposite 
spanwise row number n) 

characteristic temperature of initial crossflow or of crossflow within array 
depending on context 

characteristic temperature of jet flow 

crossflow channel mixed-mean temperature upstream of spanwise row n 
heat transfer surface temperature 


Superscript 

( — ) = overbar refers to mean value over jet plate 
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T i ’ m j 


JET PLATE 




Figure 2. Impingement cooled airfoil - midchord jet arrays 
subject to initial crossflow. 
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CONSIDER AT SPECIFIED x 

q= f(rry m c , T., T c> T s , x n , y n , z, d, L, /u., k, c pi ) 

FLOW TEMPERATURE GEOMETRY FLUID 

THEN POSSIBLE DIMENSIONLESS FORM AT SPECIFIED x/L IS 


q d 

kCTs-Tj) 


f [m c /m., Re,, (T C -T.>/(T S -T.), x n /d, y n /d, z/d, L/x n , Pr] 


Figure 3. Possible formulation of streamwise resolved impingement 
surface heat transfer characteristics in terms of array parameters. 



CONSIDER 

q= f ( G c » G j » T aw’ Ts> x n’ ^n’ z >^>F t, ^> c p^ 

WHERE 

^aw = * ^j* V x rr ^n’ z> c p^ 

DIMENSIONLESS FORM 


q d 

k ( T $ - T QW ) 




Figure 4. Possible formulation of streamwise resolved impingement surface 
heat transfer characteristics in terms of individual spanwise jet row parameters. 
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B(5, 8,3) I -*t*-B(5,4,3)I 


A Non-uniform geometry data 

V Uniform geometry data 
(Spanwise row numbers 
noted over abscissa) 

— Uniform geometry correlation 

t i i i t 


Rej = 1.60 X 10 


* x 


A Non-uniform geometry data 

V Uniform geometry data 
(Spanwise row numbers 
noted over abscissa) 

—— Uniform geometry correlation 


B(5,8,3)I^j— B(5,4,3)I 


I * t 
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Figure 5, Examples of nonuniform array geometry Mussel t numbers 
compared ifith. uniform array data and correlation. 
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Nu=hd/k -q = (T aw -Tj)/(T C -Tj) 
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ignre 6 » Examples of streamwise profiles of 
initial crossflow Nnsselt numbers (Nu) and 
dimensionless adiabatic wall temperatures (tj) 
specified in terms of array parameters. 


Figure 7. Examples of flow distributions with 
initial crossflow for array geometry of figure 6 










Figure 8. Examples of stream-wise resolved Nusselt numbers as a function of 
individual spanwise row cross-to-jet mass velocity ratios. Solid symbols denote 
data from first and second upstream rows of array. 
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8(5, 4,3) I 

m c /m j 
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□ 0.47 
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m c /m j 
A 0.20 
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Figure 9. Examples of dimensionless adiabatic wall temperature as a function 
of individual spanwise row cross-to-jet mass velocity ratio, using mixed-mean 
total temperature upstream of row as characteristic crossflow temperature. 
Solid symbols denote data from first and second rows of array. 
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igure 10, Examples of dimensionless adiabatic wall temperature as a function 
of individual spanwise row cross-to-jet velocity ratio, using adiabatic wall 
temperature at upstream row as characteristic crossflow temperature. 

Solid symbols denote data from first and second rows of array. 






COMBUSTION HOT SECTION TECHNOLOGY 



David B. Ercegovic 


Propulsion Laboratory 

AVARDCOM Research and Technology Laboratories 
Lewis Research Center 
Cleveland, Ohio 


The overall objective of the Turbine Engine Hot Section Technology Combustion 
Project is to develop and verify improved and more accurate analysis methods 
for increasing the ability to design with confidence the combustion system 
for advanced aircraft turbine engines. The analysis methods developed will 
be generically applicable to combustion systems and not restricted to one 
specific engine or manufacturer. 

This project's approach was to first assess and evaluate existing combustor 
aerothermal analysis models by means of a contracted effort initiated during 
FY '82. This evaluation effort has assessed and quantified known models' 
strengths and deficiencies. The results of this assessment will be summarized 
later in this report by the three contractors Garrett, General Electric and 
Pratt and Whitney Aircraft who conducted this independent evaluation. A 
balanced contract and in-house program will then be conducted to support, 
focus, and accelerate the development of new methods to more accurately pre- 
dict and model the physical phenomena occurring within the combustor. This 
program will include both analytical and experimental research efforts in 
the areas of aerothermal modeling and liner cyclic life. 

During FY '84 the Aerothermal Modeling Program - Phase II will be initiated, 
which is expected to have contracted model development efforts in the areas 
of improved numerical methods for turbulent viscous flows, flow interactions, 
and fuel spray flow-field interactions. A Phase III effort is planned to add- 
ress remaining model deficiencies. The primary in-house effort in this area 
will be the determination of high pressure flame radiation characteristics in 
a full annular combustor. This experiment will be conducted in the NASA LeRC 
High Pressure Facility with the results compiled into a comprehensive flame 
radiation and liner heat flux model. 

In the area of liner cyclic life, HOST is engaged in a co-operative effort 
with Pra-t and Whitney Aircraft to develop a test apparatus to economically 
determine combustor thermal strains and cyclic life. This test apparatus 
will be run in-house at NASA LeRC and will be the test vehicle for many of 
the advanced high temperature instruments developed under HOST sponsorship. 

The fundamental data generated in this project will be used to assess and 
further develop current analytical liner life programs. 
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AEROTHERMAL MODELING PROGRAM - PHASE I 


R. Srinivasan, R. Reynolds, K. Johnson, and H. Mongia 
Garrett Turbine Engine Company 


The main objective of the NASA-sponsored Aerothermal Modeling Pro- 
gram, Phase I, was to assess current aerothermal submodels used in the 
Garrett Turbine Engine Company (GTEC) analytical combustor models. 

A number of "benchmark" quality test cases were selected after an 
extensive literature survey (ref. 1). The selected test cases, includ- 
ing both nonreacting and reacting flows, were broadly divided into the 
following categories: 


o 

Simple flows 


o 

Complex nonswirling 

flows 

o 

Swirling flows 


o 

Dilution jet mixing 

in confined 


crossflows . 


These test cases were used to assess the following submodels sepa- 
rately and jointly for various combustion processes: 


o k-€ model of turbulence and algebraic stress model, with and 
without various corrections, including low Reynolds number and 
Richardson number corrections 

o Scalar transport models 

o Multistep kinetic schemes 

o Turbulence/chemistry interaction 

o Spray combustion. 


Based upon the Phase I work, the following assessments of the sub- 
models were made: 


o The k-€ model and its modifications give good correlations for 
Simple flows 

Far-field regimes of nonswirling/swirling flows involv- 
ing regions of recirculation 

Non-recirculating swirling flows 

Outer regions of strong swirling flows. 
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o These models give reasonable correlations for 

Nonswirling recirculating flows excluding the vicinity 
of any attachment point 

Confined disk flow with a central jet 

Shear layer of strong swirling flows 

Confined swirler with hub and shroud expansions. 

The k-€ model and its modifications predict only the trends in 
recirculation zones of swirling flows and confined swirler flows without 
outer expansion. 

The Reynolds stresses predicted by the algebraic stress model (ASM) 
correlate well for simple flows. The ASM correctly predicts the normal 
stresses in nonswirling recirculating flows. The ASM predictions give 
reasonable correlations for: 

o Shear stresses in nonswirling recirculating flows 

o Normal stresses in swirling flows. 

For predicting the scalar transport, the k-6 model with constant 
Prandtl number gives good correlation where the gradient diffusion 
approximation is valid. The algebraic scalar transport model is a prom- 
ising approach and more validation efforts are needed. 

On the turbulence/chemistry interaction models, the following gen- 
eral conclusions were made: 

o Both two-step and four-step kinetic schemes show promise 

o A modified eddy breakup model can easily be extended to multi- 
step kinetic schemes 

o Bilger's two-reaction- zone model gives good results for jet 
flames, but requires more work. 

When dilution jet mixing was considered, the three-dimensional ana- 
lytical model had the following characteristics: 

o Slightly underpredicted jet penetration 

o Centerline temperatures predicted well 

o Transverse mixing predictions slower than data 

o Effect of S/D, H/D, J on mixing predicted correctly (where D = 
geometric orifice diameter, H = local duct height at the sur- 
vey plane, J = jet versus mainstream momentum ratio, S = ori- 
fice spacing) 
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o Good correlation for jet injection from 
One wall 

Both walls - inline orifices 
Both walls - staggered orifices 

The following general conclusions were derived from the Phase I 
work (ref . 2) : 

o An accurate numerical scheme should be developed to minimize 
numerical diffusion in the computations of recirculating flows 

o Benchmark quality data should be generated under well defined 
environments for validating the various submodels used in gas 
turbine combustion analysis 

o Although current aerothermal models make reasonable predic- 
tions, intensive model development and validation efforts 
should continue for the following submodels: 

Algebraic stress model 

Algebraic scalar transport model 

Two-step and four-step schemes 

Probability density function approach for a two-step 
scheme 

Double-reaction-zone model. 
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A Closing Progress Report 


G. J. Sturgess 
Pratt & Whitney Engineering 


The overall objective of the two-phase Aerothermal Modeling Program is to 
develop the computational fluid dynamics tools needed to improve the design, analysis 
and development process for the gas turbine engine combustion chamber- In the first 
phase, 1982 state-of-the-art models were evaluated, shortcomings were identified, and 
improvements to these models were recommended. The models were of steady-state, 
time-averaged formulation. It is planned to implement the recommendations made in the 
first phase as part of the second phase effort. 

The results of the study highlighted the importance of inlet boundary conditions 
when making evaluations of the modeling by means of comparison with benchmark-quality 
experimental data (ref. 1). The study also revealed that the calculation accuracy was 
not independent of the numerics used in the solution procedure. It proved necessary 
to establish a hierarchy of physical models (fig. 1), whose realism must be improved 
in sequence. 

It was found that relatively simple flows could be calculated with acceptable 
accuracy (fig. 2). However, there was less success with more complicated flows (fig. 
3). Shortcomings were identified in the turbulence model with respect to swirling 
flows and flows with strong streamline curvature, and for calculation of turbulent 
mass transport (table 1) . 

The evaluation suffered from a serious lack of readily available and suitable 
benchmark quality experiments. Many experiments exist in the literature, but they 
frequently suffer from incomplete boundary conditions, lack of parametric variation 
of operating conditions, and clear delineation between the physical processes 
involved. 

It was concluded that the current models are of qualitative accuracy only, and 
are therefore suitable for scale-up, perturbation, and diagnostic studies. In these 
roles they can serve as useful engineering tools. Further work is required to improve 
quantitative accuracy for the intended use. 

Phase I has proved useful in that it has provided a detailed and systematic 
examination of the physical modeling used, and it has shown the weaknesses inherent 
in this modeling. It has also provided substantial insight into the powerful 
influences exerted on solution accuracy by the inlet boundary conditions, and the 
numerics used in the computer codes. 

The study made it apparent that serious difficulties exist with the present 
generation of computational fluid dynamics codes for application to the combustor. 
These difficulties are not perceived as being fatal however. The economic and 
competitive driving forces in the industry that gave rise to the codes are still real 
and urgent, and the codes have indeed demonstrated great potential for productivity 
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improvement in addition to great technical capability. Ways to improve the existing 
methods exist, and others can be devised. There is no reason or need to abandon the 
effort. Improvement of the codes should be continued. 
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o • □ ■ Experimental data 


Figure 2 

MORE COMPLEX FLOWS NOT 
CALCULATED AS WELL 

Honeycomb 



Reacting flow: 2 kg/sec. air, 16 kg/hr. propane 
Figure 3 
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Table 1 


USE OF CONSTANT TURBULENT SCHMIDT 
NUMBER APPEARS INAPPROPRIATE 


Flow type 

Characteristics 

Schmidt 

No. 

Boundary layers 

Predominantly unidirectional 

0.9 

Round free jets 

2-D, external entrainment 

0.7 

Simple recirculation zones 

Simple stagnation points 

0.5 

Confined round jets 

2-D, self-entrainment 

0.7 

Complex recirculation zones 

2 or 3-D, multiple stagnation points, 
multiple recirculations 

0.2 
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AEROTHERMAL MODELING PROGRAM 


PHASE I 

GENERAL ELECTRIC COMPANY 
M. J. Kenworthy 


The objective of this Phase I Aerothermal Modeling program was to identify 
deficiencies in the aerothermal models that are currently used or planned for 
use at General Electric in the design of gas turbine engine combustors. This 
effort was accomplished by assessing the predictions of the existing models 
with "benchmark" quaility test data. This effort, together with other similar 
efforts sponsored by NASA, was intended to help assess the current state of 
the art, and to define Phase II efforts needed to evolve improved and more 
advanced aerothermal analysis methods, increasing the ability to design 
improved high performance, durable combustors for advanced aircraft gas 
turbine engines. 

This Phase I effort focused on the ability of the existing models to 
predict or analyze those features associated with durability. These primarily 
are; 1) the combustor liner metal temperature distribution, and 2) the 
combustor exit gas temperature patterns, which in turn affects the life of the 
downstream turbine components. Other features of interest to the combustor 
designer such as blowout limits and emissions do not directly affect design 
for durability. 

Figure 1 shows schematically General Electric's calculation modules 
within the overall aerothermal model. Except for the internal flow module, 
the modules are well developed and have been in use for many years. The 
internal flow module is planned for use in the design process at General Electric 
after it has been developed and/or demonstrated to have useful accuracy. Thus, 
the assessment of the existing internal flow module constituted a major 
portion of the efforts. 

The internal flow module, INTFLOW, contains two basic computing codes: 

1) the code assembled by the Garrett Turbine Company under the Army sponsored 
Combustor Design Validation Program efforts, Reference 1, and the very similar 
code assembled by the Northern Research and Engineering Company, Reference 2. 

To help evaluate accurately the physical submodules and the numerical tech- 
niques within INTFLOW or under consideration for eventual incorporation in 
INTFLOW, the two-dimensional (planar or axisymmetric) research code, GETREF 
(general JElectric JTurbulent Reacting _Flow) developed at General Electric was 
also utilized in the assessments. 

Four axisymmetric benchmark experiments were chosen to assess the 
capabilities of available two-dimensional analysis methods. GETREC was used 
in conducting these assessments. This first experiment which was performed 
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at General Electric, involved detailed studies of an axisymmetric fuel jet 
burning in a coflowing air (Reference 3) . The primary reasons for its 
selection are: 1) the removal of "numerical diffusion" errors since the 

governing equations become parabolic, and 2) the high quality LDV and 
spontaneous Raman scattering data available. An unambiguous evaluation of 
the turbulence chemistry interaction model was made which confirmed the model 
approach. In addition to this jet experiment, three other axisymmetric 
recirculating flow experiments were chosen: an isothermal annular dilution 

jet experiment (Reference 4), an isothermal double swirled pipe flow 
(Reference 5), and a bluff-body stabilized non-premixed flame (Reference 6). 
These provided checks of the model under progressively more complex circum- 
stances. The results demonstrated the value of an algebraic stress model 
for modifying the K-e turbulence model in the presence of swirl, and showed 
that the recirculation zone with heat release could be closely predicted. 

For assessments of three-dimensional analysis capabilities, the aero- 
dynamics of simple round dilution jet flows in a crossflow were studied. Data 
from Reference 7 were chosen for this examination. To provide data from more 
complex flowfields, a series of experiments were devised and conducted as 
part of this aerothermal modeling program which involved progressively more 
complex flows beginning with a row of jets in a crossflow followed by: 
opposing jets, alternate jets of different sizes, and then a second row of 
dilution jets with the two row pattern simulating the dilution pattern of one 
General Electric combustor. This was followed with experiments including 
swirl cups in the dome and heat release. 

Assessments with INTFLOW began with variations of grid density and 
turbulence input parameters to examine agreement with the dilution jet in 
crossflow. Figure 2 shows the most detailed of the grids examined. Figure 3 
compares measurements of jet penetration and temperature mixing with calcu- 
lations achieved with the two basic codes within INTFLOW using identical input 
parameters. It is significant that this excellent agreement was achieved 
even though the two codes calculated different turbulence intensity levels 
throughout the flowfield. Also, difficulties were encountered in predicting 
all jet momentum ratios well with the same turbulence input parameters. 

Figure 4 shows measurements from a F101 combustor sector compared with 
calculations. The overall pattern is reasonably predicted showing hottest 
regions concentrated between fuel nozzles. This calculation, utilized as 
input, measured flowfield conditions just downstream of the dome swirl cup. 

While plausible looking results can be achieved with the INTFLOW module 
and agreement can be achieved after some educated selection of the inputs, 
predictions sufficiently accurate for aiding in design or development for 
enhancing durability have not yet been achieved within the limits of available 
(or even reasonable) computer storage. Improved numerics, reduced computation 
time, and better input information/ techniques for the swirl cup fuel injector 
region are needed in addition to the improved physical modeling features 
demonstrated during this program. 
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Figure 1. General Electric’s Aerothermal Model. 


141 









ANNULUS HEIGHT 


Y - Z PLANE 



Figure 2. 


10 cm H 

X - Z PLANE 


High Grid Density, 41,600 Nodes, for Dilution Jet Study. 


TEST DATA 


NEEC CALCULATION 


G AR RETT C ALC U LATION 





TANGENTIAL DIRECTION 


Figure 3. Comparison of Jet Penetration Temperature Contours cf Two 

Basic Codes within INTFLOW with Measurements at 4 X/D Down- 
Stream of Jet. 


142 







Page intentionally left blank 


DILUTION JET MIXING PROGRAM - PHASE II 


R. Srinivasan, E. Coleman, K. Johnson, and H. Mongia 
Garrett Turbine Engine Company 


This program is a continuation of the ongoing program under NASA 
Contract NAS3-22110. In Phase I of this program, the dilution jet mix- 
ing study was limited to single-sided jet injections (ref. 1) . The 
objectives of Phase II of the dilution jet mixing program are; 

o Extend the data base on mixing a single-sided row of jets with 
a confined cross flow 

o Collect data base on mixing a two-sided row of jets with a 
confined cross flow 

o Develop empirical jet mixing correlations 

The parameters investigated during Phase II of the program are; 

o Circular versus square orifices, two-dimensional slot 

o J, H/D, S/D (where D = geometric orifice diameter, H = local 
duct height at the survey plane, J = jet versus mainstream 
momentum ratio, S = orifice spacing) 

o In-line and staggered orifice configurations 

o Nonuniform cross-stream temperature and velocity profiles 

o Cross-stream flow area convergence 

Figure 1 shows a flow schematic of the dilution jet test rig. The 
mainstream temperature and velocity profiles can be tailored by adjust- 
ing the flow to the profile generator upstream of the test section. 
Dilution air enters the test section through sharp-edged orifices in the 
test section walls from the top, bottom, or both. 

Figure 2 provides more detail on the test sections and orifice con- 
figurations used in this study. The height of the test section at the 
injection plane, Ho, was 10.16 cm for all tests. Orifice plate open 
areas were varied from 2.5 to 20 percent of the mainstream cross section 
at the injection location. The primary independent geometric variables 
are the orifice size and the spacing between adjacent orifices. These 
are conveniently expressed in dimensionless form as the ratio of the 
duct height to orifice diameter, Ho/D, and the ratio of the orifice 
spacing to duct height, S/Ho. The product of these ratios is the orifice 
spacing-to-diameter ratio, S/D. 

Tests were performed with two-sided injection in convergent ducts 
having a convergence rate of 0.5 cm/cm. The tests were made in symmetric 
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as well as asymmetric convergent ducts with isothermal mainstream condi- 
tions . 


The dilution jet mixing characteristics were determined by measur- 
ing temperature and pressure distributions with a vertical rake probe, 
positioned at different axial and lateral stations. This probe had 20 
thermocouple elements, with a 20-element total pressure rake, and a 20- 
element static pressure rake located nominally 5 mm (0.05 Ho) on each 
side of the thermocouple rake. The center-to-center spacing between 
sensors on each rake was also 0.05 Ho. 

The vertical rake probe was traversed over a matrix of 44 to 60 Z-X 
plane survey locations (where X is parallel to the duct axis? Y is paral- 
lel to the orifice centerline, radial direction; and Z is normal to the 
duct axis, transverse direction) . The flow-field mapping in the Z-direc- 
tion was done over a distance equal to one or one and a half times the 
hole spacing, S, at intervals of S/10. For most tegts, the X-Y plane 
containing the orifice centerline (centerplane) was at the center of the 
span surveyed; i.e., data surveys were from midplane to midplane. 

Measurements in the X-direction were made at up to 5 planes with 

0. 25. < X/Ho<2. Note that because the designers' objective was to iden- 
tify dilution zone configurations to provide a desired mixing pattern 
within a given combustor length, the downstream stations of interest 
were defined in intervals of the duct height at the injection location. 
Ho, rather than at the orifice diameter, D. 

The correlations developed by Holdeman and Walker (ref. 2) were 
found to be accurate for single-sided injections in Phase I of this pro- 
gram. However, these correlations were not adequate for two-sided 
injections. Improved accuracy was obtained with the new correlations 
de veloped in this program. Predictions for profiled mainstream test 
cases were obtained by superposition of mainstream temperature profiles 
and the NASA/Garrett correlations. These predictions were found to be 
good to first-order accuracy. 

The NASA/Garrett correlations also include the effects of flow area 
convergence. However, further effort is needed to improve the accuracy 
of the correlations. 
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MASS AND MOMENTUM TURBULENT TRANSPORT EXPERIMENTS 


B. V. Johnson 

United Technologies Research Center 

An experimental study of mixing downstream of nonswirling (Ref. 1) and swirling 
coaxial jets discharging into an expanded duct is being conducted to obtain data for 
the evaluation and improvement of turbulent transport models currently used in a 
variety of computational procedures throughout the propulsion community for combustor 
flow modeling. The study uses laser velocimeter (LV) and laser induced fluorescence 
(LIF) techniques to measure velocities and concentration, hot film techniques to 
qualitatively determine the time dependent characteristics of the flow and the scale 
of the turbulent structure. 

The combined LV and LIF measurement techniques are used to acquire simultaneous 
measurements of velocity and concentration of a trace material inserted into one 
fluid (Fig. 1) . A two-color LV system was used to simultaneously measure two velo- 
city components. Velocities in three directions, concentrations, momentum turbulent 
transport in three directions are obtained from a series of six measurements. The 
simultaneous data pairs are used to determine the mean r.m.s. deviation from the 
mean, the skewness and the kurtosis for the probability density functions of each 
data set (Fig. 2). These results can be compared with results from numerical codes 
or with a variety of existing or proposed turbulent transport models (Fig. 3). 

Flow visualization studies and detailed laser velocimeter/laser induced fluores- 
cence measurements have been conducted for nonswirling and swirling coaxial jets dis- 
charging into expanded ducts (Fig. 4). Typical results for the nonswirling flow con- 
dition are shown in Figs. 5 through 7. Turbulent transport results showed the exis- 
tence of counter gradient and mass trasnfer and nonisotropic turbulent transport dif- 
fusivity coefficients. 

The current effort is directed toward determining length scale and dissipation 
rates used in the turbulent transport models. Hot film measurements of the axial 
velocity component will be used to obtain a velocity spectrum and an autocorrelation 
upstream of the inlet plane, near the inlet plane and at selected locations in seve- 
ral shear layers. Dissipation rates and integral length scales will be determined 
from these measurements using conventional approximations. These results can be com- 
pared with those calculated by the numerical codes. LV and LIF measurements will 
also be made to determine the effect of changing the inlet duct geometry on the near 
inlet turublent transport characteristics. 


149 




Full Reynolds stress model 



















Page intentionally left blank 





HOST STRUCTURAL ANALYSIS PROGRAM 
OVERVIEW 
by 

R. H. Johns, Subproject Manager 


Hot section components of aircraft gas turbine engines are subjected to severe 
thermal -structural loading conditions, especially during the start-up and 
take-off portions of the engine cycle. The most severe and damaging stresses 
and strains are those induced by the steep thermal gradients induced during 
the start-up transient. These transient stresses and strains are also the 
most difficult to predict, in part because the temperature gradients and 
distributions are not well known or readily predictable, and also because the 
cyclic elastic-viscoplastic behavior of the materials at these extremes of 
temperature and strain are not well known or readily predictable. 

A broad spectrum of structures-related technology programs is underway to 
address these deficiencies. The problems are being addressed at the basic as 
well as the applied level, including participation by industry and 
universities as well as in-house at NASA Lewis. In addition to the HOST 
structural analysis program, some related program elements are being supported 
through our Base R&T program. 

One element of the structures program is developing improved time-varying 
thermal -mechanical load models for the entire engine mission cycle from 
start-up to shutdown. The thermal model refinements will be consistent with 
those required by the structural code including considerations of mesh-point 
density, strain concentrations, and thermal gradients. Models will be 
developed for the burner liner, turbine vane and turbine blade. One aspect of 
this part of the program is a thermal data transfer module recently developed 
which automates the transfer of temperatures from available heat transfer 
codes or experimental data sets to the structural analysis code. Another part 
of the program is an automated component-specific geometric modeling 
capability which will produce 3-D finite element models of the components. 
Self-adaptive solution strategies will be developed and included to facilitate 
selection of appropriate elements, mesh sizes, etc. 

Another major part of the program is the development of new and improved 
nonlinear 3-D finite elements and associated structural analysis programs, 
including the development of temporal elements with time-dependent properties 
to account for creep effects in the materials and components. Two contracts 
were recently signed to accomplish these developments. Improved constitutive 
modeling methods to facilitate improved prediction of cyclic thermomechanical 
viscoplastic material behavior are also under development. Experimental 
facilities to aid in developing and verifying theories and models are 
currently being established in-house at Lewis. 

Further explanation and details about the various structures program elements 
mentioned above are given in the following write-ups. 
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STRUCTURAL ANALYSIS JTS ROLE IN HOST 

GOAL: 

DEVELOP AND VALIDATE INTEGRATED TIME-VARYING THERMOMECHANICAL LOAD MODELS, 
COMPONENT-SPECIFIC AUTOMATED GEOMETRIC MODELING AND SOLUTION STRATEGY 
CAPABILITIES, AND ADVANCED INELASTIC ANALYSIS METHODS INCLUDING PLASTICITY AND 
CREEP EFFECTS FOR NONLINEAR ANISOTROPIC FINITE ELEMENT STRUCTURAL ANALYSIS 
AND DESIGN COMPUTER CODES FOR TURBINE ENGINE HOT SECTION COMPONENTS 

PROGRAM INTEGRATION PROGRAM ELEMENTS: 

• THERMAL/STRUCTURAL DATA 
TRANSFER MODULE 


• THERMAL/MECHANICAL LOAD/MISSION 
AND COMPONENT-SPECIFIC STRUCTURAL 
MODELS 


• 3-D INELASTIC ANALYSIS METHODS 

• BURNER LINER CYCLIC RIG (STRUCTURAL 
ASPECTS) 

• HIGH TEMPERATURE STRUCTURES 

LABORATORY cs-83-o<543 



HOST 

STRUCTURAL ANALYSIS PROGRAM 


OBJECTIVE: 

TO DEVELOP AND VALIDATE INTEGRATED TIME- VARYING THERMAL/MECHANICAL 
LOAD MODELS, COMPONENT-SPECIFIC AUTOMATED GEOMETRIC MODELING AND 
SOLUTION STRATEGY CAPABILITIES, AND ADVANCED INELASTIC ANALYSIS 
METHODS AND CONSTITUTIVE MODELS, INCLUDING PLASTICITY AND CREEP 
EFFECTS, FOR NONLINEAR, ANISOTROPIC, FINITE ELEMENT STRUCTURAL 
ANALYSIS AND DESIGN COMPUTER CODES, 
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BURNER LINER THERMAL/STRUCTURAL LOAD MODELLING 
R. J. Maffeo 

General Electric Company 
Aircraft Engine Business Group 

The objective of this program is to develop a thermal data transfer computer 
program module for the Burner Liner Thermal -Structural Load Modelling Program. 

This will be accomplished by (1) reviewing existing methodologies for thermal data 
transfer and selecting three heat transfer codes for application in this program, 

(2) evaluating the selected codes to establish criteria for developing a computer 
program module to transfer thermal data from the heat transfer codes to selected 
stress analysis codes, (3) developing the automated thermal load transfer module, 
and (4) verifying and documenting the module. 

In aircraft turbine engine hot section components, cyclic thermal stresses 
are the most important damage mechanism. Consequently, accurate and reliable 
prediction of thermal loads is essential to improving durability. To achieve this 
goal, a considerable effort over the past 20 years has been devoted to the 
acquisition of engine temperature test data, as well as the development of accurate, 
reliable, and efficient computer codes for the prediction of steady-state and 
transient temperatures and for the calculation of elastic and inelastic cyclic 
stresses and strains in hot section components. There is a need for continued 
development of these codes, because the availability of more accurate analysis 
techniques for complex configurations has enabled engine designers to use more 
sophisticated designs to achieve higher cycle efficiency and reduce weight. 

It has become apparent in recent years that there is a serious problem of 
interfacing the output temperatures and temperature gradients from either the heat 
transfer codes or engine tests with the input to the stress analysis codes. 

With the growth in computer capacity and speed and the development of input 
preprocessors and output postprocessors, the analysis of components using hundreds 
and even thousands of nodes in the heat transfer and stress models has become 
economical and routine. This has exacerbated the problem of manual transfer of 
output temperatures from heat transfer nodes to stress analysis input to where 
the engineering effort required is comparable to that required for the remainder 
of the analysis. Furthermore, a considerable amount of approximation has been 
introduced in an effort to accelerate the process. This tends to introduce errors 
into the temperature data which negates the improved accuracy in the temperature 
distribution achieved through use of a finer mesh. There is, then, a strong need 
for an automatic thermal interface module. 

The overall objectives of this thermal transfer module are that it handle 
independent mesh configurations, finite difference and finite element heat 
transfer codes, perform the transfer in an accurate and efficient fashion and 
that the total system be flexible for future applications. 

Based on our study of existing thermal transfer modules, and our experience 
with our two-dimensional in-house transfer code, three levels of program development 
criteria were identified. 


Work done under NASA contract NAS3-23272. 
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Level I contains the general criteria which must be satisfied for a usable 
product. These include: 

-Independent Heat Transfer and Stress Model Meshes 
-Accurate Transfer of Thermal Data 
-Computationally Efficient Transfer 
-User Friendly Program 
-Flexible System 

Level II contains specific criteria which must be satisfied to meet require- 
ments associated with gas turbine applications, such as: 

-Internal Coordinate Transformations 
-Automatic Exterior Surfacing Techniques 
-Geometrical and Temporal Windowing Capability 

Level III has criteria which are desirable but not necessary. Total automation 
of these could be accomplished in future enhancements of the transfer code. Items 
that fall in this area include: 

-Automatic Scaling of Temperatures Based on Engine Power Setting 
-Altered Stress Geometry 

-Automatic Handling of Temperature Discontinuities 

All of the Level I and Level II criteria have been developed and implemented 
into the thermal load transfer code. This code is being used at General Electric 
and has been used in conjunction with a three-dimensional model of a combustor 
liner for the verification phase of this contract. The existing transfer module 
can process heat transfer results directly from the MARC and SINDA programs and 
will output temperature information in the forms required for MARC and NASTRAN. 

The input and output routines in the module are very flexible and could easily be 
modified to except data from other heat transfer codes and format data to other 
stress analysis codes. 

This thermal load transfer module has been shown to efficiently and accurately 
transfer thermal data from dissimilar heat transfer meshes to stress meshes. 

The fundamental part of the code, the 3D search, interpolation and surfacing 
routines, have much more potential. They form an outstanding foundation for 
automatic construction of embedded meshes, local element mesh refinement, and the 
transfer of other mechanical type loading. 
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Problem 


• Large Heal Transler/Stress Models 

• Different Mesh Densities 

• Finite Difference vs Finite Element Codes 

• Thermal Transfer Time Consuming and Error Prone 


Objectives 

• Transfer Temperatures From a Heat Transfer Study 
to a Stress Analysis 

- Independent Meshes 

- Accurate/Efficient Transfer 

- Flexible 


Objective 
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Thermal T ransler Module 



Overall System 


Enhanced Program Features 


Exterior Point Method Ixoparametric Mapping 



162 







Verification Study 


® 413 Heel Trsnstef Element* 
9 837 Corner Hodes 



Temperature Verification Plots 
Convex Side Mid Section 
Large Model (Airfoil) 


2100 


2000 

“ Cl Heat Transfer Result 


O Hut Transfer fteeutt 

isooC 

t 

V A Thermal Transfer Module 

1000 


Temperature 

cr n 1 

1700 

>wDt rfi 

1000 


1300 


1400 

1300 

1 1 ! 1 i i i i i i i 


0 0.1S 0.30 0.45 0.80 .0.7$ 0.90 1JH 1.20 143 I SO IjW 


Distance 


Future Applications 

Automatic Boundary Condition Assignment 



1J0 


163 



Page intentionally left blank 



COMPONENT-SPECIFIC MODELING* 

M. L. Roberts 
General Electric Company 


Introduction 


Modern jet engine design imposes extremely high loadings and temperatures on 
hot section components. Fuel costs dictate that minimum weight components be used 
wherever possible. In order to satisfy these two criteria, designers are turning 
toward improved materials and innovative designs. Along with these approaches, 
however, they must also have more accurate, more economical, and more comprehensive 
analytical methods. 

Numerous analytical methods are available which can, in principle, handle any 
problem which might arise. However, the time and expense required to produce 
acceptable solutions is often excessive. This program addresses this problem by 
setting out a plan to create specialized software packages, which will provide the 
necessary answers in an efficient, user-oriented, streamlined fashion. Separate 
component-specific models will be created for burner liners, turbine blades, and 
turbine vanes using fundamental data from many technical areas. The methods 
developed will be simple to execute, but they will not be simple in concept. The 
problem is extremely complex and only by a thorough understanding of the details 
can the important technical approaches be extracted. The packaging of these 
interdisciplinary approaches into a total system must then conform to the modular 
requirements for useful computer programs. 

Objective 

The overall objective of this program is to develop and verify a series of 
interdisciplinary modeling and analysis techniques which have been specialized to 
address three specific hot section components. These techniques will incorporate 
data as well as theoretical methods from many diverse areas, including cycle and 
performance analysis, heat transfer analysis, linear and nonlinear stress analysis, 
and mission analysis. Building on the proven techniques already available in these 
fields, the new methods developed through this contract will be integrated to 
provide an accurate, efficient, and unified approach to analyzing combustor burner 
liners, hollow air-cooled turbine blades, and air-cooled turbine vanes. For these 
components, the methods developed will predict temperature, deformation, stress, 
and strain histories throughout a complete flight mission. 

Background 

This program, to a great extent, will draw on prior experience. This base of 
experience is invaluable for understanding the highly complex interactions among 
all the different technical disciplines as well as for estimating the importance of 
different engine parameters. In particular, there are four specific areas in which 
experience will be especially beneficial. 


*This work was done under NASA Contract NAS3-23687. 
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First, with the recent increases in fuel costs, greater emphasis has been 
placed on more accurate solutions for stresses and strains in order to understand 
and improve the durability and life of hot section components. Conventional linear 
elastic analyses are no longer sufficient; instead, they now provide the boundary 
values for more refined creep and plasticity calculations. These nonlinear 
analyses are now performed routinely as part of the design process at General 
Electric. This extensive experience with these plasticity and creep methods will 
contribute directly to developing component specific models. 

Second, advances in 3-D modeling capability are being achieved by the concepts 
developed under the NASA-supported ESMOSS program. ESMOSS concepts will provide 
the basis to develop an efficient modeling system for geometric and discretized 
models of engine components. 

Third, the NASA-funded Burner Liner Thermal/Structural Load Modeling Program 
will contribute strong support to this program. The specific area addressed, 
transfer of data from a 3-D heat transfer analysis model to a 3-D stress analysis 
model, will provide the background and framework for the data interpolation 
required for all thermomechanical models in this contract. 

Fourth, over the past 10 years. General Electric has developed internally a 
family of computer programs: LASTS, OPSEV, and HOTSAM. These programs all have the 
common thread of using selected points from cycle data, heat transfer, and stress 
analyses, and a decomposition/synthesis approach to produce accurate values of 
temperature, stress, and strain throughout a mission. These programs are totally 
consistent with the overall objectives of this program, and represent a proven 
technology base upon which the component specific models will be developed. 
Significant advances to be made are the inclusion of nonlinear effects and the 
introduction of improved modeling and data transfer techniques. 

Approach 

The program is organized into nine tasks which can logically be separated into 
two broadly parallel activities (Figure 1). On the right of Figure 1 we have the 
Component Specific Thermo-Mechanical Load Mission Modeling path. Along this path a 
Decomposition/Synthesis approach will be taken. In broad terms this means 
developing methods to generate approximate numerical models for the engine cycle 
and the aerodynamic and heat transfer analyses needed to provide the input 
conditions for hot parts stress and life analysis. 

The left path. Component Specific Structural Modeling will provide the tools to 
develop and analyse finite element nonlinear stress analysis models of combustor 
liners and turbine blades and vanes. These two paths are shown in more detail in 
Figures 2 and 3. 

Software Development, Task IV consists of planning and writing the computer 
programs for both paths, with the necessary interconnections, using a structured, 
top down approach. 

In the Thermomechanical Load Mission Modeling portion of the program (Figure 2) 
we will develop, in Task III, a Thermodynamic Engine Model which will generate the 
engine internal flow variables for any point on the operating mission. The method 
for doing this is described below. Task V will develop techniques to decompose 
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flight missions into characteristic mission segments. In Task VII a Thermo- 
mechanical Mission Model will be developed. This will use the flow variables from 
the Thermodynamic Model to determine metal temperature and pressure distributions 
for a representative combustor liner and turbine blade and vane. 

Individual tasks for the Structural Modeling activity are shown in Figure 3. 

The requirements of Software Design, Task II, will be factored into Task VI, the 
evaluation of the structural analysis methods which were selected for evaluation in 
Task I. Task VIII will provide the capability for structurally modeling current 
state-of-the-art combustor liners and hollow turbine blades and vanes, given the 
defining dimensional parameters. These parameters will be chosen to facilitate 
parametric studies. 

The component specific models will be developed in two steps. In the first a 
geometric model will be defined. In the application of the Component Specific 
Modeling Program this data will then be transferred to the Thermomechanical Load 
Mission Model to provide the geometry for determining component pressures and 
temperatures. Thus, a data transfer link will be developed to do this in Task IV, 
Software Development. The capability for generating from the geometric model a 
discretized, finite element model will also be part of Task VIII. At this point 
another link between the two paths will be needed to transfer the component 
temperatures and pressures from the Thermomechanical Load model to the finite 
element model, interpolating the data as needed to define nodal temperatures and 
pressures. This also will be completed in Task IV. 

The final function in Task VIII will be the development of component-specific 
stress analysis models for the three components to perform cyclic elastic, plastic 
and creep analyses using loading conditions defined by the Thermomechanical Load 
Models Progress. 

At this time considerable progress has been made on the Thermodynamic Engine 
Model. The model is being developed as a simple calculational tool which will take 
as inputs the three variables, altitude (h), Mach number (M) and power level (PL) 
for the allowed flight map of an engine, as shown in Figure 4. In addition, 
ambient temperature deviations from the standard atmosphere, airframe bleed air 
requirements and engine deterioration can also be included as part of the input to 
the Thermodynamic Model. For each input condition, specified by h, M and PL the 
Thermodynamic Engine Model will calculate gas weight flow (w), temperature (t) and 
pressure (p) at selected aerodynamic engine stations, as needed to determine 
component thermal loadings. These stations are shown in Figure 5. 

The technique for developing a Thermodynamic Engine Model is shown in Figures 6 
and 7. The engine to be analyzed must be defined thermodynamically by an engine 
cycle deck (computer program) which can be run to generate the internal flow 
variables at the chosen aerodynamic stations (Figure 5). To encompass the complete 
engine operating map (Figure 4), 148 operating points are chosen and w, t and p are 
calculated using the cycle deck for the selected stations, as well as Nj and N 2 , 
the fan and core speeds. From this station data an Engine Performance Cycle Map is 
constructed. This is essentially a set of three-dimensional data arrays which map 
the station data (w, t, p, Nj and N 2 ) on to the engine operating map (Figure 4). 
Given an arbitrary operating point defined by h, M and PL it is then, in principle, 
possible to interpolate on the Engine Performance Cycle Map to determine station 
data. In practice the station parameters are nonlinear functions of the input 
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parameters and considerable effort was needed to develop these multi-dimensional 
interpolations. The computer program used to generate the Engine Performance Cycle 
Map from the engine cycle desk output has been developed as part of Task III. The 
functioning of the Thermodynamic Engine Model is shown in Figure 7. Given an 
engine mission, as shown schematically in Figure 8 it can be defined by values of 
the input variables h, M and PL at selected times through the mission. Using these 
input variables and the Engine Performance Cycle Map an Interpolation Program, now 
being developed in Task III of this program, will calculate engine station 
parameters throughout the mission (Figure 7). These are then used to define 
Station Mission Profiles of w, t, p, Nj, and N 2 , as functions of time at each 
aerodynamic station. These Station Mission Profiles are become the input to the 
Thermomechanical Engine Model. 

The Thermomechanical Model is less well developed at this time than the 
Thermodynamic Model. Its form will be based on types of correlations previously 
developed within General Electric. Figure 9 shows a representative correlation for 
a turbine vane. Metal temperatures at various points on the vane Tva are 
correlated in terms of a vane overall cooling effectiveness, n, and station gas 
temperatures T 3 at compressor discharge, and T 4 at combustor discharge. Using the 
Station Mission Profiles it will be possible to calculate the temperatures at 
selected locations on each component as functions of time, given the input 
parameters h, M and PL that define the engine mission. These then will provide the 
boundary conditions for the component stress analyses. 

On the Component Specific Structural Modeling path, concepts have been defined 
and are being implemented. Additional evaluations are needed, however, before they 
can be presented for discussion. 

Conclusion 


When completed this program will provide a stress analysis system for hot 
section parts that will allow the component designers to evaluate quickly the 
effects of mission variations, be easy to use, cost effective, and make a 
significant contribution to assessing hot section durability. 
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FIGURE 1 . COMPONENT SPECIFIC MODELLING - BASE PROGRAM 



FIGURE 2. COMPONENT SPECIFIC 

THERMOMECHANICAL LOAD MISSION MODELLING 
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FIGURE 3. COMPONENT SPECIFIC STRUCTURAL MODELING 



FIGURE 4. ENGINE OPERATING MAP 
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FIGURE 7, THERMODYNAMIC ENGINE MODEL 



Tise - Minutes 


FIGURE 8, TYPICAL FLIGHT CYCLE 


172 











T* Design 



Van* Overall Cooling EHactlvanou, T) =( T «-Ty A , n , s , x VtT^-T 3 ) 


FIGURE 9. TURBINE VANE COOLING EFFECTIVENESS 


173 


Page intentionally left blank 



3D INELASTIC ANALYSIS METHODS FOR HOT SECTION COMPONENTS 

L.T. Dame 
R.L. McKnight 
General Electric Company 
Aircraft Engine Business Group 

The objective of this research is to develop an analytical tool capable of 
economically evaluating the cyclic time-dependent plasticity which occurs in hot 
section engine components in areas of strain concentration resulting from the com- 
bination of both mechanical and thermal stresses. The techniques developed must 
be capable of accommodating large excursions in temperatures with the associated 
variations in material properties including plasticity and creep. 

The overall objective of this proposed program is to develop advanced 3-D 
inelastic structural/stress analysis methods and solution strategies for more 
accurate and yet more cost-effective analysis of combustors, turbine blades, and 
vanes. The approach will be to develop four different theories, one linear and 
three higher order with increasing complexities including embedded singularities. 

The objective will be achieved through a four-phase program consistent with 
the NASA Statement of Work. 

In Task I, the linear formulation theory will be developed. These will con- 
sist of three linear formulation models in which stress, strain, and temperature 
are linear functions of the spatial coordinates; and the increments in loading, 
temperature, and time are linear. Three constitutive relations will be developed 
for these linear formulation models each capable of predicting elastic, plastic, 
thermal, and creep strains and cyclic effects. One constitutive relation will 
be approximate, one will be of the current genre, and one will be a unified 
theory . 


In Task II, the polynomial formulation theory will be developed. These 
will consist of three polynomial formulation models in which stress, strain, and 
temperature are polynomial functions of the spatial coordinates, and the incre- 
ments in loading temperature and time are quadratic. They will also accommodate 
two-intersecting embedded discontinuities. Three constitutive relations will 
be associated with these polynomial formulation models. 

In Task IV, the special functions theory will be developed. These will 
consist of three special function formulation models in which stress, strain, 
and temperature are special functions of the spatial coordinates and the incre- 
ments in loading, temperature, and time are special functions. These models 
will accommodate eight intersecting embedded similar discontinuities and have 
three associated constitutive relations . 

In Task V, the general functions theory will be developed. These will 
consist of three general function formulation models in which stress, strain, and 
temperature are general functions of the spatial coordinates and the increments 
in loading, temperature, and time are general functions. 

These models will accommodate eight intersecting embedded different dis- 
continuities and have three constitutive relations associated with them. One of 
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the constitutive relations will be more complex than those used for the special 
functions theory. 


Task III and VI are reporting requirements . 

Work began on this program in April 1983. A literature survey was conducted 
in order to identify various solution techniques and constitutive models . In 
developing an economical computer tool it is essential to consider the inter- 
actions of solution technique, constitutive model, integration scheme, load 
incrementing, temperature incrementing, and time incrementing procedures. 

For time dependent plasticity (or creep) problems the solution technique 
most commonly used in finite element codes is a constant stiffness scheme. Any 
economical scheme for thermo-mechanical cycling will of necessity require automatic 
time incrementing. Initially, the constitutive equation subroutines are being 
developed in the context of a finite element code with a constant stiffness itera- 
tion procedure, automatic time incrementing, linear variation of loads and boundary 
conditions and isothermal conditions . The introduction of time varying tempera- 
tures will require a modified solution procedure. The isothermal Bodner formula- 
tion has been implemented in a single element computer program and compared with 
test data and calculations in Reference 1. 

In view of the increasing availability of vector processors, work is also 
underway to develop more efficient numerical procedures for solving sets of 
linear algebraic equations on such machines. Iterative procedures such as the 
Jacobi or Gauss-Seidel methods are more feasible with vector machines (Reference 
2-3) , especially if the stiffness matrix is stored in column form. Improvements 
in forming stiffness matrices may also be realized (Reference 4) which would 
greatly enhance the calculation time in tangent stiffness schemes . Eigenvalue 
calculation procedures may also be more economical using vector processors 
(Reference 5) . 
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3-D INELASTIC ANALYSIS METHODS 
FOR HOT SECTION COMPONENTS (BASE PROGRAM) 


E. So Todd 

Pratt & Whitney Engineering 


The objective of this program is to produce a series of new computer codes that 
permit more accurate and efficient three-dimensional inelastic analysis of selected 
hot section components - combustor liners, turbine blades and turbine vanes. The 

computer codes embody a progression of mathematical models and are streamlined to 
take advantage of geometrical features, loading conditions, and forms of material 
response that distinguish each group of selected components. 

Software in the form of stand-alone codes and modules for use in General Purpose 
Structural Analysis (GPSA) programs will be provided by Pratt & Whitney Aircraft 
(P&WA) with assistance from three uniquely well qualified subcontractors: MARC 
Analysis Research Corporation (MARC), State University of New York at Buffalo 
(SUNY-B), and United Technologies Research Center (UTRC). Primary development of 
special finite element models will be accomplished by MARC, while mechanics of 
materials models and constitutive formulations will be assembled by UTRC. Development 
of advanced formulation (boundary element) models will be shared by P&WA and SUNY-B. 
Verification of the various analysis packages will be done by P&WA. 

The technical effort of this program will be conducted over a period of 
twenty-four (24) months and will involve three (3) distinct tasks: 

Task I ~ Linear Theory 


Under this task a linear theory will be developed consisting of three linear 
formulation models: 1) mechanics of materials, 2) special finite elements, and 3) 
advanced formulation (boundary element) models. The linear theory is defined to mean 
that: 1) the stresses or strains and temperatures in the generic modeling region are 
linear functions of the spatial coordinates, and 2) the solution increments for load, 
temperature and/or time are extrapolated linearly from previous information. The 
emphasis in the linear theory for three-dimensional inelastic analysis will be on a 
methodology using a large number of generic modeling regions and a large number of 
increments, but simple calculations per increment. 

Task II - Polynomial Theory 


This task will involve the development of a polynomial theory consisting of three 
polynomial formulation models: 1) mechanics of materials, 2) special finite elements, 
and 3) advanced formulation (boundary element) models. The polynomial theory is 
defined to mean that: 1) the stresses or strains and temperatures in the generic 
modeling region are at least quadratic functions of the spatial coordinates, and 2) 
the solution increments for load, temperature and/or time are extrapolated at least 
quadratically from previous information. The emphasis in the polynomial theory for 
three-dimensional inelastic analysis will be on a methodology using a smaller number 
of generic modeling regions with two intersecting, embedded discontinuities and a 
smaller number of increments, but more complex calculations per increment than for 
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the linear theory (Task I). 


Task III - Reporting Requirements 

This task includes periodic technical, financial and schedular reporting and, at 
the conclusion of the technical effort, the submittal of a final Contractor Report. 

Optional Program 

An optional program, to be exercised at the discretion of the Government, will 

extend the base computer codes to include higher order representations of strain in 
space and time and to deal effectively with more complex collections of 

discontinuities such as cooling holes and coating cracks. 
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VALIDATION OF STRUCTURAL ANALYSIS METHODS USING 
BURNER LINER CYCLIC RIG TEST DATA 


R. Thompson 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


INTRODUCTION 

The objectives of the HOST Burner Liner Cyclic Ri q Test Program are 
basically threefold: fl) to assist in developing predictive tools needed to 

improve design analyses and procedures for the efficient and accurate predic- 
tion of burner liner structural response; (?) to calibrate, evaluate and 
validate these predictive tools by comparing the predicted results with the 
experimental data generated in the tests; and (3) to evaluate existing as well 
as advanced temperature and strain measurement instrumentation, both contact 
and non-contact, in a simulated engine cycle environment. The data generated 
will include measurements of the thermal environment (metal surface tempera- 
tures) as well as structural (strain) and life (fatigue) responses of simu- 
lated burner liners and specimens under controlled boundary and operating 
conditions. These data will be used to calibrate, compare and validate 
analytical theories, methodologies and design procedures, as well as improve- 
ments in them, for predicting liner temperatures, stress-strain responses and 
cycles to failure. Comparison of predicted results with experimental data 
will be used to show where the predictive theories, etc. need improvements. 

In addition, as the predictive tools, as well as the tests, test methods, and 
data acquisition and reduction techniques, are developed and validated, a 
proven, integrated analysis/experiment method will be developed to determine 
the cyclic life of a simulated burner liner. 


TEST RIGS 

Figure 1 includes a list of the test rigs under construction, with 
anticipated delivery dates and the basic specimens and burner liner components 
to be analyzed and tested in each rig. Flat plate specimens, with and without 
holes, will be tested in the Box Rig while tubular specimens, with and without 
holes, and subelements of burner liners will be tested in the Annular Rig. 

Each test will be increasingly more complex than the preceding one, both from 
the standpoint of geometry as well as the data acquisition and reduction 
requirements. Correspondingly, the structural analysis will become more com- 
plex in the progression of test configurations. For example, figure 2 shows 
the progression of flat plate specimens to be tested and structurally analyzed. 

A schematic of the Box Rig, being fabricated in-house, is shown in fig- 
ure 3. Four quartz lamps are used to heat the 8 x 5 x 0.05 inch flat plate 
specimens. Both water and air cooling are used. A viewing port provides for 
visual inspection of the specimen and data acquisition. Provisions for 
instrumentation include 30 thermocouples and 10 strain gage connections. 
Initially, the Box Rig will provide temperature and limited strain measure- 
ments using commercially available thermocouples and wire resistance strain 
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gages. But as advanced temperature and strain measurement instrumentation is 
developed and becomes available under the HOST program and elsewhere, the Box 
Rig will serve to evaluate the advanced instrumentation, for example, infrared 
thermal imaging camera for temperature mapping, thinfilm thermocouples and 
strain gages, laser speckle strain measurements, etc. The instrumentation 
with the greatest potential will then be installed and used on the Annular Rig. 

A schematic of the Annular Rig is shown in figure 4. The test section, 
being built by P&W as part of a cooperative effort between NASA Lewis and 
P&W, has 112 quartz lamps to heat the 21-inch annular test specimen. The 
length of the heated portion of the test specimen is 10 inches. Both water 
and air cooling are used. Three active viewing ports are provided for visual 
inspection of the specimen and data acquisition. Seven quick disconnect 
instrumentation panels are shown in the figure. These panels provide for 100 
thermocouple, 28 strain gage and 21 pressure transducer connections. Both 
the Box Rig and Annular Rig Test Programs will be conducted at NASA Lewis in 
ECRL-1 . 

Test conditions and variables to be considered in each of the test rigs 
and test configurations, and also used in the structural analysis, for valida- 
tion of the predictive theories and tools will include thermal and mechanical 
load histories (simulating an engine mission cycle), different boundary condi- 
tions (fixed, free, etc.) a variety of specimen and subelement geometries 
(including advanced burner liner structural design concepts), different 
materials (initially Hastelloy-X will be used), various cooling schemes and 
cooling hole configurations, and the simulation of hot streaks. Based on 
these test conditions and test variables, test matrices for each rig and con- 
figuration will be developed with the intent to verify the predictive tools 
over as wide a range as is feasible, using the simplest of possible tests. 

Table I is a sample test matrix for a flat plate without holes. 


EXPERIMENTS 

To verify the feasibility of the tests and test rig designs, and identify 
potential test problems and analysis/experiment complications, a limited 
number of preliminary experiments were conducted and structural analyses per- 
formed using 5 x 8 x 0.05 inch Hastelloy-X flat plate specimens. The speci- 
mens were tested in the box rig shown in figure 5. Three quartz lamps, placed 
about 2 inches apart, center-to-center, about 1-1/2 inches from the reflector 
and about 1-1/2 inches below the test specimen, were used to heat the speci- 
men. The lamps were air cooled through a manifold located at the bottom of 
the box. Air pressure for cooling the lamps ranged from 4 psi at the lowest 
power setting (30 W) to 60 psi at the maximum power (18 000 W) setting. The 
power was controllable in steps from 30 to 18 000 watts. No water cooling was 
provided in this rig. The top of the plate could be viewed through a cut-out 
in the top cover. Some of the more salient results of these tests are pre- 
sented below. 

The test specimen was held fixed between two identical frames by tighten- 
ing 10 bolts, as shown in figure 6. The frames and specimen were positioned 
horizontally above the lamps and were held in place by leg supports (long 
bolts) which rested on the bottom of the box. 
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Both the frame and test specimen were instrumented. A total of '22 
thermocouples (0. 032-in. -diam. Chromel-Al umel ) were mounted on the frame. 

Ten thermocouples were mounted on the specimens; three on the outer face (cool 
side) and seven on the inner face (hot side) exposed to the quartz lamps. 

The thermocouples were mounted along the two center axes of the plate. High- 
temperature resistance strain gages, BLH Model HT-12/2-5A, were mounted on two 
specimens. The strain gage is made of platinum-tungsten alloy and works for 
temperatures up to 1200 F. The locations of the thermocouples and strain 
gages on the frame and test specimen are shown in figure 7. 

Three flat plate specimens without holes were tested. The first series 
of tests were conducted to verify that the rig was working properly, to verify 
the rigs performance capability and durability, and to observe the failure 
mode of the plate (buckling). The frame was instrumented with thermocouples 
but the plate was not. Temperature measurements of the frame indicated that 
depending on the rate at which the power was increased, the temperature 
difference between the top and bottom half of the frame could vary from 0 to 
100° F at temperatures below 500° F (maximum frame temperature) to a maximum 
variation of 0 to 400° F at 1000° F and above. By increasing the power 
slowly, this difference could be controlled and kept to a minimum. The con- 
clusion reached is that control of the frame temperature is needed to prevent 
buckling of the plate for a fixed edge boundary condition. 

In the second series of tests, both temperature and strain measurements 
were obtained. Temperature variations on the plate and frame were recorded 
for several power settings. For a single power setting, representative plots 
of the plate and frame temperatures are shown in figures 8 and 9. The power 
level was stepped up from 375 volts (12 100 W) to 435 volts (15 500 W) and 
held constant for about 7 minutes. The delta increase in frame temperature 
ranged from 50° to 100° o F whereas the delta increase in plate temperature 
ranged from 100° to 230° F. The initial maximum/minimum temperatures for 
the frame and plate were 510° (TC 12)/260° F (TC 4) and 1270° (TC 4)/830 F 
(TC 1), respectively. The final maximum/minimum temperatures for the frame 
and plate were 600°/310° F and 1500°/1000° F respectively. The frame tempera- 
ture increased almost linearly with time while the plate temperature reached 
almost steady state conditions at about 2 minutes into the step increase in 
power. Contrary to our expectations, a non-uniform plate temperature was 
observed for this and all tests in the longitudial direction, while a more 
uniform temperature was observed across the plate (see figs. 8 and 9) for o all 
tests. Through the thickness temperatures varied widely, ranging from 20° to 
200° F. Unexpl ainably, the data show that the front half o of the plate was 
much hotter than the back half of the plate (by about 530° F for data shown in 
fig. 9), and the front half of the frame (both bottom and top) was hotter than 
the back half (by about 100° F for the data shown in fig. 8) for all tests. 

The most obvious reasons for this anomaly were ruled out, such as a slopping 
plate, distance between lamps and test specimen, cooling air flow patterns, 
plugged air holes, and misalignment of the cooling air manifold. These tests 
showed: it was possible to control the plate temperature, both transient and 

steady state, by varying power settings and air flows; it was possible to 
achieve desired maximum/minimum temperatures in the plate, approximately 
1700° F temperature at a power setting of 440 volts and a 900 F temperature 
at a power setting of 300 volts; it was possible to control the frame tempera- 
ture to a degree; it was possible to control strains by varying both frame and 
plate temperatures for fixed edge boundary conditions; it was not possible to 
achieve uniform surface temperatures in the plate either in the horizontal or 
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longitudinal directions; it was not possible to achieve the desired or uniform 
through the thickness temperature gradients; and strain measurements were not 
obtained above 1200° F. 

The primary objective of the third series of tests was to demonstrate 
whether or not a simulated engine mission cycle could be obtained. The test 
results showed that controlled complex, cyclic thermal cycles could be imposed 
on a flat plate specimen. Other than this, the test results obtained were 
similar to those obtained from the second series of tests. 


ANALYSIS 

A structural analysis of a flat plate specimen without holes was per- 
formed using MARC,, a general purpose, nonlinear finite element structural 
analysis computer program. Walker's viscoplastic constitutive model was 
incorporated into MARC to account for the interaction between creep and plas- 
tic deformations. 

The plate was analyzed using a four node plane stress element. The mesh 
configuration was set up so that direct comparisons could be made between 
models with a single hole and without a hole. The mesh conf igurations for 
these two models are shown in figure 10. The configurations are doubly symme- 
tric, thus only a quarter of the plate is modelled. 

A parametric study was conducted to examine different boundary conditions 
for a flat plate without holes. The results indicated that both the test 
specimen and frame would have to be analyzed together and that controlled, 
variable boundary conditions were necessary to reproduce the desired stress 
levels in the plate. Controlling the boundary conditions was done by pro- 
viding truss members at the edges of the plate as shown in figure 11. By 
appropriately adjusting the temperature difference between the plate and the 
edge members, the stress levels can be varied in the plate. Thus, boundary 
conditions ranging from free edge (with the plate and edge bars at the same 
temperature) to fixed edge (with the edge bars kept at a reference temperature 
while the plate temperature is increased), as well as all intermediate condi- 
tions is possible. Further, depending on whether the edge bar temperatures 
lag behind the plate temperatures or vice versa, compression or tension may be 
produced in the plate. This reversal of stress would make it possible to 
trace hysteresis loops and thus study thermal ratcheting, compare experiment 
with predictions, etc. 

A simulated engine mission cycle was studied assuming uniform tempera- 
tures in the plate and the edge bars, with a 5-percent temperature difference 
between the plate and frame as shown in figure 12. The initial temperature 
was set at 940° F and the cycle consisted of a 25-second ramp up to a maximum 
temperature of 1740° F, a 40 second hold time at this temperature level and a 
30-second ramp down to a temperature of 940° F. Figure 12 shows the time 
variation of stress in element 18 (see fig. 10). The stress-strain variation 
for this element is plotted in figure 13. No direct comparisons could be made 
between this analysis and the experiments described earlier because of the 
nonuniformity of plate and frame temperatures obtained from the experiment. 
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CONCLUSIONS AND FUTURE RESEARCH 


The preliminary experiments demonstrated this feasibility of the rig to 
achieve the desired plate temperatures, temperature gradients and cyclic 
thermal load histories, although uniform plate temperatures were not achieved. 
Finite element models of plates, with and without holes, and frames were gen- 
erated and structural analyses were performed. However, a direct comparison 
between experiment and analysis could not be done because of the nonuniform 
plate temperatures obtained in the experiments. 

The new Quartz Lamp Box Rig under construction has four quartz lamps, 
with a maximum power of 24 000 watts. This rig is expected to produce much 
more uniform plate temperatures than those obtained in the rig just described. 
In addition, computer controlled operation and data acquisition will allow for 
more accurate and extensive tests. The ESCORT II Data Acquisition System, 
available at Lewis, will be used for data acquisition and control of the 
experiment. The data obtained, both temperature and strain measurements, will 
be stored on tape for processing, reduction and analysis at a later date. 

Also, an infrared thermal video imaging system will be integrated into the 
system to obtain temperature maps of the plate. 

The Annular Rig, as described earlier, under construction, will provide 
data on cylindrical shells and subelements of burner liners. The same fea- 
tures available for the Box Rig will be carried over to the Annular Rig. 
Preliminary structural analyses of plates and shells will be used to guide 
the direction of the experimentation in both the Box and Annular Rigs. Other 
analyses will be conducted to predict the material stress-strain response 
using the measured temperature distributions obtained from the tests. These 
analyses will provide a basis for comparing analytical predictions, for exam- 
ple, using several constitutive models, with the experimental data for valida- 
tion and subsequent selection of improved analysis methods to predict, more 
efficiently and accurately, the structural response of burner liners. 
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TABLE I - SAMPLE TEST MATRIX FOR FLAT PLATE 
TEST CONFIGURATION - FLAT PLATE WITHOUT HOLES 
MATERIAL - HASTELLOY X 


TEST 

NO. 

TEMPERATURE 

STATE 

TEMPERATURE HISTORY 

RAMP UP HOLD TIME 

TIME, sec. TEMP .,°F TIME, sec.TEMP.,°F 

RAMP DOWN 
TIME, sec.TEMP., °F 

BOUNDARY COND. 
FOR 

PLATE EDGES 

NUMBER OF 
CYCLES 

(TIME BET. CYCLES 

1 

STEADY STATE- 
UNIFORM TEMP. 

20 

RT-600 

300 

600 

— 

FIXED EDGES 

— 

2 

II 

20 

RT-800 

300 

800 

- 

ii 

-- 

3 

II 

20 

RT-1000 

300 

1000 

- 

ii 

-- 

4 

II 

20 

RT-1200 

300 

1200 

- 

ii 

.. 

5 

II 

20 

RT-1400 

300 

1400 

- 

u 

- 

6 

II 

20 

RT-1600 

300 

1600 

- 

II 

- 

7 

II 

20 

RT-1800 

300 

1800 

- 

II 

- 

8 

II 

40 

RT-1000 

300 

1000 

- 

It 

— 

9 

II 

40 

RT-1400 

300 

1400 

- 

II 

— 

10 

II 

40 

RT-1800 

300 

1800 

- 

II 

- 

11 

II 

60 

RT-1000 

300 

1000 

- 

II 

— 

12 

II 

60 

RT-1400 

300 

1400 

- 

II 

- 

13 

II 

60 

RT-1800 

300 

1800 

- 

II 

- 

14 

It 

20 

RT-1000 

500 

1000 

-- 

II 

- 

15 

II 

40 

RT-1000 

500 

1000 

. - 

II 

- 

16 

It 

60 

RT-1000 

500 

1000 

- 

II 

- 

17 

II 

20 

RT-1400 

500 

1400 

- 

It 

- 

18 

II 

40 

RT-1400 

500 

1400 

- 

II 

- 

19 

II 

60 

RT-1400 

500 

1400 

- 

II 

- 

20 

II 

20 

RT-1800 

500 

1800 

— 

II 

-- 

21 

• 1 

40 

RT-1800 

500 

1800 

*- 

II 



22 

II 

60 

RT-1800 

500 

1800 

- 

II 

— 

23 

CYCLE- 

20 

RT-1400 

80 

1400 

20 1400-800 

It 

5 

24 

VARIABLE TEMP. 

20 

RT-1600 

80 

1600 

20 1600-800 

II 

5 

25 

II 

20 

RT-1800 

80 

1800 

20 1800-800 

II 

5 

26 

II 

20 

RT-1400 

120 

1400 

20 1400-800 

II 

5 (60) 

27 

II 

20 

RT-1600 

120 

1600 

20 1600-800 

II 

5 (60) 

28 

II 

20 

RT-1800 

120 

1800 

20 1800-800 

II 

5 (60) 

29 

II 

30 

RT-1400 

80 

1400 

30 1400-800 

II 

5 (60) 

30 

It 

30 

RT-1600 

80 

1600 

30 1600-800 

II 

5 (60) 

31 

II 

30 

RT-1800 

80 

1800 

30 1800-800 

II 

5 (60) 

32 

II 

20 

RT-1400 

80 

1400 

20 1400-800 

II 

25 (60) 

33 

It 

20 

RT-1600 

80 

1600 

20 1600-800 

II 

25 (60) 

34 

35 

II 

HOT STREAK (+50°F) 

20 

RT-1800 

80 

1800 

20 1800-800 

II 

25 (60) 

36 

UNIFORM TEMP. 

20 

RT-1000 

300 

1000 

- 

II 

- 

37 

II 

20 

RT-1400 

300 

1400 

- 

II 

- 

38 

39 

II 

HOT STREAK (+50°F) 

20 

RT-1800 

300 

1800 


II 

— 

40 

CYCLIC TEMP. 

20 

RT-1400 

80 

1400 

20 1400-800 

II 

5 (60) 

41 

II 

20 

RT-1600 

80 

1600 

20 1600-800 

II 

5 (60) 

42 

II 

20 

RT-1800 

80 

1800 

20 1800-800 

II 

5 (60) 

43 

HOT STREAK (+100° F) 20 


80 

1400 

20 1400-800 

It 

5 (60) 

44 

CYCLIC TEMP. 

20 


80 

1600 

20 1600-800 

II 

5 (60) 

45 

II 

20 


80 

1800 

20 1800-800 

II 

5 (60) 
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COMPONENTS TO BE ANALYZED AND TESTED 
IN HOST BURNER LINER CYCLIC RIGS 

1. QUARTZ LAMP BOX RIG (DELIVERY DATE - SEPTEMBER '83) 

FLAT PLATES WITH AND WITHOUT HOLES 

2. QUARTZ LAMP ANNULAR RIG (DELIVERY DATE - APRIL '84) 

CYLINDRICAL SHELLS WITH AND WITHOUT HOLES 
SUBELEMENTS OF COMBUSTOR LINERS 


FLAT PLATE SPECIMEN 

CONFIGURATIONS TO BE ANALYZED AND TESTED 



o o 




o o o 
o o o 
o o o 
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SCHEMATIC OF QUARTZ LAMP BOX RIG 


^-VIEW PORT 


■SPECIMEN COOLING AIR 


TEST 

SPECIMEN 


HEATER 

LAMPS- 



1 

-■« & 

b 



3 

g=fl ^ an Q.ga— 

=1 




WATER COOLING 


LAMP COOLING AIR 


SCHEMATIC OF QUARTZ LAMP ANNULAR RIG 


LINER COOLING AIR 


QUICK DISCONNECT 

INSTRUMENTATION 

PANELS 



LINER AIR 
AND HEATER 
DISCHARGE 
ATMOSPHERIC 
PRESSURE 
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DIMENSIONS 
ARE IN in. 



FRAME AND PLATE 

THERMOCOUPLE AND STRAIN GAGE LOCATIONS 



STRAIN GAGES 











TIME VARIATION OF STRESS IN ELEMENT 18 
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HOST LINER CYCLIC FACILITIES 


FACILITY DESCRIPTION 
Donald Schultz 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland s Ohio 44135 


The HOST Liner Cyclic Program is utilizing two types of test apparatus, 
rectangular box rigs and a full annular rig. To date two quartz lamp cyclic 
box rigs have been tested and a third is to begin testing in late October 
1983. The box rigs are used to evaluate 5x8 inch rectangular linear samples. 

A 21 inch diameter outer liner simulator is also being built up for testing 
beginning in April 1984, This annular test rig is being procured through a 
cooperative agreement with P&W-E, Hartford. Linder this agreement they supply 
the test rig and two test liners while we supply the test facility. The data 
are shared by both. All rigs are atmospheric rigs. 

The first box rig, a three 6-kVA lamp installation, was operated under 
adverse conditions to determine feasibility of using quartz lamps for cyclic 
testing. This work was done in December 1981 and looked promising. 

The second box rig, again using three 6-kVA lamps, was operated to obtain 
instrumentation durability information and initial data input to a Finite 
Element Model. This limited test program was conducted in August 1983. Five 
test plates were run. Instrumentation consisted of strain gages, thermocouples 
and thermal paint. The strain gages were found to fail at 1200° F as expected 
though plates were heated to 1700° F. 

The third box rig, containing four 6-kVA lamps, is in build-up for testing 
to beging in late October 1983. In addition to 33 percent greater power input, 
this rig has provision for 400° F backside line cooling air and a viewing port 
suitable for IR camera viewing. The casing is also water cooled for extended 
durability. 

The 21 inch diameter outer liner simulator uses 112 6-kVA lamps to 
cyclicly heat the test liner. Power levels will be adjusted to simulate 
typical linear heat loadings. Air will be supplied to provide typical 
liner film and backside cooling. In addition to the test liners being supplied 
by P&W, which are production configurations, several simple liners will be 
tested to obtain confidence in the Finite Element Model. 

This apparatus in addition to 672 kVA of 480 volt power, requires 
7.5 lb/sec of 500° F air at 35 psia, 3.5 lb/sec of ambient temperature air at 
5 psig, 1.5 lb/sec of ambient temperature air at 1 psig and 70 gal/min of 
cooling water. This apparatus is scheduled to go under testing in April 1984. 


Special Test Instrumentation 

Liner cold side temperatures will be measured using an IR-TV camera system 
with thermocouple measurements for varif i cat ion. The IR-TV camera system will 
permit several hundred temperature measurements to be made in a relatively 
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small area. Linear hot side temperatures will be measured with thin film 
thermocouples. New technology high temperature strain gauges will be used to 
obtain local strain measurements. 


Facility Preheater Test 

The proposed natural gas fired vitiated preheater for the annular test 
rig was tested in February 1983. As part of this test, quartz tubing was 
placed in the preheater exhaust stream to determine if exhaust smoke would 
deposit on the quartz. The tubing remained soot free but an iron oxide deposit 
did appear. The rust is thought to have originated in the carbon steel labo- 
ratory air supply system. As a result a filter and stainless steel pipe will 
be used in the final rig installation. Additional preheater testing in the 
actual installation is scheduled for December 1983. 
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HOST LINER CYCLIC FACILITIES 


I. QUARTZ LAMP BOX RIGS 

A. FIRST 

B. SECOND 

C. THIRD 

II. QUARTZ LAMP ANNULAR RIG 


FIRST QUARTZ LAMP BOX RIG 



CS-82-2381 
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FIRST QUARTZ LAMP BOX RIG 


•OBTAINED 2000° F TEST PLATE TEMPERATURE 
•LIMITED LAMP LIFE 


SECOND QUARTZ LAMP BOX RIG 


OBJECTIVE; 

•VERIFY PERFORMANCE AND DURABILITY OF INSTRUMENTATION 
•PROVIDE INITIAL DATA IMPUT TO FINITE ELEMENT MODEL 
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TEST HARDWARE 


•5x8 in. FLAT HASTELLOY-X PLATES 

•TWO UN INSTRUMENTED 

•ONE INSTRUMENTED WITH THERMAL PAINT 

•ONE WITH THERMAL PAINT, NINE C/A-t/c's AND ONE STRAIN GAGE 

•ONE WITH THERMAL PAINT, NINE C/A-t/c's AND EIGHT STRAIN GAGES 

•TEST PLATE SUPPORT FRAME WITH TWENTY-TWO C/A-t/c's 


TEST RIG AND HARDWARE USED ON INITIAL FLAT PLATE TESTS 




TEST RESULTS 


* SUPPORT FRAME HEATED TO 1500° F 

• SAMPLE PLATE TO 1700° F 

•STRAIN GAGES FAILED AS EXPECTED AT 1200 P F 

•STRAINS WERE CONTROLLABLE BY VARYING LAMP POWER 
AND COOLING AIR FLOW RATE 


QUARTZ LAMP BOX RIG 


SILICON WINDOW— \ 


QUARTZ 
WINDOW 
FOUR 6-kW 
QUARTZ LAMPS 


LAMP REFLECTOR PLATE-^ 



SAMPLE PLATE 

BACKSIDE SAMPLE 
COOLING AIR IN 


SAMPLE PLATE 
SUPPORT SHELF 
'SPI ^-QUARTZ WINDOW 
SUPPORT SHELF 
LAMP COOLING 
AIR IN 
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COOLING WATER 


NATURAL GAS 


40-150 psig 
AIR SUPPTY 


VITIATED 

PREHEATER 


MUFFLER 


VIEW -I 

PORTS! 













QUARTZ LAMP CYCLIC COMBUSTOR TEST RIG 


QUICK DISCONNECT 
INSTRUMENTATION 



LINER AIR 
AND HEATER 
DISCHARGE 
ATMOSPHERIC 
PRESSURE 


SPECIAL INSTRUMENTATION 

I. IR - TV MONITORING OF LINER TEMPERATURES 

II. THIN-FILM THERMOCOUPLES 

III. LASER STRAIN GAUGE 

IV. HIGH TEMPERATURE STRAIN GAUGES 
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FACILITY PREHEATER TEST 


OBJECTIVE: 

1. DETERMINE IF A NATURAL GAS FIRED VITIATED PREHEATER 
IS CLEAN ENOUGH FOR USE WITH THE ANNULAR LAMP RIG 

2. EVALUATE PERFORMANCE OF LOW PRESSURE LINER CYCLIC 
CAN RIG VITIATED PREHEATER 


TEST RESULTS 

1. RUST FROM AIR SYSTEM PIPING DEPOSITED ON QUARTZ 
TUBES AT SOtf 5 F 

2. NO SMOKE 

3. COMBUSTION EFFICIENCY OVER 95 PERCENT 
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DEPOSITS ON QUARTZ TUBING AFTER FACILITY PREHEATER TEST 



HOST CYCLIC LINER PROGRAM 


QUARTZ LAMP BOX RIGS 
FIRST 
SECOND 
THIRD 

QUARTZ LAMP ANNULAR RIG 
PREHEATER TESTING 


SCHEDULE 


FY 82 

FY83 

FY84 

FY85 












FY86 


DESIGN, PROCUREMENT, INSTILLATION 
TEST 
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LIFE PREDICTION AND CONSTITUTIVE BEHAVIOR 


G. R. Halford 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


One of the primary drivers that prompted the initiation of the HOST Pro- 
gram was the recognized need for improved cyclic durability of costly hot 
section components. All too frequently, fatigue in one form or another was 
directly responsible for the less than desired durability, and prospects for 
the furture weren't going to improve unless a significant effort was mounted 
to increase our knowledge and understanding of the elements governing cyclic 
crack initiation and propagation lifetime. Certainly one of the important 
factors is the ability to perform accurate structural stress-strain analyses 
on a routine basis to determine the magnitudes of the localized stresses and 
strains since it is these localized conditions that govern the initiation and 
crack growth processes. Developing the ability to more accurately predict 
crack initiation lifetimes and cyclic crack growth rates for the complex load- 
ing conditions found in turbine engine hot sections is of course the ultimate 
goal of the life prediction research efforts. 

It has been found convenient to divide the research efforts into those 
dealing with nominally isotropic and anisotropic alloys; the latter for appli- 
cation to directionally solidified and single crystal turbine blades. 

I do want to emphasize the underlying thrust of these programs: The 

development and verification of WORKABLE engineering methods for the calcula- 
tion, in advance of service, of the local cyclic stress-strain response at the 
critical life governing location in hot section components, and the resultant 
cyclic crack initiation and crack growth lifetimes. 

In attacking the problems, we are utilizing the talents of industry, 
universities, and the government. Figures 1 to 4 display the specific pro- 
grams that have been initiated and funded by the HOST Project Office. Those 
programs that have been on-going long enough to permit progress to be reported 
will be presented by the respective principal investigators. The remaining 
programs will be briefly touched upon either by the NASA Lewis Technical Pro- 
gram Manager or by myself. 
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CRACK INITIATION LIFE PREDICTION AND CONSTITUTIVE MODELING 

(NOMIMALLY ISOTROPIC MATERIALS) 

9 CREEP-FATIGUE LIFE PREDICTION FOR ENGINE HOT SECTION MATERIALS (ISOTROPIC) 
G. HALFORD. LERC TECHNICAL PROGRAM MANAGER 
NAS 3-23288, PRATT & WHITNEY AIRCRAFT 
V. MORENO, PRINCIPAL INVESTIGATOR 
STARTED JUNE 1982 

2 YEAR BASE PROGRAM/3 YEAR OPTIONAL PROGRAM 

• CONSTITUTIVE MODELING FOR ISOTROPIC MATERIALS 
A. KAUFMAN, LERC TECHNICAL PROGRAM MANAGER 
DUAL AWARDED PROGRAM 
NAS 3-23927, GENERAL ELECTRIC COMPANY 
J. LAFLEN, PRINCIPAL INVESTIGATOR 
STARTED MAY 1983 

2 YEAR BASE PROGRAM/3 YEAR OPTIONAL PROGRAM 

NAS 3-23925, SOUTHWEST RESEARCH INSTITUTE 
U. LINDHOLM, PRINCIPAL INVESTIGATOR 
STARTED MAY 1983 

2 YEAR BASE PROGRAM/3 YEAR OPTIONAL PROGRAM 

IN-HOUSE HOST PROGRAMS 

9 FATIGUE LABORATORY UPGRADING PROGRAM 
M, MCGAW, LERC PROGRAM MANAGER 
MULTIYEAR EFFORT 

- REMODELING OF LABORATORY 

- MASTER/SLAVE COMPUTER SYSTEM 

- BIAXIAL FATIGUE FACILITY 

- LCF/HCF FACILITY 

9 SINGLE CRYSTAL DURABILITY PROGRAM 
G. HALFORD, LERC PROGRAM MANAGER 
MULTIYEAR EFFORT 

JOINT EFFORT WITH MATERIALS DIVISION 
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CRACK PROPAGATION LIFE PREDICTION MODELING 


* HIGH TEMPERATURE CRACK PROPAGATION 

T. ORANGE, LERC TECHNICAL PROGRAM MANAGER 


2 YEAR BASE PROGRAM/2 YEAR OPTIONAL PROGRAM 

e CRACK GROWTH MECHANISMS 

J. SHANNON, LERC TECHNICAL PROGRAM MANAGER 
NAG 3-348, SYRACUSE UNIVERSITY 
PROF. H. LIU, PRINCIPAL INVESTIGATOR 
STARTED OCTOBER 1982 

3 YEAR GRANT PROGRAM 

CRACK INITIATION LIFE PREDICTION AND CONSTITUTIVE MODELING 

(ANISOTROPIC MATERIALS) 

• LIFE PREDICTION AND MATERIAL CONSTITUTIVE BEHAVIOR FOR ANISOTROPIC MATERIALS 
R. BILL, LERC TECHNICAL PROGRAM MANAGER 


5 YEAR BASE PROGRAM/2 OPTIONAL PROGRAMS PARALLEL WITH BASE 

9 THEORETICAL CONSTITUTIVE MODELING OF ANISOTROPIC MATERIALS 
R. THOMPSON, LERC TECHNICAL PROGRAM MANAGER 
DUAL AWARDED PROGRAM 
NAG 3-XXX, UNIVERSITY OF CINCINNATI 
PROF. D. STOUFFER, PRINCIPAL INVESTIGATOR 
STARTED AUGUST 1983 
3 YEAR GRANT PROGRAM 

NAG 3-XXX, UNIVERSITY OF CONNECTICUT 
PROF. E. JORDAN, PRINCIPAL INVESTIGATOR 
STARTED SEPTEMBER 1983 
3 YEAR GRANT PROGRAM 
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CREEP FATIGUE LIFE PREDICTION FOR ENGINE HOT SECTION 
MATERIALS (ISOTROPIC)* 

Vito Moreno 

Pratt & Whitney Engineering 
United Technologies Corporation 


The presentation will summarize the activities performed during the first 
year of the NASA HOST Programs "Creep Fatigue Life Prediction for Engine Hot 
Section Materials (Isotropic)", being conducted by Pratt & Whitney Aircraft. 
The program is a 5-year, two part effort aimed at improving the high 
temperature crack initiation prediction technology for gas turbine hot section 
components. The two-year base program comprises the following tasks: 

Task I - Material/Coating/ Component Selection and Acquisition 

Task II - Screen Candidate Life Prediction Approaches 

Task III - Evaluate Best Candidate Life Prediction Approach 

Task IV - Reporting 

Significant results of the program produced thus far are listed below. 

Task I - Materia 1/Coating/ Component Selection and Acquisition 

1. Cast B1900 + Hf and wrought IN 718 were selected as the base and 
alternate materials, respectively. 

2. A single heat of B1900 + Hf was obtained and test specimens 
fabricated. 

3. The material was characterized with respect to grain size, y' size, 
carbide distribution, and dislocation density. 

4. Monotonic tensile and creep testing has shown engineering properties 
within anticipated scatter for this material. 

5. Examination of the tensile tests has shown a transition from 
inhomogeneous "planar" slip within the grains at lower temperatures 
to more homogeneous matrix deformation. 

6. Examination of the creep tests has shown a transgranular failure mode 
at 1400°F and an intergranular failure mode at 1600°F and 1800°F. 


* NASA Contract NAS3-23288 
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Task II - Screen Candidate Life Prediction Approaches 

1. A study was conducted to investigate the effects of test specimen 
geometry and fabrication process on fatigue life. As a result, axial 
strain controlled specimens were designed with a smooth (no 
extensometer ridges) gage section and fabricated using centerless 
grinding followed by light electropolishing. 

2. A fatigue test matrix was established to provide baseline data to 
define crack initiation life as a function of major variables and for 
life prediction model evaluation. A total of 43 fully reversed strain 
controlled fatigue tests have been completed. Major variables 
investigated were temperature [871°C(1600°F) vs 538°C(1000°F)], 
strain range and strain rate. 

3. Examination of specimens during testing indicated that measurable 
76mm((.030 in.) surface cracks appear early in the specimen life 
i.e. , 15% of total life at 871°C(1600°F) and 50% of life at 
538°C(1000°F) . This has been used as the definition of crack 
initiation for the initial model evaluation work. 

4. Observed crack initiation sites are all surface initiated and 
associated with either grain boundary carbides or local porosity. The 
initiation life is not significantly affected by the character of the 
site. 

5. Transgrannular cracking is observed at the initiation site for all 
conditions tested. 

6. Grain dislocation structure is significantly less than that observed 
in monotonic tensile or creep tests. 

7. A ranking procedure for evaluation of the prediction models has been 
established. The procedure assigns a numerical score based on the 
amount of data required, the predictive capability and the 
adaptability to engine relevant loading conditions for each model 
considered. 

Task III - Evaluate Best Candidate Life Prediction Approach 

1. Life prediction models representative of macroscopic (Coffin-Manson) 
and microscopic (Damage-Rate) approaches were selected for 
preliminary evaluation using the model ranking procedure developed in 
Task II. 

2. Using limited data obtained at 871°C(1600°F)-(2 strain rates) and 
538°C(1000°F)-(1 strain rate), the macroscopic approach obtained a 
higher overall score on the basis of data requirements and predictive 
capabil ity. 
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CONSTITUTIVE MODELING FOR ISOTROPIC MATERIALS 


A. Kaufman 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 


Under HOST Contract NAS3-23288, crack initiation life prediction methods 
will be developed for hot section components fabricated from isotropic mate- 
rials. To apply these methods it is first necessary to determine the compo- 
nent structural response, specifically the stress-strain history at the 
critical cracking location. The structural analysis method must be capable of 
accounting for cyclic thermomechanical loading, plastic flow during thermal 
transients, creep and stress relaxation during steady-state operation and 
inelastic strain ratchetting and reversal due to repeated flight cycles. 

In recent years, nonlinear finite element computer codes such as ANSYS 
and MARC have become available for cyclic analysis of components involving 
inelastic strains. These codes are based on classical plasticity theory and 
use separate creep constitutive models. The classical methods utilize simpli- 
fying assumptions for computational convenience. Among these assumptions are 
the definition of a specific yield surface with associated flow rules and 
hardening models and the uncoupling of time-independent (plastic) and time- 
dependent (creep) inelastic strain effects. That these classical methods and 
their assumptions do not realistically represent superalloy material behavior 
under cyclic loading have been demonstrated in two pre-HOST programs (the 
turbine blade durability study reported in NASA CR-165268 and the combustor 
liner durability study reported in NASA CR-165250). 

The objective of this program is to develop a unified constitutive model 
for finite-element structural analysis of turbine engine hot section compo- 
nents. This effort constitutes a different approach for nonlinear finite- 
element computer codes which have heretofore been based on classical inelastic 
methods. A unified constitutive theory will avoid the simplifying assumptions 
of classical theory and should more accurately represent the behavior of 
superalloy materials under cyclic loading conditions and high temperature 
environments. Model development will be directed toward isotropic, cast 
nickel-base alloys used for aircooled turbine blades and vanes. The Con- 
tractor will select a Base Material for model development and an Alternate 
Material for verification purposes from a list of three alloys specified by 
NASA. The candidate alloys represent a cross-section of turbine blade and 
vane materials of interest to both large and small size engine manufacturers. 
Material stock for the Base and Alternate Materials will be supplied to the 
Contractor by the Government. 

The contractual effort will be conducted in two phases, a Basic Program 
of two years duration and an optical follow-on program also of two years dura- 
tion. In the Basic Program, a unified constitutive model will be developed 
for the prediction of the structural response of isotropic materials for the 
temperatures and strain ranges characteristic of cooled turbine vanes in ad- 
vanced gas turbine engines. The constitutive model will be implemented in an 
existing finite-element computer code or a new code specifically developed for 
unified models. An evaluation will be made of the capability of the analyti- 
cal method to predict structural response for multiaxial stress states and 
nonisothermal conditions by conducting thermomechanical loading and benchmark 
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notch verification experiments and analyses. As a final evaluation of the 
analytical methods, a structural analysis will be performed for a hot section 
component fabricated of the Base Material for simulated engine operating 
conditions. In the Optical Program further development will be undertaken to 
consider thermal history effects and to correct any deficiencies indicated in 
the constitutive model or in the computational algorithms in the code. In 
addition, the constitutive model development will be verified for an alternate 
materi al . 

This program will be conducted under a dual effort by the General Elec- 
tric Company and Southwest Research Institute. The principal investigators 
are Dr. J. H. Laflen for General Electric and Dr. U. S. Lindholm for Southwest 
Research Institute. Successful completion of this program will result in a 
significant advance in structural analysis capabilities for gas turbine engine 
hot section components. These improvements in the structural analysis tools 
should result in better life prediction capability and greater hot section 
durability for the next generation of gas turbine engines. 


212 



CONSTITUTIVE MODEL REQUIREMENTS 


• AVOID SIMPLIFICATIONS OF CLASSICAL THEORY 
(YIELD SURFACE, UNCOUPLED CREEP-PLASTIC STRAINS) 

• APPLICABLE TO STRAIN-TEMPERATURE RANGES AT 
CRITICAL LOCATIONS OF TURBINE BLADES AND VANES 

• CAPABLE OF REPRESENTING INELASTIC RESPONSE DURING 
TRANSIENT AND STEADY-STATE ENGINE OPERATION 

• REASONABLE MATERIAL TESTING REQUIREMENTS 

• PRACTICAL FOR IMPLEMENTATION IN NONLINEAR 
STRUCTURAL, FINITE-ELEMENT COMPUTER CODES 


CONSTITUTIVE MODELING FOR ISOTROPIC MATERIALS 
OBJECTIVE : 

TO DEVELOP A UNIFIED CONSTITUTIVE MODEL FOR REPRESENTING 
CYCLIC INELASTIC BEHAVIOR OF ISOTROPIC CAST NICKEL-BASE 
ALLOYS USED FOR AIRCOOl£D GAS TURBINE BLADES AND VANES 
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CONSTITUTIVE MODELING FOR ISOTROPIC MATERIALS 


R.H. Van Stone 
L.T, Dame 
R.L. McKnight 
J.H. Laflen 

Aircraft Engine Business Group 
General Electric Company 

The inelastic finite element analytic model that is chosen for a particular 
type of numerical analysis of material behavior represents a compromise between 
(1) the requirements of physical versimilitude , (2) mathematical accuracy and 
stability, and (3) computational convenience and economy. In the past, the mathe- 
matical and computational problems have been so dominant that most analysts have 
been content to adopt simple uncoupled material models. Plasticity has been 
represented by an engineering stress-strain curve, monotonic or cyclic, and creep 
by a relatively uncomplicated power law equation. 

It is well known that the simplest uncoupled material models ignored many 
of the more complex, well documented, types of high temperature material behavior: 
inelastic recovery, cyclic creep, strain rate effects, and thermomechanical 
(simultaneous temperature and load variation) effects on material response. Recent 
work has evolved constitutive theories which contain representations for various 
combinations of these effects. Frequently these models do not make the classical 
assumption which separates the analytical treatment of time dependent and time 
independent inelastic strains. These methods are generally referred to as unified 
theories. It will be the purpose of this program to thoroughly evaluate such 
methods for application to typical isotropic cast nickel base superalloys used for 
air-cooled turbine blades and vanes. 

This objective will be accomplished through a two year combined analytical 
and experimental program. During this nine-task program, a nonlinear finite 
element program will be developed which contains at least two selected constitu- 
tive theories. These models will be evaluated using smooth thin-walled uniaxial 
Rene* 80 specimens. These experiments will be selected for both determining 
material constants and for verification purposes . Further experimental work on 
more complicated specimen geometries will be used to select the most promising 
constitutive model and to develop efficient numerical algorithims for the nonlinear 
finite element code. These additional experiments will involve multiaxial states 
of stress and a notch specimen geometry subjected to a variety of loading condi- 
tions. The finite element computer program will be additionally demonstrated by 
modeling an actual hot path component such a turbine blade or vane. 

Work on this program commenced in May of 1983, and as such most of the 
efforts have gone toward defining the experimental plan. This includes the 
selection of the base (Rene' 80) and alternate (Mar-M-247) materials; development 
of the test matrix for all of the base program tasks; and selection of the 
specimen geometries for the uniaxial, multiaxial, and notched specimen experiments. 
Various unified constitutive theories which were identified through a literature 
search are currently under review. A series of criteria have been developed for 
selecting the preferred theories which will be further investigated during later 
tasks . 
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CONSTITUTIVE MODELING FOR ISOTROPIC MATERIALS* 


Ulric S. Lindholm 
Southwest Research Institute 
San Antonio, TX 78284 


The objective of this program is to develop a unified constitutive 
model for finite-element structural analysis of turbine engine hot section 
components. This effort constitutes a different approach for nonlinear 
finite-element computer codes which have heretofore been based on classical 
inelastic methods. A unified constitutive theory will avoid the simplifying 
assumptions of classical theory and should more accurately represent the 
behavior of superalloy materials under cyclic loading conditions and high 
temperature environments. Model development will be directed toward iso- 
tropic, cast nickel -base alloys used for air-cooled turbine blades and vanes. 
Recent studies have shown that this approach to the modeling of material 
behavior is particularly suited for determining the cyclic behavior of 
superalloy type blade and vane materials and is entirely compatible with 
three dimensional inelastic finite element formulations. More efficient 
and accurate inelastic analysis of hot section components - turbine blades, 
turbine vanes, combustor liners and seals - fabricated from "age hardenable" 
isotropic superalloy materials will be realized as the result of this 
activity. The final constitutive model, programmed in a stand-alone module 
suitable for incorporation in a general purpose structural analysis program, 
will be provided by Southwest Research Institute (SwRI) with the assistance 
of Pratt and Whitney Aircraft (PWA) as a principal subcontractor. 

The program will be conducted in two phases. A basic program (Tasks A 
through I) and an optional follow-on program (Tasks J through M) . In the 
Basic Program of twenty six months' duration a unified constitutive model will 
be developed for the prediction of the structural response of isotropic 
materials for temperatures and strain range characteristics of cooled turbine 
vanes in advanced gas turbine engines (Task A). A data base of uniaxial and 
multi axial material properties required for the constitutive model develop- 
ment will be obtained for the base material (Tasks C and E). The constitutive 
model will then be incorporated into a finite-element computer code (Task D). 
An evaluation will be made of the capability of the analytical method to 
predict the structural response for multi axial stress states (Task E) and 
nonisothermal conditions by conducting thermomechanical loading and bench- 
mark notch verification experiments and analysis (Task F). As a final evalu- 
ation of the analytical model, a structural analysis will be performed for 
a hot section component fabricated from the base material for simulated engine 
operating conditions (Task G). In the optional program material property 
test procedures will be further developed to minimize the amount of testing 
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required, and to study the possibility for estimating the material model 
constants from conventional property data (Task J). Further development 
of the model will be undertaken to consider thermal history effects and 
to correct any deficiencies indicated in the model or the computational 
algorithms in the code (Task K). In addition the constitutive model 
development will be verified for an alternate material (Task L). 

Technical progress to date has concentrated on a review and screening 
of candidate constitutive models (Task A) and the fabrication of test speci- 
mens (Task B). The discussion below will concentrate on the formulation of 
unified constitutive theories. Test specimen fabrication is underway with 
the two alloys selected for study; PWA Alloy B1900+Hf and MAR M-247. 


The unified theories considered will be based on two primary premises: 
First, that the total strain or strain rate may be partitioned into two 
components, 




( 1 ) 


where all inelastic deformation is included in e. - p which is additive with 
the elastic strain rate. There is no further partitioning, e.g., between 
plasticity and creep. A second important premise is that both components 
of strain rate are continuous, non- zero functions of stress; i.e., there is 
no partitioning of stress space between purely elastic and elastic-plastic 
or elastic-viscoplastic regions as is typical of classical theories in 
plasticity. The latter premise results in computational efficiency since 
a single set of continuous equations can be integrated throughout loading, 
unloading, reverse loading and all other complex histories including strain 
rate and temperature variations. The classical yield criteria (plastic 
potential) and associated flow rule are replaced by evolutionary equations 
for a finite set of internal variables describing the progression of work 
hardening and recovery processes. 


Most theories utilize one of four variations on the basic Prandtl- 


and 


These are 



e.. p = 
1J 

A 3f 

1 

(2) 

• P 
£ • • = 
1J 

X 2 ^ s i j " n ij* 

(3) 

• P 
£ . . = 
1J 

X 3 s ij 

(4) 

£.- P = 

^jkl s kl 

(5) 


Eq. 2 is in the form utilizing a flow potential f = f(s ii5 K, T) with a 
limiting surface, f = constant, bounding elastic behavior. Eqs. 3-5 are unified 
theories without a yield surface. In Eqs. 3 and 4, A is a scalar variable 
encorporating both isotropic and anisotropic variables. Most theories also 
employ the kinematic hardening tensor ft j j , often referred to as a "back" stress 
or "equilibrium" stress. An alternate approach to describe anisotropic harden- 
ing is the generalization of the isotropic flow law, Eq. 4, to the anisotropic 
case, Eq. 5. This latter form is relatively undeveloped. 
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A second basic feature of most constitutive models is the functional 
form selected for the kinematic relationship between plastic strain rate and 
stress or their invarients, D?P = e-j jP e-jjP,and J 2 = (s-jj-ftij) (s-jj-Oji). 
Three typical examples found in the literature are given below and plotted 
in Figure 1. 



= A (J 2 /K 2 ) 
= D q 2 exp 



= B 


si nh (J 2 /K 2 ) m 


( 6 ) 

(7) 

( 8 ) 


where K is a hardening variable and A, B, D , m and n are material constants 
(possibly temperature dependent). The values of m and n are the dominant 
variables governing strain rate sensitivity. 


An additional general concern is the choice of a tracking parameter 
to follow hardening and recovery processes. This parameter must include 
history dependence. The two most commonly used are / ePdt and / Wpdt, 
integrals over the plastic strain rate or the plastic work rate history. 
Differences lie in whether hardening occurs in the direction of the strain 
increment or the stress increment. 


Comparison and evaluation of existing models is now underway with the 
objective of clarifying limitations and recommending improvements. Some 
examples of computations using unified theory for cyclic uniaxial and biaxial 
strain histories are given in Figures 2 and 3. Biaxial, nonproportional strain 
paths are particularly critical to the development of the evolutionary equations 
for hardening and recovery. 



FIGURE 1. A COMPARISON OF THREE BASIC FORMS OF THE KINETIC EQUATION. 
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FIGURE 2. CYCLIC HARDENING RESPONSE FOR HASTELLOY X 
AT 1600F COMPUTED WITH UNIFIED THEORY. 




FIGURE 3. BIAXIAL HARDENING OF HASTELLOY X AT ROOM TEMPERATURE 
UNDER 90° OUT-OF-PHASE TENSION-TORSION STRAIN 
CYCLING, EXPERIMENT (LEFT) AND UNIFIED THEORY (RIGHT). 
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LIFE PREDICTION AND CONSTITUTIVE MODELS 
FOR ANISOTROPIC MATERIALS 


Robert C. Bill 
Propulsion Laboratory 

AVRADCOM Research and Technology Laboratories 
Lewis Research Center 
Cleveland, Ohio 

The trend toward improved engine efficiency and durability is placing 
increased demands on gas turbine materials, especially in the hot section. 

New materials and coatings are being developed to meet these demands. Turbine 
airfoil components present a particular challenge because they must survive in 
an especially hostile thermal, mechanical, and chemical environment. Here 
single crystal (SC) and directionally solidified or recrystallized (DSR) 
polycrystall ine materials are finding application because of their inherently 
superior resistance to high temperature deformation and failure, compared to 
conventionally cast materials. Difficulties impeding the full implementation 
of SC or DSR materials include the limited knowledge and understanding of 
failure (crack initiation) mechanisms and constitutive behavior. 

The objective of this program is to develop practical life prediction and 
constitutive behavior models for anisotropic SC and DSR materials employed in 
turbine airfoils. Though emphasis will be placed on the development of widely 
applicable models (not material specific), experimental work will focus 
primarily on SC and DSR materials that are likely candidates for service in 
the near future. 

The program is divided into a base program and two options. The base 
program is concerned with single crystal materials under conditions pertinent 
to first stage turbine airfoil applications. Coating effects on airfoil 
durability and constitutive response are an important part of the base 
program. Option one will address the problem of single crystal blade root 
attachment, including the effects of attachment geometric features. 
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temperatures in the blade root area, and the type of loading encountered. The 
second option includes airfoil and blade root attachment considerations for 
DSR materials. 

Upon completion of the program, generally applicable analytical models 
describing the time dependent cyclic stress-strain response, and models needed 
for life prediction and durability assessment of anisotropic materials will be 
available to the gas turbine design community. These models will be 
incorporated into a sophisticated structural analysis code available to the 
public. 
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CONTRACT; 

CYCLIC CONSTITUTIVE MODELING AND LIFE PREDICTION METHODS FOR 
ANISOTROPIC MATERIALS 

OBJECTIVE: 

DEVELOP AND VERIFY CYCLIC MATERIALS CONSTITUTIVE MODELS AND 
LIFE PREDICTION METHODS FOR COMPONENT SPECIFIC ANISOTROPIC 
MATERIALS FOR USE IN STRUCTURAL ANALYSIS COMPUTER PROGRAMS 

DURATION: 

5 YEAR, 35 MAN-YEAR EFFORT 


PROGRAM STRUCTURE: THREE PHASE PROGRAM 

•BASE PROGRAM - COATED SC AIRFOILS 
•OPTION 1 - UNCOATED SC BLADE ROOTS 

•OPTION 2 - DSR MATERIALS 


APPROACH: 


•SELECT MATERIALS AND COATINGS 

•SCREEN ADVANCED CONSTITUTIVE AND LIFE 
PREDICTION MODELS ; SELECT BEST 

•INTEGRATE MODELS WITH STRUCTURAL ANALYSIS 
PROGRAMS 

•VERIFY IN SIMULATED COMPONENT TEST 


PROGRAM WILL ADDRESS: 

•CREEP-FATIGUE 
•THERMOMECHANICAL FATIGUE 
•ORIENTATION EFFECTS 
•COATING/SUBSTRATE INTERACTIONS 
•BIAXIAL LOADING 

•ATTACHMENT STRESS CONCENTRATIONS 
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BASE PROGRAM 


LIFE PREDICTION AND CONSTITUTIVE MODELS FOR COATED 
SINGLE CRYSTAL MATERIALS AT AIRFOIL TEMPERATURES 

1 

2 

3 

4 

5 

TASK I. MATERIAL/COATING SELECTION AND ACQUISITION 

TASK II. SELECTION OF CANDIDATE LIFE PREDICTION AND 

CONSITITUTIVE MODELS 

TASK III. LEVEL I EXPERIMENTS FOR SC MATERIALS 

EXERCISE OF OPTION 1 

TASK IV. CORRELATION OF MODELS WITH LEVEL 1 SC EXPERIMENTS 

TASK V. LEVEL 2 SC EXPERIMENTS 

PROVISIONAL EXERCISE OF OPTION 2 
TASK VI. FINAL SELECTION OF LIFE PREDICTION AND CONSTITUTIVE 

MODELS 

TASK VII. SUBCOMPONENT VERIFICATION FOR PRIMARY SC MATERIAL 

TASK VIII. ALTERNATE SC MATERIAL CHARACTERIZATION FOR AIRFOIL 
APPLICATIONS 

TASK IX. MODEL VERIFICATION ON ALTERNATE SC MATERIAL 

TASK X. DELIVERY OF COMPUTER CODE TO NASA 

TASK XI. REPORTING REQUIREMENTS 

A L 
AAAl 

A 

iAAAi 


~A 

A 

A 
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OPTION 1 


LIFE PREDICTION AND CONSTITUTIVE MODELS FOR UNCOATED 
SC MATERIALS AT ROOT ATTACHMENT TEMPERATURES 

1 

2 


■ 

5 

TASK XII. SPECIMEN PREPARATION 

TASK XIII. SELECTION OF CANDIDATE CONSTITUTIVE AND LIFE PREDICTION 

MODELS FOR UNCOATED SC MATERIALS AT ROOT ATTACHMENT 
TEMPERATURES 

TASK XIV. CYCLIC LIFE AND CONSTITUTIVE BEHAVIOR OF UNCOATED SC 

MATERIALS AT ROOT ATTACHMENT TEMPERATURES 

EXERCISE OF OPTION 2 

TASK XV. FINAL SELECTION OF CONSTITUTIVE MODEL AND LIFE PREDICTION 

MODEL FOR SC MATERIAL AT ROOT ATTACHMENT TEMPERA- 
TURES 

TASK XVI. MODEL VERIFICATION ON PRIMARY SC MATERIAL FOR BLADE 

ROOT ATTACHMENT 

TASK XVII. ALTERNATE SC MATERIAL CHARACTERIZATION FOR BLADE ROOT 

ATTACHMENT 

TASK XVIII. MODEL VERIFICATION ON ALTERNATE SC MATERIAL FOR BLADE 
ROOT ATTACHMENT 

TASK XIX. DELIVERY OF COMPUTER CODE TO NASA 

TASK XX. REPORTING REQUIREMENTS 


- 

A 

A L 

— A 
A 

.AAA 
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OPTION 2 


LIFE PREDICTION AND CONSTITUTIVE MODELS FOR DSR MATERIALS 

1 

Z 

3 

4 

5 

TASK XXI. DSR SPECIMEN PREPARATION AND MATERIAL CHARACTERIZATION 

TASK XXII. SELECTION OF CANDIDATE LIFE PREDICTION AND CONSTITUTIVE 

MODELS FOR DSR MATERIALS 

TASK XXIII. CREEP-FATIGUE EXPERIMENTS FOR COATED DSR MATERIAL AT 

AIRFOIL TEMPERATURES 

TASK XXIV. STRESS-STRAIN EXPERIMENTS FOR UNCOATED DSR MATERIALS AT 

ROOT ATTACHMENT TEMPERATURES 

TASK XXV. FINAL SELECTION OF LIFE PREDICTION AND CONSTITUTIVE MODELS 

FOR DSR MATERIALS 

TASK XXVI. MODEL VERIFICATION FOR DSR MATERIALS AT AIRFOIL TEMPERA- 

TURES 

TASK XXVII. MODEL VERIFICATION FOR UNCOATED DSR MATERIALS AT ROOT 

ATTACHMENT TEMPERATURES 

TASK XXVIII. DELIVERY OF DSR MATERIAL COMPUTER CODES TO NASA 

TASK XXIX. REPORTING REQUIREMENTS 
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Surface Protection Overview 


Stanley R. Levine 

National Aeronautics and Space Administration 
Lewis Research Center 


Among the factors that must be considered in the design of turbine engine hot 
section components is the ability of the materials to withstand the environment at 
the design conditions of life, temperature and loads. In some instances the design 
envelope is governed by the environmental resistance of available alloys and coatings. 
The environmental attack processes of concern, oxidation and hot corrosion are illu- 
strated schematically in Figure 1. The progress of these processes can by acceler- 
ated by particulate erosion. In the oxidation process, the less noble alloy consti- 
tuents react to form oxides. These oxide scales may then spall as a result of ther- 
mal cycling. To combat the process, slow growing, adherent, non-volatile oxide scales 
such as alumina, AI 2 O 3 , are most desirable. In the hot corrosion process condensed 
salts formed from air or fuel impurities cause accelerated attack by penetration of 
the oxide scale via defects or by fluxing the oxide scale. Again, the best defense 
is the formation of a slow-growing, adherent scale that is not reactive with the con- 
densed salt. As illustrated in Figure 2, the hot corrosion process can by more des- 
structive than the oxidation process in spite of the fact that hot corrosion gener- 
ally occurs below temperatures at which oxidation is a concern. The environmental 
attack processes can by quite destructive as illustrated in Figure 3. Here, the first 
stage vanes of a turboprop aircraft have been severely hot corrosion attacked. At 
earlier stages of this attack aerodynamic efficiency was no doubt impaired. Environ- 
mental attack can also disrupt the effectiveness of airfoil cooling designs by 
changing the dimensions or plugging passages via deposition or scale formation. In 
gas turbines where the gas temperature can exceed the alloy melting point, the con- 
sequences of cooling system breakdown can be severe. 

Environmental attack can also be insidious, but equally destructive. For example, 
oxidation or hot corrosion can play a significant role in crack initiation and growth. 
Finally, the design envelope expanding metallic or ceramic coatings display mechani- 
cal properties that are very different from the structural alloys they protect. Thus, 
the mechanical behavior of these coatings, particularly with regard to creep and crack 
initiation and growth, must be considered by the designer. 

To summarize, to adequately predict the consequences of environmental attack we 
must be able to predict salt deposition, the kinetics of oxidation and hot corrosion 
attack, the depletion of coatings via reaction with substrates and the effects of 
environmental attack and coatings on system mechanical response. 

The objective of the HOST surface protection effort is to develop a first-cut 
integrated environmental attack life prediction methodology for hot section components. 
The program elements and the expected results are described in Figure 4. Environmen- 
tal/mechanical property interaction is a joint effort involving the Fatigue and Frac- 
ture, Surface Protection, and Metallic Materials Branches. The airfoil deposition 
model effort resides in the High-Temperature Chemistry Section while the coating life 
prediction effort is primarily concentrated in the Surface Protection Section, both 
part of the Materials Durability Branch. 
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Further details of the surface protection effort are given in Figure 5. Those 
efforts which have been active through most of the previous year are reported by 
the principal investigators in this workshop. The primary new efforts for the coming 
year are the initiation of thermal barrier coating life prediction and dual cycle 
attack and procurement initiation for metallic coating life prediction verification. 
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EXPECTED RESULTS 


MODEL FOR EFFECTS OF ENVIR. 
ATTACKS COATINGS ON CRACK 
INITIATION 

MODEL TO PREDICT THE LOCATION & 
POTENTIAL SEVERITY OF CORROSION/ 
ATTACK OF TURBINE AIRFOILS 

CAPABILITY TO PREDICT COATING 
DEGRADATION ON BLADES, VANES, 
COMBUSTORS 


231 















Page intentionally left blank 


TURBINE AIRFOIL DEPOSITION MODELS* 


Daniel E. Rosner 
Yale University 

Department of Chemical Engineering 

1. INTRODUCTION, OUTLINE 

Gas turbine failures associated with sea-salt ingestion and sulfur-containing fuel 
impurities have directed attention to alkali sulfate deposition and the associated "hot 
corrosion" of gas turbine (GT) blades under some GT operating conditions. These salts deposit 
and form thin molten films which can undermine the protective metal oxide coating normally 
found on GT blades. This research project (fig. 1) deals with the prediction of molten salt 
deposition, flow and oxide dissolution, and their effects on the lifetime of critical turbine 
blades. Our goals include rationalizing and predicting corrosion patterns on operational GT 
rotors and stators, and, ultimately providing some of the "tools" required to design 
laboratory simulators and future corrosion-resistant high-performance engines. 

2. DEPOSITION RATE THEORY 

A necessary step in the complex sequence of events leading to hot corrosion failures is 
the deposition of the alkali salts from the combustion products onto turbine stator and rotor 
blades. A comprehensive but tractable method for predicting the rates of chemical vapor 
deposition (CVD) of, say, alkali sulfates from multicomponent salt-laden combustion products 
has accordingly been developed, illustrated, and tested (refs. 4,5,8), The theory predicts 
important effects of mu 1 t ic amp onent diffusion, and thermal (Soret) diffusion (ref. 6) on dew 
points, CVD-rates and deposit compositions (ref. 7). For our present purposes, an important 
feature of the deposition rate is its dependence on surface temperature, shown in figure 2 for 
the case of K^SO^ vapor deposition (ref. 1). Note the existence of a "dew point" temperature 
(above which a macroscopic surface condensate is not stable) and a temperature-insensitive 
"plateau" deposition rate at temperatures sufficiently far below the dew point — i.e.« 
features included in the abovementioned CVD-rate theory (see, e.g,, refs. 4,5,7,8), While in 
its most general form allowance is made for the transport of each chemical species (i = 1,2, 
...) containing the elements of interest (Na, S, ...) across the diffusion boundary layer (BL) 
(fig. 3), in order to study the dynamics of the resulting molten salt condensate layer 
(sections 3,4, below) we have exploited the single-component limiting form of the CVD theory 
to economically provide the spatial distribution of condensate "arrival" (see fig. 8, below). 
However, our approach also applies to the deposition rate of fine particles (refs. 9,10) — 
i.e., mist droplets small enough to behave like "heavy molecules" in the prevailing carrier 
gas flow field. 


3. CONDENSATE LAYER DYNAMICS 

An important process in the molten salt attack of GT blades is evidently the dissolution 
("fluxing") of the normally protective metal oxide, and this, in turn, implies dynamic 
processes which allow dissolution/reaction to occur without "saturating" the molten salt 
"solvent." For these reasons, we have initiated studies of the dynamics of thin condensate 
liquid layers (ref. 11), allowing for the interplay of arrival rate (Section 2) and "run-off" 
induced by aerodynamic shear stress,^ (fig. 4) and blade rotation rate, A (fig. 5). These 


*NAG3~201 : "Theory of Mass Transfer from Combustion Gases" 
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studies should ultimately prove useful in the interpretation of observed hot corrosion ! Wps” 
on blades removed from operational engines (Section 4). They have already guided the 
development and interpretation of laboratory "simulator" experiments (eg., refs, 3,4,8) in 
which molten salt deposition and runoff necessarily occur. 

In reference 11 the necessary liquid layer theory was developed and used to predict 
steady-state laminar condensate layers on smooth non-rotating, isothermal targets, with 
emphasis on the circular cylinder in high Reynolds number cross-flow. In our current studies, 
this approach has been extended to include the treatment of condensate layer flow along smooth 
non-isothermal rotor blades (fig. 5). Illustrative calculations are being made for a test 
turbine using the aerodynamic design data provided by Haas and Kofskey in NASA TP 1018 (1977). 
From the inviscid stream velocity data (fig, 6 ) we compute, using efficient integral methods, 
the corresponding distributions of gas-side momentum? heat- and mass- transfer coefficients, and 
the distribution of blade "recovery" temperature (used to estimate this root-cooled blade 
temperature distribution (shown in fig, 7 for the suction surface)). The corresponding 
condensate arrival rate distribution (fig. 8 ) and liquid viscosity distributions are then 
inserted into the partial differential equation (PDE) governing the liquid layer thickness 
5 1 (x,z) (fig. 5). This nonlinear first order PDE is then numerically solved (by the method- 
of-characteristics) to provide the condensate layer streamline pattern (fig. 9), as well as 
the corresponding normalized liquid layer thickness 6 ^ (x, z)/8^ (0,0 ) (fig. 10). The 
absolute value of the root-nose liquid layer thickness, 5 ^ ( 0 , 0 + ), is, of course, environment- 
(e.g., alkali-level) dependent, and is given by a simple explicit formula (derived in ref. 11 ) 
involving -m"(0,0 + ) and (d‘ r w /3 x) x _q. Note that two byproducts of these calculations are a) 
the solvent "inventory" on the blade, and b) relative tip- and trailing edge salt runoff 
rates. The writer is indebted to his students, R. Nagarajan and S. Burke, for their 
assistance in generating these preliminary numerical results (figs. 6 - 10 ). 

4. METAL OXIDE "DISSOLUTION" RATE DISTRIBUTIONS 

Nonstoichiometric molten sodium sulfate is known to be a solvent for AI 2 O 3 (s) or C^O^s) 
— the oxides ordinarily relied upon for corrosion protection of the underlying blade alloys. 
Accordingly, the next step in this research program is prediction of the 
dissolution rate distributions associated with "solvent" flow patterns of the type described 
in Section 3. Indeed, it should be instructive to compare such predictions with hot corrosion 
rate patterns observed on blades removed from operational GT engines. Toward this end, 
estimates are first being made of the oxide diffusion coefficient in the solvent, equilibrium 
solubility, and maximum (kinetic) rate of dissolution. This information is then used to 
predict the solute diffusion-limited dissolution rate along each streamline (cf. fig. 9), 
using a generalization of Leveque-Levich BL-theory. Preliminary dissolution rate maps, now 
being obtained, will be checked, generalized, and discussed in terms of their parametric 
dependences, and agreement with operational experience. 

5, IMPLICATIONS, FUTURE WORK 

While much of the required input information remains to be tightened up, the present 
formalism is reaching the point where very instructive parametric studies can now be 
economically performed for both GT stator and rotor blades over a wide variety of 
environmental conditions (salt level, turbine inlet temperature, stagnation pressure level, 
blade contour and cooling, tip speed , etc.). Moreover, many aspects of the theory can now 
be generalized, including possibly mixed alkali deposits (ref. 7), partially "dry" and/or 
rough surfaces, condensate arrival by other (than vapor diffusion) mechanisms (refs. 9,10), 
and transient condensate flows. It may also be useful to examine the effects on localized 
dissolution rates cf secondary flows, produced in part by surface tension gradients. Feedback 
and guidance from the participants of this interdisciplinary HOST Workshop will be especially 
welcome. 
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DEPOSITION RATE (normalized) 


AIRFOIL DEPOSITION MODEL 


GRANT NAG 3-201: "THEORY OF MASS TRANSFER FROM COMBUSTION GASES," 

WITH PROFESSOR D.E. ROSNER, ChE DEPARTMENT, YALE UNIVERSITY 

EMPHASIS: TRACE SALT VAPOR DEPOSITION AND CORRESPONDING SHEAR- 

DRIVEN CONDENSATE LAYER FLOW 

OBJECTIVE: 

• OVERALL- DEVELOP MODEL TO PREDICT CORRODANT DEPOSITION ON 

TURBINE AIRFOILS 

• 1ST YEAR - MODEL DEPOSITION RATE FOR SEVERAL SIMPLE GEOMETRIES 

• 2ND YEAR - PREDICT AND DISPLAY LIQUID LAYER EVOLUTION ON TURBINE 

VANES AS A RESULT OF VAPOR DEPOSITION AND LIQUID LAYER FLOW 

Figure i 



TARGET TEMPERATURE 


Figure 2 
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Figure 5 



Figure 6 
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TEMPERATURE PROFILE ON SUCTION SURFACE 



Figure 7 



Figure 8 
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Figure 9 


LIQUID LAYER THICKNESS PROFILE ON SUCTION SURFACE 
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DEPOSITION MODEL VERIFICATION 1 


Suleyman A. Gokoglu 
Analex Corporation 
Cleveland* Ohio 

This research is motivated by gas turbine failures associated with deposition, 
erosion, and accelerated "hot" corrosion. A comprehensive yet tractable theoretical 
framework of deposition from combustion gases has recently been developed covering 
the spectrum of various mass delivery mechanisms including vapor (refs. 1 and 2), 
thermophoretically-enhanced small particle (ref. 3), and inertially impacting 
large particle (ref. 4) deposition. Rational yet simple correlations have been 
provided to facilitate engineering surface arrival rate predictions (refs. 5 and 
6) . The goal of the program at NASA Lewis Research Center is to experimentally 
verify the deposition theory. Toward this end, a Mach 0.3 burner rig has been 
designed to measure deposition rates from salt-seeded combustion gases on an 
internally-cooled cylindrical collector (ref. 7) . 

Although very good agreement is obtained between the NaCl -seeded experiments 
and the chemically-frozen boundary layer (CFBL) vapor deposition theory, consis- 
tently higher values have been measured for sea-salt, NaNOg, and Na2SC>4-seeded 
experiments. Vaporization time calculations of the sprayed salt solution drop- 
lets during their residence time inside the combustor have shown the possibility 
of the presence of particles capable of inertially impacting the collector. In 
fact, deposition rate calculations based on the mean droplet size have agreed 
well with the experimental observations. The disagreement between the theory 
and the experiments above the melting point of the deposit is attributed to the 
shear-driven molten deposit layer flow and run-off from the smooth collector sur- 
face (ref. 8). 

Future experiments will focus on the determination of parameters that will 
assure deposition by vapor diffusion only. Under these conditions, the deposition 
rate predictions of the CFBL theory will be verified using different salts 
(e.g., K2 so 4> Na 2 S0 4> NaCl, etc.). Above the melting point of the deposit, the 
theoretical molten layer thickness distribution prediction will be experimentally 
confirmed on a stationary cylindrical collector in crossflow. The experimental 
verification of the Mach number and the pressure dependence of the CFBL vapor 
deposition theory is also planned for the future. 
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1 Work done under NASA contract NAS3-23293 at the Lewis Research Center. 
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DROPLET SIZE DISTRIBUTION 
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PERTINENT TEST PARAMETERS 


• SALT INJECTION SYSTEM: 

SALT COMPOSITION 
SALT CONCENTRATION 
SALT SOLUTION PUMP SPEED 
ATOMIZING AIR PRESSURE (FLOWRATE) 
RESIDENCE TIME IN THE COMBUSTOR 
SALT VAPOR PRESSURE 

® BURNER: 

MASS FLOW RATE 

AIR TEMPERATURE INTO THE COMBUSTOR 
COMBUSTION GAS TEMPERATURE 
FUEL-TO-AIR RATIO 
SULFUR CONCENTRATION 
COMBUSTION GAS PRESSURE 
COMBUSTOR LINER DESIGN 
EXIT NOZZLE DESIGN 

• COLLECTOR: 

GEOMETRY 

TEMPERATURE AND TEMPERATURE GRADIENTS 
AT BETWEEN THE COLLECTOR AND GAS STREAM 
SURFACE ROUGHNESS 
ROTATING - NOT ROTATING 
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FUTURE WORK 


« DETERMINATION OF DROPLET SIZE DISTRIBUTION AS A FUNCTION OF: 

(a) SALT SOLUTION CONCENTRATION 

(b) SALT SOLUTION FLOWRATE 

(c) ATOMIZING AIR PRESSURE (FLOWRATE) 

(d) MACH NUMBER 

AT CONSTANT Na (OR K) ELEMENT CONCENTRATIONS IN BURNER 
FOR DIFFERENT SALTS (Na^NaCI, K^O^ ETC.) 

• DETERMINATION OF PARAMETERS THAT ASSURE DEPOSITION BY 
VAPOR DIFFUSION ONLY 

• DETERMINATION OF DISCHARGE COEFFICIENTS OF DIFFERENT NOZZLES 
USED IN DEPOSITION EXPERIMENTS 

• VERIFICATION OF ENTIRE DEPOSITION RATE CURVE (VIA CFBL THEORY) 
USING K 2 SO 4 SEEDED EXPERIMENTS 

• VERIFICATION OF CFBL PREDICTED PLATEAU DEPOSITION RATES USING 
VARIOUS Na-SALT-SEEDED EXPERIMENTS 

• VERIFICATION OF MOLTEN Na2S0 4 DEPOSIT THICKNESS DISTRIBUTION 
AROUND STATIONARY CYLINDER IN CROSSFLOW 

• VERIFICATION OF PRESSURE AND MACH NUMBER DEPENDENCE OF 
DEPOSITION THEORY 
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EFFECTS OF SURFACE CHEMISTRY ON HOT CORROSION LIFE* 

R, E, Fryxell 
General Electric Company 

This new program has as its primary objective the development of hot corrosion 
life prediction methodology based on a combination of laboratory test data and evalu- 
ation of field service turbine components which show evidence of hot corrosion. This 
program is divided into five tasks. All burner rig testing will be performed by TRW. 
The discussion will describe the overall program with a brief summary of activity to 
date. 


Task I involves the evaluation of six hot corroded components, having known 
operating history, to establish the degradation mechanisms in the corroded areas. 
Correlation of the corrosion with operating conditions will be made to the extent 
possible. Evaluation will include chemical and X-ray diffraction analyses of surface 
scales/deposits as well as metallography, scanning electron microscopy, and electron 
microprobe examination. Both uncoated and coated hardware will be evaluated and 
emphasis will be placed on the cause and effect parameters which are associated with 
hot corrosion. Task I will be completed by the end of this year. 

Task II will establish a hot corrosion baseline for the program alloys and 
coatings in the as-processed condition and will involve up to 1000 hours exposure of 
duplicate specimens in a Mach 0.3 burner rig using Jet A fuel (0.045-0.065% sulfur) 
and in which 0.5 ppm sodium by weight, as NaCl, is added to the combustion air. The 
test cycle will be one hour at 900C followed by at least six minutes of forced air 
cooling to ambient temperature. Tests are scheduled to begin in December 1983. 

At approximately 20 cycle intervals, specimens will be visually examined, photo- 
graphed, and inductance readings taken. Specimens will be removed when evidence of 
hot corrosion is noted in three successive test intervals. Unexposed specimens will 
then be inserted, as space is available, for time periods of 100, 300, and 500 hours 
if these time periods do not exceed 2/3 the time the original specimen of the same 
material was exposed. 

Specimens to be tested in this task include U700 and Rene" 80, both uncoated and 
with the following coatings: 

pack aluminide RT21 (Chromalloy) and Codep (General Electric) 
low pressure plasma NiCoCrAlY (Ni-23Co-18Cr-12Al-0. 3Y) 

Evaluation of tested specimens will be as in Task X. 

In Task III, triplicate specimens of the above coated alloys will be aged at 
1100C under a variety of conditions to determine effects of coating^alloy inter- 
diffusion and/or oxidation on hot corrosion behavior. 


*Contract NAS3-23926 
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These aging treatments will be isothermal inert atmosphere for 100 hours s iso- 
thermal air furnace oxidation for 100 , 300, 600 hours, one hour air furnace cycles 
for 100 hours and cyclic burner rig oxidation for 100, 300 and 600 hours . 

Weight changes will be determined on a regular basis (every 20 hours for cyclic 
exposures and at the end of exposure intervals for other cases). Also, specimen in- 
ductance shall be determined on the same basis. 

One specimen representing each of the above conditions will be evaluated using 
the methods outlined in Task I; the remaining duplicate specimens representing each 
condition shall be used in Task IV which is a hot corrosion burner rig test identical 
with Task II, 

At the conclusion of the Task IV test, all the data generated in Tasks I through 
IV will be evaluated and an empirical hot corrosion life prediction model based on 
these data will be proposed. Also, recommendations will be made for other test 
parameters to be evaluated, and evaluation methodology necessary to permit prediction 
of hot corrosion life. 

Task V will be a hot corrosion high velocity cyclic burner rig test designed to 
check the validity of the proposed life prediction model. Duplicate specimens of a 
maximum of four alloys and five coatings will be tested up to a maximum of 1000 hours 
and evaluated as in Tasks II and IV, The alloys and coatings will be those from this 
program plus additional alloy-coating systems selected by the NASA Project Manager. 

The Task V experiment will be designed to test two aspects of the life prediction 
model developed in Task IV. The ability of the model to predict hot corrosion life 
under different rig test conditions than those previously used in Task II and IV will 
be tested in one portion of the experiment. The second portion of the experiment will 
test the ability of the model to predict hot corrosion life of new alloys and coatings. 

The ability of the model to predict hot corrosion life under different rig test 

conditions will be evaluated using alloys and coatings previously tested in Task II 
and Task IV, The hot corrosion tests for this portion of the program will use one or 
more variations of the test parameters used in Task II, III and IV, Test parameter 
variations that will be considered include the following: 

o A different salt level 
o A different test temperature 
o Intermittent salt injection 
o A different sulfur level in the fuel 
o Additional or modified aging cycles 

The specific testing parameters will be derived from the results of Task II and 
Task IV, previous experience in evaluating field service hardware and previous experi- 
ence in conducting hot corrosion testing under a variety of test conditions. 

The ability of the model to predict hot corrosion resistance of new coatings and 
alloys will be evaluated using alloys and coatings that have not been evaluated pre- 
viously in Task II and Task IV, Rig test conditions will be identical to those used 
in Task II and Task IV, The aging cycle will be chosen from the results of Task IV. 

Life prediction for both portions of the Task V experiment will be based on the 
model developed in Task IV as well as extensive experience gained from evaluating 
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field service components and conducting hot corrosion tests. Knowledge of the effect 
of variations in both engine operating conditions and rig testing conditions will be 
used in predicting hot corrosion life under different rig test conditions in the first 
portion of the experiment with the same alloys and coatings tested in Task II and 
Task IV. Typical variations in the relative performance of four Ni base superalloys 
as a function of salt flux in a one atmosphere burner rig are shown in Figure 1. With 
respect to testing of additional alloys and coatings, use will be made of extensive 
empirical correlations of composition and hot corrosion performance in predicting hot 
corrosion life. 


Current activity on this urogram is primarily in Task X. Six comnements have 
been selected, and include the following: 

4 stage one high pressure turbine blades from CF6-50 engines, representing 
four different carriers. Codep coated Rene* 80. 

1 stage one low pressure turbine nozzle from a CF6-50 engine. Uncoated Rene’ 77, 


1 stage one high pressure turbine blade from a J79 engine. 
Evaluation is in progress and preliminary results will be described, 



Sea Salt, PPM 


FIGURE 1 
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Coating Life Prediction 


Michael A. Gedwili 

National Aeronautics and Space Administration 
Lewis Research Center 

At elevated temperatures* high strength turbine alloys are life 
limited because of severe oxidation attack. To improve life, the hot 
section components of gas turbine engines are coated to take advantage of 
the oxidation resistance of metallic coatings. The oxidation resistance 
of these coatings is derived from the formation of a protective oxide 
(generally AI 2 O 3 ) scale and reformation of the scale if it spalls during 
temperature excursions. During thermal cycling, reformation of the scale 
is very important since all scales eventually experience spallation; some 
significantly more than others. There are some scales (for example O 2 O 3 ) 
that volitalize rapidly in a dynamic (high gas velocity) environment at 
elevated temperatures. Thus, the coating serves as a reservoir for the 
formation/reformation of the protective scale. Interactions of coating 
and substrate, however, reduce the tendency for reformation of the protec- 
tive scale. Although vast improvements in life have been made and will be 
made in the future by coating turbine components, it should be obvious 
that coated components in themselves are life limited at high temperatures 
and during thermal cycling. 

The engine designer/manufacturer needs reliable, accurate and precise 
information to properly design turbine engines to coating life limits. 

Such engines would be cost effective because they would be capable of 
operating either for longer times or at higher temperatures. The end-user 
would also benefit from less costly engine inspections and component 
replacements which now are many times made on a conservative basis. A 
methodology for predicting life of coated turbine airfoils would, there- 
fore, benefit al 1 . 

The study to be discussed herein is a first-cut at an approach 
directed toward developing the life prediction methodology for metallic 
coatings on turbine alloys in cyclic oxidation. A cyclic oxida- 
tion/spalling computer model was developed at LeRC which predicts cyclic 
furnace oxidation behavior of bulk cast/wrought alloys over a range of 
cycling frequencies. The input data for the model are obtained from 
short-time isothermal oxidation tests. The required data are measurements 
of the oxide growth rate(s) and a spall constant, Q 0 , The model predicts 
specific weight changes, specific metal consumption, oxide growth rate(s), 
oxide spall fractions, and an effective time. The agreement between the 
observed and predicted cyclic oxidation behavior of bulk alloys has been 
shown to be within ±10%. In the present study, the computer model is 
being applied to an aluminide coating on 11-700, a low pressure plasma 
sprayed (LPPS) NiCoCrAl Y coating on U-700, and a bulk LPPS NiCoCrAl Y. To 
predict coating life in cyclic furnace oxidation, an empirical diffusion 
model to account for coating degradation (oxidation and coating/substrate 
interdiffusion) will be integrated with the oxidation/spall ing model. The 
integrated model will then be verified/adjusted to predict cyclic burner 
oxidation. Further verification/adjustment will lead to life prediction 
model for coated turbine airfoils. Results of isothermal and cyclic 
furnace oxidation of aluminide coated U-700, in particular, will be 
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presented and discussed. The results will also include morphological 
changes observed during furnace oxidation. 
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COATING LIFE PREDICTION 


OBJECTIVE: 

DEVELOP IMPROVED METHODOLOGY FOR PREDICTING CYCLIC 
OXIDATION LIFE OF METALLIC COATINGS ON GAS TURBINE 
AIRFOILS 


V-1721 


ENVIRQMENTAL AND SUBSTRATE REACTIONS DEGRADE COATINGS 



LPPS NiCoCrAlY COATED 
U-700 ANNEALED 
1080P C/4 hr/Ar 



REACTION 


- - - REACTION 


AFTER 100-1 hr CYCLES 
AT HOOP C 


V-1728 
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ISOTHERMAL VS CYCLIC OXIDATION 


SAMPLE 

SPECIFIC 

WEIGHT 

CHANGE 



ISOTHERMAL OR NONSCALE SPALLING 
SPECIFIC WEIGHT CHANGE CURVE 


TIME- 


NET SAMPLE WEIGHT CHANGE 

ACTUAL OXYGEN UPTAKE 



V-1723 


BASIC SPALL EQUATION 

w s ■ Qo (wi)» + 1 

WHERE: 

W $ = SPECIFIC WEIGHT OF SPALL ON COOLING FROM i th CYCLE 
Q 0 - SPECIFIC SPALL CONSTANT 

W r ' ■ SPECIFIC WEIGHT OF OXIDE PRIOR TO COOLING FROM i th 
CYCLE 

a = CONSTANT 


V ~ 17 19 
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CYCLIC FURNACE OXIDATION OF ALUMINIDE COATED U-700 


mm 


Wm~ mlii 

300 CYCLES 


OXIDE SCALE AT 100 CYCLES 


SPECIFIC 

WEIGHT 

CHANGE, 

mn/rm'- 


1000 CYCLES 


TO-10. 6 AT 1000 CYCLES 


NUMBER OF 1-hr CYCLES AT 1100° C 


300 CYCLES 


SPECIFIC 

WEIGHT 

CHANGE, 

mn/rm*- 


TO - 10.6 AT 
1000 CYCLES 


100 CYCLES 


0 


200 400 600 

NUMBER OF 1-hr CYCLES AT 1100° C 


800 





COATING: ALUMINIDE 
(TWO STEP) 


LPPS NiCoCrAiY 
(ANNEALED 4 hr/ 1080° C/Ar) 


\ 


OXIDATION OF ALUMINIDE COATED U-70G 

1150° C 













APPROACH 


NOMOLITHIC NiCoCrAlY 


NICoCrAIY/U-700 


ALUMINIDE/U-700 


SHORT TIME ISOTHERMAL 
FURNACE OXIDATION 
AT 1000° TO 1200° C 


PARAMETERS FOR 
CYCLIC oxidation; 
SPALLING MODEL 


EMPIRICAL 

DIFFUSION 

MODa 


CYCLIC FURNACE OXIDATION 
FOR MODEL VERIFICATION 


PRELIMINARY COATING 
LIFE PREDICTION MODEL 


VERIFICATION/AD JUSTMENT(S) 
OF MODEL FOR CYCLIC 
BURNER OXIDATION 


LIFE PREDICTION MODa 
FOR COATED TURBINE AIRFOILS 
















CYCLIC FURNACE OXIDATION OF ALUMINIDE COATED U-700 
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CONCLUDING REMARKS: HOST 2nd ANNUAL WORKSHOP 


Daniel E. Sokolowski 


National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, OH 44035 


ms reported at this workshop, the HOST Project activities are well underway 
and are starting to produce results. Contractor annual and final reports are 
starting to become available and their number should increase in the future. 
Workshops such as this will continue on an annual basis; the HOST 3rd Annual 
Workshop is tentatively scheduled for October 23-24 of next year. 

Figure 1 shows the objectives from the Aeronautics Long Range Plan of NASA. A 
review of the HOST activities, how they are being conducted, and their impact 
indicates that the HOST Project meets all five objectives of the plan. In 
addition, the HOST Project is now being recognized for the value of focused 
and interdependent research when compared with generic, independent base R&T 
activities. The nature of the problem being addressed, durability in this 
case, has much influence in advocating and successfully implementing such a 
project, however. 

To date, $11.5 million has been invested in HOST. In FY '84 another $6.1 
million will be spent. The present plan for FY '85- ‘86 is for $9.0 to 11.1 
million per year. Such funding will support the future directions listed in 
Figure 2. 

Finally, I want to say "Thank You" to the HOST Project Team for a job well 
done in conducting this workshop. In particular, I want to thank the 
contractor speakers, the Subproject Managers for being session chairmen, and 
my Assistant Manager, Bob Ensign and Administrative Assistant, Joanne Flowers 
for helping organize and coordinate the multitude of efforts required. 


265 



AERONAUTICS LONG RANGE PLAN OBJECTIVES 


1. ESTABLISH RECOGNITION OF IMPORTANCE OF NASA AERONAUTICS TO BOTH CIVIL 
AND MILITARY AVIATION 

2. PROVIDE U.S. R&T CAPABILITY BY MAINTAINING RESEARCH CENTERS IN 
POSITIONS OF UNDISPUTED EXCELLENCE IN FACILITIES, COMPUTATIONAL 
CAPABILITY AND TECHNICAL STAFF 

3. RESTORE A BALANCED AERONAUTICS PROGRAM CONSISTING OF DISCIPLINE 
RESEARCH, SYSTEMS RESEARCH AND SELECTED PROOF OF CONCEPT ACTIVITIES 
WITH EMPHASIS ON COMMON APPLICABILITY TO BOTH CIVIL AND MILITARY 
AVIATION 

4. STRENGTHEN NASA UNIVERSITY PARTNERSHIP IN AERONAUTICS R&T 

5. STRENGTHEN USER INTERFACES TO PROMOTE TECHNOLOGY TRANSFER 


FUTURE DIRECTIONS 

• CONTINUING HOT SECTION DURABILITY RESEARCH 

• SELECTED INTERDISCIPLINARY GRANTS 

• STRONG NASA LEWIS IN-HOUSE RESEARCH EFFORTS 
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